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IX  ABSTRACT  (Hartman 300  man*) 

SAP  IV  is  a  general  purpose  finite  element  analysis  program  for  linear  structural  analysis.  The  program  can  analyze 
structures  that  are  subject  to  either  static  or  dynamic  loads.  The  structures  can  be  described  in  terms  of  trass,  frame, 
plate,  shell,  or  brick  elements  or  combinations  of  elements.  The  static  loads  include  concentrated  loading  applied  to 
discrete  points  an  the  structural  models;  and  loads,  concentrated  loading  applied  to  discrete  points  on  the  structural 
models;  and  loads,  concentrated  or  uniform,  applied  to  individual  elements  in  the  model.  The  loads  can  be  combined 
into  load  cases.  The  dynamic  analysis  features  include  modal  analysis  to  obtain  mode  shapes  and  frequencies,  response 
history  analysis  using  either  modal  superposition  or  direct  integration  techniques,  and  response  spectrum  analysis. 
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PREFACE 


The  information  contained  in  this  report  is  technical  in  nature. 

It  does  not  contain. policy  of  the  Naval  Civil  Engineering  Laboratory, 
the  Naval  Facilities  Engineering  Command,  the  Assistant  Commander  of 
Engineering  and  Design,  or  the  Director  of  Engineering  Systems  Manage* 
roent  Division.  This  report  is  one  of  a  series  which  is  designed  to 
provide  technical  information  and  documentation  for  studies  conducted  in 
support  of  the  Naval  Facilities  Engineering  Command  Graphics  Engineering 
and  Mapping  Systems,  Computer-Aided  Design,  Cost  Engineering  System,  and 
Guide  Specifications  System  as  they  relate  to  the  Naval  Facilities  System 
and  life  cycle  facilities  management. 
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INTRODUCTION 


SAP  IV  is  a  general  purpose  finite  element  analysis  program  for 
linear  structural  analysis.  The  program  can  analyze  structures  that  are 
subject  to  either  static  or  dynamic  loads.  The  structures  can  be  described 
in  terms  of  truss,  frame,  plate,  shell,  or  brick  elements  or  combinations 
of  elements.  The  static  loads  include  concentrated  loading  applied  to 
discrete  points  on  the  structural  models;  and  loads,  concentrated  or 
uniform,  applied  to  individual  elements  in  the  model.  The  loads  can  be 
combined  into  load  cases.  The  dynamic  analysis  features  include  modal 
analysis  to  obtain  mode  shapes  and  frequencies,  response  history  analysis 
using  either  modal  superposition  or  direct  integration  techniques,  and 
response  spectrum  analysis. 


CAPABILITIES 


t 

GEOMETRIC  MODEL  1 

The  structure  to  be  analyzed  is  described  as  a  mathematical  model  of 
three  dimensional  coordinates.  These  coordinates  are  numerically 
identified  as  nodes.  Each  node  holds  a  discrete  location  in  three  dimen¬ 
sional  space.  The  node  can  be  free  to  move  with  translational  or 
rotational  displacement,  or  it  can  be  restrained  from  movement.'  The 
coordinate  system  can  be  either  a  cartesian  or  cylindrical  systems.  A 
limited  node  generation  feature  is  included. 


ELEMENT  LIBRARY 

There  are  nine  element  types  that  can  be  used  separately  or  combined 
to  complete  the  mathematical  description  of  the  structure  to  be  analyzed. 
Each  element  type  is  used  to  define  the  stiffness  properties  of  the  struc¬ 
ture.  The  element  types  include: 

Truss  Element 

The  truss  element  is  a  three  dimensional  element.  Axial  forces  and 
stress  are  calculated  for  each  member .  A  uniform  temperature  change  and 
inertia  loads  in  three  directions  can  be  considered  as  the  basic  element 
load  conditions. 

Beam  Element 

The  beam  element  is  a  three  dimensional  element.  Forces  (axial  and 
shear)  and  moments  (bending  and  torsion)  are  calculated  with  respect  to  the 
local  beam  coordinate  system  for  each  beam  element.  The  beams  are  assumed 
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to  be  symmetrical  about  and  loaded  in  the  plans  of  their  major  and  minor 
axes.  Provision  for  unsymroetrical  beams  are  not  currently  in  the  program, 
however,  this  feature  has  been  developed.  The  effect  of  shear  deformation 
and  member  end  releases  can  be  included  in  the  analysis.  The  element  loads 
include  concentrated  loads  (no  moments)  due  to  gravity,  and  fixed-end 
forces.  Member  uniform  and  concentrated  loads  and  thermal  loads  are  speci¬ 
fied  using  fixed-end  forces  (including  moments). 

Plane  Stress  Membrane  Elements 

Quadrilateral  and  triangular  elements  can  be  used  for  plane  stress 
membrane  elements  of  specified  thickness  which  are  oriented  in  an  arbitrary 
plane.  These  elements  have  temperature  dependent  orthotropic  material  pro¬ 
perties.  Incompatible  displacement  modes  can  be  included  to  improve  the 
bending  properties  of  the  element.  Gravity,  inertia,  pressure,  and  temper¬ 
ature  loads  are  included. 

Two-Dimensional  Finite  Elements 

Quadrilateral  and  triangular  elements  can  be  used  as:  1)  axisymmetric 
solid  elements  symmetrical  about  the  z-axis,  2)  plane  strain  elements  of 
unit  thickness  in  the  Y-Z  plane,  and  3)  plane  stress  elements  of  specified 
thickness  in  the  Y-2»  plane.  All  of  these  elements  have  temperature  depend¬ 
ent  orthotropic  material  properties.  Incompatible  displacement  modes  can 
be  included  to  improve  the  bending  properties  of  the  element.  Gravity 
inertia,  pressure,  and  temperature  loads  are  Included. 

Three-Dimensional  Solid  Elements 

General  three-dimensional,  eight-node,  isoparametric  elements  with 
three  translational  degrees  of  freedom  per  node  can  be  used.  The  elements 
use  isotropic  material  properties.  The  element  load  cases  include  surface 
pressure,  hydrostatic  loads,  inertia  loads,  and  thermal  loads.  The  six 
components  of  stress  and  the  three  principal  stresses  are  computed  at  the 
center  of  each  element.  Surface  stresses  are  also  evaluated.  No  incompat¬ 
ible  displacement  modes  are  assumed  in  the  formation  of  the  element 
stiffness. 

Plate  and  Shell  Elements 

A  quadrilateral  element  is  used  to  describe  the  plate  or  shell  element. 
Anisotropic  as  well  as  isotropic  material  properties  are  possible.  The 
element  load  cases  include  surface  pressure,  inertia  loads,  and  thermal 
loads.  For  the  analysis  of  shallow  shell,  rotational  constraints  normal  to 
the  surface  may  be  imposed  by  using  boundary  elements  at  the  nodes.  The 
results  include  moments  per  unit  length  and  membrane  stresses. 

Boundary  Element 

This  element  is  used  to  constrain  nodal  displacements  to  specified 
values,  to  compute  support  reactions,  and  to  provide  linear  elastic 
supports  to  nodes,  The  boundary  element  is  essentially  a  spring  which  has 
a  translational  and  rotational  stiffnesses. 
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General  Three-Dimensional  Isoparametric  Elements 

A  minimum  of  8  and  a  maximum  of  21  nodes  are  used  to  describe  a 
general  three  dimensional  isoparametric  element.  The  element  is  used  to 
represent  orthotropic  elastic  media.  The  use  of  nodes  9  through  21  is 
optional.  Element  loads  include:  surface  pressure,  hydrostatic  loads, 
inertia  loads,  and  thermal  loads.  Six  global  stresses  are  computed  at  up 
to  7  user  selected  locations  within  the  element. 

Three-Dimensional  Pipe  Elements 

Straight  or  curved  cylindrical  pipe  elements  may  be  used  to  describe 
pipe  networks.  Axial  and  shear  forces,  torque  and  bending  moments  are 
calculated  for  each  member.  Gravity  loadings  in  the  global  (X,  Y,  Z) 
directions,  uniform  temperature  changes,  and  extensional  effects  due  to 
internal  pressure  form  the  basic  member  loading  conditions. 


STATIC  ANALYSIS 

The  static  analysis  features  include  analysis  for  nodal,  element, 
gravity,  and  temperature  loads.  The  element  loads  have  been  described  with 
the  element  library  where  each  element  type  has  particular  element  loads 
that  can  be  considered.  The  nodal  loads  can  be  defined  for  each  of  tfce 
nodes.  These  loads  include  concentrated  forces  and  moments  for  each  of  the 
three  degrees  of  freedom.  The  nodal  loads  can  be  assigned  to  a  load  case. 
The  load  cases  are  completed  by  assigning  the  appropriate  element  load 
sets,  of  which  there  are  four.  The  element  loads  are  assigned  to  the  load 
case  using  a  multiplier  factor.  In  this  manner  any  number  of  loading 
conditions  can  be  studied  separately  or  combined. 


DYNAMIC  ANALYSIS 

There  are  four  types  of  dynamic  analysis  options  available.  The  modal 
analysis  option  computes  the  mode  shapes  and  associated  modal  frequencies. 
Response  history  analysis  or  time  history  analysis  can  be  accomplished  using 
two  methods.  The  first  uses  the  modal  superposition  technique  which 
includes  computing  the  mode  shapes  and  frequencies.  The  second  uses  the 
direct  integration  method.  The  last  dynamic  analysis  option  is  response 
spectrum  analysis  which  also  uses  the  mode  shapes  and  frequencies  to 
compute  the  R.M.S  (root -mean-square)  stresses  and  deflections  due  to  an 
input  displacement  of  acceleration  spectrum. 
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SOLUTION  METHOD 


The  numerical  methods  used  in  the  SAP  IV  program  are  well  documented 
in  many  text  books  dealing  with  the  finite  element  method.  The  primary 
reference  for  this  program  is  a  report  that  is  available  through  the 
Earthquake  Engineering  Research  Center,  College  of  Engineering,  University 
of  California,  Berkeley,  California  (Reference  1).  The  original  version  of 
thi.*  computer  program  was  also  obtained  from  the  Earthquake  Engineering 
Research  Center.  Additional  information  on  the  techniques  used  in  the  pro¬ 
gram  can  be  obtained  from  Reference  2.  Discussion  on  the  formulation  of 
the  element  stiffness  equations  similar  to  those  used  in  the  program  can  be 
obtained  in  Reference  3. 


* 
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PROBLEM  DESCRIPTION  DATA  PREPARATION  INSTRUCTIONS 


HEADING  LINE  (A72) 

Notes  Columns  Variables  Entry 

(1)  1-72  HED  Enter  the  heading  information  to 

be  printed  with  the  output 

NOTES: 

(1)  Begin  each  new  data  set  with  a  new  heading  line.  Several 
problems  can  be  run  at  the  same  time  by  adding  each  problem  data  set  one 
after  another.  To  stop  the  program  without  an  error  requires  two  blank 
lines  at  the  end  of  the  last  data  set. 


MASTER  CONTROL  LINE  (815) 


Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

NUMNP 

Total  number  of  nodal  points  (joints) 
in  the  model 

(2) 

6-10 

NELTYP 

Number  of  element  groups 

(3) 

11  -  15 

LL 

Number  of  structure  load  cases; 

6E.  1;  static  analysis 

EQ.  0;  dynamic  analysis  . 

i 

(A) 

16  -  20 

NF 

Number  of  frequencies  to  be  found 
in  the  eigenvalue  solution; 

EQ.  0;  Static  analysis 

GE.  1;  dynamic  analysis 

(5) 

21  -  25 

NDYN 

Analysis  type  code: 

EQ.  0;  static  analysis 

EQ.  1;  Eigenvalue/vector  solution 

EQ.  2;  forced  dynamic  response  by 
mode  superposition 

EQ.  3;  response  spectrum  analysis 

EQ.  A;  direct  step-by-step 
integration 

(6) 

26  -  30 

MODEX 

Program  execution  mode: 

EQ.  0;  problem  solution 

EQ.  1;  data  check  only 

(7) 

31  -  35 

NAD 

Total  number  of  vectors  to  be  used  in 

a  SUBSPACE  ITERATION  solution  for 
eigenvalues/vectors : 

EQ.  0;  default  set  to: 

MIN  [2*  NF,NF±41J 
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Notes 

Columns 

Variables 

Entrv 

(8) 

36  -  40 

KEQB 

Number  of  degrees  of  freedom 
(equations)  per  block  of  storage: 

EQ.  0;  calculated  automatically  by 
the  program 

(9) 

40  -  45 

N10SV 

Element  stress  porthole  flag 

EQ.  0;  no  element  stresses  saved 
GT.  0;  Save  element  stresses  on 
TAPE 10 

NOTES: 

(1) 

Nodes  are 

labeled  with 

integers  ranging  from  "l"  to  the  total 
_ ••  _ 

number  of  nodes  in  the  system,  "NUMNP" .  The  program  exits  with  no  diag¬ 
nostic  message  if  NUMNP  is  zero.  Thus,  two  blank  lines  are  used  to  end  the 
last  data  case  in  a  rtm;  i.e.,  one  blank  heading  line  and  one  blank  line 
for  the  master  control  line. 

(2)  For  each  different  element  type  (TRUSS,  BEAM,  etc)  a  new  element 
group  needs  to  be  defined.  Elements  within  groups  are  assigned  integer 
labels  ranging  from  "l"  to  the  total  number  of  elements  in  the  group. 

Element  groups  are  defined  in  the  section  dealing  with  elements,  below. 

Element  numbering  must  begin  with  one  in  each  different  group.  It  is 
possible  to  use  more  than  one  group  for  an  element  type.  For  example,  all 
columns  (vertical  beams)  of  a  building  may  be  considered  on  group  and  the 
girders  (horizontal  beams)  may  be  considered  another  group. 

(3)  J^t  least  one  load  condition  must  be  specified  for  a  static 
(NDYN  equals  0)  analysis.  If  the  data  case  calls  for  one  of  the  dynamic 
analysis  options  (NDYN  equals  1,  2,  3,  or  4)  no  load  cases  can  be  requested 
(i.e.,  LL  is  input  as  "0").  The  program  always  processes  the  Concentrated 
Load/Mass  Data  and  Element  Load  Multiplier  data  and  expects  to  read  some 
data.  For  the  case  of  a  dynamic  analysis  (NDYN.GE.  1)  only  mass  coefficients 
can  be  input  in  the  Concentrated  Load/Mass  Data  and  one  blank  element  load 
multiplier  data  line  is  expected. 

(4)  For  a  static  analysis,  NF  equals  0.  If  NDYN  equals  1,  2,  or  3, 
the  lowest  NF  eigenvalues  are  determined  by  the  program.  Note  that  a  dynamic 
solution  may  be  restarted  after  eigenvalue  extraction  (providing  a  previous 
eigenvalue  solution  for  the  model  was  saved  as  described  in  Dynamic  Analysis 
Restart  Instructions  section).  NF  for  the  original  and  restart  runs  must 

be  the  same. 

(5)  If  NDYN  equals  2  or  NDYN  equals  3  the  program  first  solves  for  NF 
eigenvalues/vectors  and  then  performs  the  forced  response  solution  or  the 
response  spectrum  analysis).  Thus,  the  program  expects  to  read  the  control 
line  governing  the  eigensolution  (mode  shapes  and  frequencies)  before  read¬ 
ing  data  in  either  the  Response  History  or  Response  Spectrum  Analysis  section. 
For  the  case  NDYN  equals  1,  the  program  solves  for  NF  eigenvalues/vectors, 
prints  the  results  and  proceeds  to  the  nex*.  data  case.  The  results  for  the 
eigen-value  solution  phase  (NDYN  equals  1)  may  be  saved  for  later  use  in 
automatic  restart  (refer  to  the  two  restart  sections  for  the  procedure 
required  to  affect  this  save  operation  for  a  subsequent  restart),  i.e.,  a 
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dynamic  solution  may  be  restarted  without  repairing  the  solution  for  modes 
and  frequencies.  If  this  data  case  is  a  restart  job,  set  NDYN  equals  -2  for 
forced  response  solution,  or  set  NDYN  equals  -3  for  a  response  spectrum 
analysis.  Note  that  the  solution  may  be  restarted  a  multiple  of  times  (to 
run  different  ground  spectra  or  different  timedependent  forcing  functions) 
because  the  program  does  not  destroy  the  contents  of  the  restart  file. 

If  NDYN  equals  4  the  program  performs  the  response  solution  by  direct 
step-by-step  integration  and  no  eigenvalue  solution  control  line  should  be 
provided. 

(6)  In  the  data-check-only  node  (MODEX  equals  1),  the  program  writes  only 
one  file,  "TAPE8",  and  this  file  may  be  saved  for  use  as  input  to  special 
purpose  programs  such  as  mesh  plotters,  etc.  TAPE8  contains  all  data  input 

in  its  completely  generated  form.  If  MODEX  equals  1,  most  of  the  expensive 
calculations  required  during  normal  (MODEX  equals  0)  execution  are  passed. 
TAPE8,  however,  is  also  written  during  normal  problem  solution. 

Note  that  a  negative  value  for  NDYN  ("-2"  or  "-3"),  when  executing  in 
the  data-check-only  mode,  does  not  cause  the  program  to  read  the  restart 
file  which  contains  the  eigensolution  information;  instead,  the  program 
jumps  directly  from  this  line  to  the  Response  History  or  Response  Spectrum 
section  and  continues  reading  and  checking  data  lines  without  performing 
the  solution. 

(7)  If  the  program  is  to  solve  eigenvalues  using  the  subspace 
iteration  algorithm,  the  entry  in  column  31-35  can  be  used  to  change  the 
total  number  of  Iteration  vectors  to  be  used  from  the  default  minimum  of 
2*NF  or  NF+8  (whichever  is  smaller)  to  the  value  "NAD."  The  effect  of 
increasing  NAD  over  the  default  value  is  to  accelerate  convergence  in  the 
calculations  for  the  lowest  NF  eigenvalues.  NAD  is  principally  a  program 
testing  parameter  and  should  normally  be  left  blank. 

(8)  KEQB  is  a  program  testing  parameter  that  allows  the  user  to  test 
multiple  equation  block  solutions  using  small  data  cases  that  would  other¬ 
wise  be  a  one  block  problem.  KEQB  is  also  used  to  reduce  the  I/O  files  in 
the  solution  procedure.  During  the  solution  complete  equation  blocks 
(i.e.,  bandwidth  times  the  number  of  equation  in  a  block)  are  written  and 
read  from  slow  speed  storage.  If  there  are  limitations  on  the  size  that 
can  be  transferred  then  KEQB  can  be  used  to  reduce  the  block  size.  KEQB  is 
normally  left  blank. 

(9)  N10SV  is  a  flag  to  determine  whether  element  stresses  are  to  be 
written  in  binary  form  onto  TAPE10.  This  porthole  file  can  then  be  used 
with  a  post  processor. 


7 


NODAL  POINT  DATA  (Al, 

I4,6I5,3F10. 

.  0, 15,F10. 0) 

Notes 

Columns 

Variables 

Entry 

(1) 

1 

CT 

Symbol  describing  coordinate  system 
for  this  node; 

EQ.  ;  (blank)  cartesian  (X,Y,Z) 

EQ.  C;  cylindrical  (R,Y,0) 

(2) 

2-5 

N 

Node  number 

(3) 

6  -  10 

IX(N, 1) 

X-translation  boundary  condition  code 

11  -  15 

IX(N, 2) 

Y-translation  boundary  condition  code 

16  -  20 

IX(N,3) 

Z-translation  boundary  condition  code 

21  -  25 

IX(N,4) 

X- rotation  boundary  condition  code 

26  -  30 

IX(N,5) 

Y-rotation  boundary  condition  code 

31  -  35 

IX(N,6) 

Z- rotation  boundary  condition  code 

EQ.  0;  free  (loads  allowed) 

EQ.  1;  fixed  (no  load  allowed) 

GT.  1;  master  node  number  (beam  nodes 

only) 

(A) 

36  -  45 

X(N) 

X  (or  R)  -ordinate 

46  -  55 

Y(N) 

Y  -ordinate 

56  -  65 

Z(N) 

Z  (or  8)  -ordinate  (degrees) 

(5) 

66  -  70 

KN 

Node  number  increment 

(6) 

NOTES: 

71  -  80 

T(N) 

Nodal  temperature 

(1)  A  special  cylindrical  coordinate  system  is  allowed  for  the  global 
description  of  nodal  point  locations.  If  a  "c"  is  entered  in  line  column 
1,  then  the  entries  given  in  columns  36-65  are  taken  to  be  references  to  a 
global  (R,Y,0)  system  rather  than  to  the  standard  (X,Y,Z)  system.  The 
program  converts  cylindrical  coordinate  reference  to  cartesian  coordinates 
using  the  formulas: 

X  -  R  sin  0 

Y  *  Y 

Z  =  R  cos  6 

Cylindrical  coordinate  input  is  merely  a  user  convenience  for  locating 
nodes  in  the  standard  (X, Y, Z)  system,  and  no  other  references  to  the 
cylindrical  system  are  implied;  i.e.,  boundary  condition  specifications, 
output  displacement  components,  etc.  are  referenced  to  the  (X,Y,Z)  system. 

(2)  Nodal  point  data  must  be  defined  for  all  (NUMNP)  nodes.  Node 
data  may  be  input  directly  (i.e.,  each  node  on  its  own  individual  line), 
may  be  entered,  or  the  generation  option  may  be  used  if  applicable  (see 
note  5,  below). 
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Admissible  nodal  point  numbers  range  from  "l"  to  the  total  number  of 
nodes  "NUMNP".  Illegal  nodal  point  numbers  are:  N.S.O  or  N.>.NUMNP. 

(3)  Boundary  condition  codes  can  only  be  assigned  the  following 
values  (M  =  1,2, — ,6) : 

IX(N,M)  =  0;  unspecified  (free)  displacement  (or  rotation) 
component 

IX(N,M)  =  1;  deleted  (fixed)  displacement  (or  rotation) 
component 

IX(N,M)  =  K;  node  number  "K"  (1  <  K  £  NUMNP  and  K  *  N)  is  the 

"master"  node  to  which  the  Mth  degree  of  freedom  at 
node  "N"  is  a  "slave" 

An  unspecified  (IX(N,M)  =  0)  degree  of  freedom  is  free  to  translate  or 
rotate  as  the  solution  dictates.  Concentrated  forces  (or  moments)  may  be 
applied  (see  instructions  pertaining  to  concentrated  load/mass  data)  in 
this  degree  of  freedom.  One  system  equilibrium  equation  is  required  for 
each  unspecified  degree  of  freedom  in  the  model.  The  maximum  number  of 
equilibrium  equations  is  always  less  than  six  times  the  total  number  of 
nodes  in  the  model. 

Deleted  (IX(N,M)  =  1)  degrees  of  freedom  are  removed  from  the  final 
set  of  equilibrium  equations.  Deleted  degrees  of  freedom  are  fixed  (points 
of  reaction),  and  loads  applied  to  these  degrees  of  freedom  are  ignored  by 
the  program.  Nodes  that  are  used  for  geometric  reference  only  (i.e.,  nodes 
not  assigned  to  any  element)  must  have  all  six  degrees  of  freedom  deleted. 
Nodal  degrees  of  freedom  having  undefined  stiffness  (such  as  rotations  in 
an  all  TRUSS  model,  out -of -plane  components  in  a  two-dimensional  planar 
model,  etc.)  should  be  deleted.  Deletions  have  the  beneficial  effect  of 
reducing  the  size  of  the  set  of  equations  that  must  be  solved.  The  table 
below  lists  the  types  of  degrees  of  freedom  that  are  defined  by  each  differ¬ 
ent  element  type.  The  table  was  prepared  assuming  that  the  element  has 
general  orientation  in  (X,Y,Z)  space. 

DEGREES  OF  FREEDOM  WITH  DEFINED  STIFFNESS 


ELEMENT  TYPE 

6X 

6Y 

6Z 

S9X 

«0¥ 

6ez 

1. 

TRUSS 

X 

X 

X 

2. 

BEAM 

X 

X 

X 

X 

X 

X 

3. 

MEMBRANE 

X 

X 

X 

4. 

2D/ QUADRILATERAL 

X 

X 

5. 

3D/BRICK 

X 

X 

X 

6. 

PLATE/SHELL 

X 

X 

X 

X 

X 

X 

7. 

BOUNDARY 

X 

X 

X 

X 

X 

X 

8. 

THICK  SHELL 

X 

X 

X 

9. 

3D/PIPE 

X 

X 

X 

X 

X 

X 
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Hence,  for  an  all  3D/BRICK  model,  only  the  X,Y,Z  translations  are 
defined  at  the  node,  and  the  number  of  equations  can  be  cut  In  half  be 
deleting  the  three  rotational  components  at  every  node.  If  a  node  is 
common  to  two  or  more  different  element  types,  then  the  nontrivial  degrees 
of  freedom  are  found  by  combination.  For  example,  all  six  components  are 
possible  at  a  node  common  to  both  BEAM  and  TRUSS  elements;  i.e.,  the  BEAM 
governs . 

A  "master/slave"  option  is  allowed  to  model  rigid  links  in  the  system. 
For  this  case,  IX(N,M)  =  K  means  that  the  Mth  degree  of  freedom  at  node  "N" 
is  "slave"  to  (dependent  on)  the  same  (Mth)  degree  of  freedom  at  node  "K"; 
node  "K"  is  said  to  be  the  master  node  to  which  node  N  is  slave.  Note  that 
no  actual  beam  need  to  run  from  node  X  to  node  N,  however  the  following 
restrictions  hold: 

(a)  Node  number  1  cannot  be  a  master  node;  i.e.,  K  i  1. 

(b)  Nodes  "N"  and  "K"  must  be  beam-only  nodes;  i.e.,  no  other 
element  type  may  be  connected  to  either  node  N  or  K. 

(c)  A  node  "N"  can  be  slave  to  only  one  master  node,  "K"; 
multiple  nodes,  however,  can  be  slave  to  the  same  master. 

(d)  If  the  beam  from  "N"  to  "K"  is  to  be  a  rigid  link 
arbitrarily  oriented  in  the  X,Y,Z  space,  then  all  six 
degrees  of  freedom  at  node  "N"  must  be  made  slaves  to 
node  "K". 

Displacement/rotation  components  for  slave  degrees  of  freedom  at 
node  "N"  a  not  recovered  for  printing;  i.e,  zeroes  appear  as  output  for 
8 lave  degrees  of  freedom. 

(4)  When  CT  (Column  1)  is  equal  to  the  character  "C”,  the  values 
input  in  Columns  36-65  are  interpreted  as  the  cylindrical  (R,Y,0) 
coordinates  of  node  "N".  Y  is  the  axis  of  symmetry.  R  is  the  distance  of 
a  point  from  the  Y-axis.  The  angle  6  is  measured  clockwise  from  the 
positive  Z-axis  when  looking  in  the  positive  Y  direction.  The  cylindrical 
coordinate  values  are  printed  as  entered  on  the  line,  but  imsiediately  after 
printing  the  global  cartesian  values  are  computed  from  the  input  entries. 
Note  that  boundary  condition  codes  always  refer  to  the  (X,Y,Z)  system  even 
if  the  node  happens  to  be  located  with  cylindrical  coordinates. 

(5)  Nodal  Point  lines  need  not  be  input  in  node-order  sequence; 
eventually,  however,  all  nodes  in  the  integer  set  {1,  NUMNP)  must  be 
defined.  Joint  data  for  a  series  of  nodes, 

JNj,  Nj  +  1  x  KN2,  Nj  +  2  x  KN2,...,N2) 
may  be  generated  from  information  given  on  two  data  lines  in  sequence: 

LINE  1  {N1,IX(N1,1),...,IX(N1,6),X(N1),...,KN1,T(N1)) 

LINE  2  {N2,IX(N2,1),...,IX(N2,6),X(N2),...,KN2,T(N2)) 
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KN2  is  the  mesh  generation  parameter  given  on  the  second  line  of  a  sequence 
The  first  generated  node  is  N^+  (lx  KN2);  the  second  generated  node  is  N  + 
(2  x  KN2),  etc.  Generation  continues  until  node  number  N_-  KN  is  estab¬ 
lished.  Note  that  the  node  difference  N2~  N.  must  be  evenly  divisible  by 
KN2.  Intermediate  nodes  between  and  are  located  at  equal  intervals 
along  the  straight  line  between  the  two  points.  Boundary  condition  code 
for  the  generated  data  are  set  equal  to  the  values  given  on  the  first  line. 
Node  temperatures  are  found  by  linear  interpolation  between  T(N.)  and  T(N  ) 
Coordinate  generation  is  always  performed  in  the  (X,Y,Z)  system}  and  no  2 
generation  is  performed  if  KN2  is  zero  (blank). 

(6)  Nodal  temperatures  describe  the  actual  (physical)  temperature 
distribution  in  the  structure.  Average  element  temperatures  established 
from  the  nodal  values  are  used  to  select  material  properties  and  to  compute 
thermal  strains  in  the  model  (static  analysis  only). 
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THREE-DIMENSIONAL  TRUSS  ELEMENTS,  ELEMENT  TYPE  1 

Truss  elements  are  identified  by  the  number  1.  Axial  forces  and 
stresses  are  calculated  for  each  member.  A  uniform  temperature  change  and 
inertia  loads  in  three  directions  can  be  considered  as  the  basic  element 
load  conditions.  The  truss  elements  are  described  by  the  following 
sequence  of  lines: 

A.  Control  Line  (315) 

Notes  Columns  Variables  Entry 

1  -  5  The  number  1 

6-10  Total  number  of  truss  elements 

11-15  Number  of  material  property  lines 

B.  Material  Property  Lines  (I5,5F10.0) 

There  need  be  as  many  of  the  following  lines  as  are  necessary  to 
define  the  properties  listed  below  for  each  element  in  the  structure. 

Notes  Columns  Variables  Entry 

1-5  Material  identification  number' 


6-15 

Modulus  of  elasticity 

16  -  25 

Coefficient  of  thermal  expansion 

26  -  35 

Mass  density  (used  to  calculate  mass 
matrix) 

36  -  45 

Cross -sectional  area 

46  -  55 

Weight  density  (used  to  calculate 

gravity  loads) 


C.  Element  Load  Factors  (4F10.0)  Four  lines 

Three  lines  specifying  the  fraction  of  gravity  (in  each  of  the  three 
global  coordinate  directions)  to  be  added  to  each  element  load  case. 
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Line  1:  Multiplier  of  gravity  load  in  the  +X  direction 


Notes 

Columns 

Variables  Entry 

1  -  10 

Element  load  case  A 

11  -  20 

Element  load  case  B 

21  -  30 

Element  load  case  C 

31  -  40 

Element  load  case  D 

Line  2: 

As  above  for  gravity  in  the  +Y  direction 

Line  3: 

As  above  for  gravity  in  the  +Z  direction 

Line  4: 

As  above  to  indicate  the  fraction  of  the  thermal 
load  to  be  added  to  each  of  the  eleaent  load  cases. 

D. 

Trass  KI< 

»ent  Data  Lines  (4I5,F10.0,I5) 

One 

line  per 

element  in  increasing  numerical  order  starting  with  one. 

Notes 

Columns 

Variables  Entrv 

1  -  5 

Eleaent  nuaber 

6-10 

Node  nuaber  I 

11  -  15 

Node  nuaber  J 

\ 

16  -  20 

Material  property  nuaber 

(2) 

21  -  30 

Reference  teaperature  for  zero  stress 

(1) 

NOTES: 

31  -  35 

Optional  parameter  k  used  for  autoaatic 
generation  of  eleaent  data. 

(1)  If  a  series  of  elements  exist  such  that  the  eleaent  nuaber,  N., 
is  one  greater  than  the  previous  eleaent  nuaber  (i.e.,  +  1)  and 
the  nodal  point  nuaber  can  be  given  by 


*i  *  Ii-i  ♦  k 


Ji  «  Ji-1  +  k 


then  only  the  first  eleaent  in  the  series  need  be  provided.  The  eleaent 
identification  nuaber  and  the  teaperature  for  the  generated  eleaents  are 
set  equal  to  the  values  on  the  first  line.  If  k  (given  on  the  first  line) 
is  input  as  zero  it  is  set  to  1  by  the  program. 


13 


(2)  The  element  temperature  increase  AT  used  to  calculate  thermal 
loads  is  given  by 

T  =  (Tf  +  Tj)/2.0  -  Tr 

Where  (T,  +  T.)/2.0  -  is  the  average  of  the  nodal  temperatures 
specified  on  the  nodal  point  data  lines  for  nodes  i  and  j;  and  T  is  the 
zero  stress  reference  temperature  specified  on  the  element  line.  For  truss 
elements  it  is  generally  more  convenient  to  set  T.  =  T.  =  0.0  such  that 
AT  =  -T  (note  the  minus  sign).  ^ 

Otfier  types  of  member  loadings  can  be  specified  using  an  equivalent 
AT: 

(a)  If  a  truss  member  has  an  initial  lack  of  fit  by  an  amount  d 
(positive  if  to  long)  then  AT  =  d/(oL). 

(b)  If  an  initial  prestress  force  P  (positive  if  tensile)  is  applied 
to  the  member  ends  that  is  released  after  the  member  is  connected  to  the 
rest  of  the  structure  then  AT  =  -P/(aAE). 

In  the  above  formulas  A  =  cross  section  area,  L  =  member  length  and  a  = 
coefficient  of  thermal  expansion. 


THREE-DIMENSIONAL  BE  JI  ELEMENTS,  ELEMENT  TYPE  2 


Beam  elements  are  identified  by  the  number  2.  Force  (axial  and  shear) 
and  moments  (bending  and  torsion)  are  calculated  (in  the  beam  local  coordi¬ 
nate  system)  for  each  beam.  Gravity  loadings  in  each  coordinate  direction 
and  specified  fixed  and  forces  form  the  basic  element  load  conditions. 

The  beam  elements  are  described  by  the  following  sequence  of  lines: 


A. 

Control  Line  (515) 

Notes 

Columns  Variables 

Entrv 

1-5 

The  number  2 

6-10 

Total  number  of  beam  elements 

11  -  15 

Number  of  element  property  lines 

16  -  20 

Number  of  fixed  end  force  sets 

21  -  25 

Number  of  material  property  lines 

B. 

Material  Property  Lines 

(I5.3F10.0) 

Notes 

Columns  Variables 

Entry 

1-5 

Material  identification  number 

6-15 

Young's  modulus 

16  -  25 

Poison's  ratio 

26  -  35 

Mass  density  (used  to  calculate  mass 
matrix) 

36  -  45 

Weight  density  (used  to  calculate 
gravity  loads) 

C. 

Element  Property  Lines 

(15.6F10.0) 

Notes 

Columns  Variables 

Entry 

1  -  5 

Geometric  property  number 

6  -  15 

Axial  area 

16  -  25 

Shear  area  associated  with  shear  forces 
in  local  2-direction 

15 


Notes 


Columns 

Variables 

Entry 

26  -  35 

Shear  area  associated  with  shear  forces 
in  local  3-direction 

36  -  45 

Torsional  inertia 

46  -  55 

Flexural  inertia  about 

local  2-axis 

56-65 

Flexural  inertia  about 

local  3-axis 

(1)  One  line  is  required  for  each  unique  set  of  properties.  Shear 
areas  need  be  specified  only  on  short  beams  where  shear  deformations  are  on 
the  same  order  of  magnitude  as  the  bending  deformations  and  are  to  be 
included  in  the  analysis. 


(2)  Section  property  computations  should  be  done  carefully. 


(a)  The  axial  area  is  the  total  cross  section  area. 

(b)  The  shear  area  represent  the  beam  cross  sectional  area 
effective  in  shear.  This  is  used  to  account  for  shear 
deformations  due  to  shear  strains.  These  deformations 
are  usually  neglected,  however,  they  may  be  important 
in  deep  beams.  The  shear  area  is  defined  as  A  =  K  * 

A  where  K  is  the  shear  coefficient,  a  factor  depending 
upon  the  ihape  of  the  cross  section.  The  factor  is 
always  less  than  one.  Reference  4  provides  more 
information  on  this  subject  and  it  includes  formulas 
for  the  shear  coefficient  for  various  common  cross 
sections. 

(c)  The  torsional  inertia  is  dependent  upon  the  form  and 
shape  of  cross  section.  For  a  circular  cross  section 
the  torsional  inertia  is  the  polar  moment  of  inertia; 
for  other  sections  it  is  less  than  the  polar  moment  of 
inertia,  often  much  smaller.  Common  values  for  the 
torsional  inertia  have  been  tabulated  in  Reference  5. 

(d)  The  flexural  inertia  is  the  appropriate  moment  of 
inertia  for  the  major  or  minor  axes  depending  upon  the 
orientation  of  the  beam. 


16 


LOCAL  COORDINATE  SYSTEM  FOR  BEAM  ELEMENT 
D.  Elamt  Load  Factors  (4F10.0) 

Nodal  point  loads  (no  aomenta)  due  to  gravity  are  computed.  Thr*® 
lines  need  be  supplied  which  specify  the  fraction  of  these  loads  (in 
each  of  the  three  global  coordinate  directions)  to  be  added  to  each 

eleamnt  load  case. 

f.i—  l:  Multiplier  of  gravity  load  in  the  +X  direction 
Notes  Coluama  Variables  Ss£J2 


1-10 
11  -  20 
21  -  30 
31  -  40 


Eleamnt  load  case  A 
Eleamnt  load  case  B 


Element  load  case  C 
Eleamnt  load  case  D 


Line  2:  As  above  for  gravity  in  the  +Y  direction 
Hu.  3:  As  above  for  gravity  in  the  +Z  direction 
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E.  Fixed-end  Forces  ( 15 , 6F10. 0/ 15 ,6F10. 0) 


Two  lines  are  required  for  each  unique  set  of  fixed-end  forces 
occurring  in  the  analysis.  Distributed  load,  concentrated  loads 
located  between  member  ends,  and  thermal  loads  can  be  specified  using 
the  fixed-end  forces. 

Line  1: 


Notes 


Notes 


Columns 

Variables 

Entry 

1-5 

Fixed-end  force  number 

6  -  15 

Fixed-end  force  in  local  1-direction 
Node  I 

at 

16  -  25 

Fixed-end  force  in  local  2-direction 
Node  I 

at 

26  -  35 

Fixed-end  force  in  local  3-direction 
Node  I 

at 

36  -  45 

Fixed -end  moment  about  local 

1 -direction  at  Node  I 

46  •  55 

Fixed-end  moment  about  local 
2-direction  at  Node  I 

56  -  65 

Fixed-end  moment  about  local 
3-direction  at  Node  I 

Line  2: 

Columns 

Variables 

Entry 

1-5 

Blank 

6-15 

Fixed-end  force  in  local  1-direction 
Node  J 

at 

• 

16  -  25 

Fixed-end  force  in  local  2-direction 
Node  J 

at 

• 

26  -  35 

Fixed-end  force  in  local  3-direction 
Node  J 

at 

36  -  45 

Fixed-end  moswnt  about  local 

1 -direction  at  Node  J 

46  -  55 

Fixed-end  moment  about  local 
2-direction  at  Node  J 

56  -  65 

Fixed -end  moment  about  local 
3-direction  at  Node  J 

18 


(1)  Note  that  values  input  are  literally  fixed-end  values. 

Corrections  due  to  hinges  and  rollers  are  performed  within  the  program. 
Directions  1,  2  and  3  indicate  principal  directions  in  the  local  beam 
coordinates . 

(2)  The  element  loads  for  beams,  usually  referred  to  as  member  loads 
are  input  as  fixed  end  forces  and  moments.  These  are  true  fixed  end 
quantities  without  end  condition  corrections.  Formulas  for  these  fixed  end 
quantities  can  be  found  in  most  structural  engineering  handbooks  and 
textbooks . 

F.  Beam  Element  Data  Lines  (1015,216,18) 

Notes  Columns  Variables  Entry 

1-5  Element  number 


6-10 

Node  number  I 

11  -  15 

Node  number  J 

16  -  20 

Node  number  K  -  see  accompanying 

gure 

21  -  25 

Material  property  number 

26  -  30 

Section  property  number 

31  -  35 

36  -  40 

41  -  45 

46  -  50 

A-i 

j  Fixed-end  force  identification  for 

^  element  load  cases  A,  B,  C,  and  D 

^  respectively 

(1)  51  -  56  End  release  code  at  node  I 

(1)  57  -  62  End  release  code  at  node  J 

(2)  63  -  70  Optional  parameter  k  used  for  automatic 

generation  of  element  data.  This 
option  is  described  below  under  a 
separate  i.^ading.  If  the  option  is  not 
used,  the  field  is  left  blank. 

NOTES: 


(1)  The  end  release  code  at  each  node  is  a  six  digit  number  of  ones 
and/or  zeros.  The  1st,  2nd,  ....  6th  digits  respectively  correspond  to 
the  force  components  Rl,  R2,  R3,  Ml,  M2,  M3,  at  each  node. 


If  any  one  of  the  abcve  element  end  forces  is  known  to  be  zero  (hinge 
or  roller),  the  digit  corresponding  to  that  component  is  a  one. 
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(2)  If  a  series  of  elements  occurs  in  which  each  element  number  NE^ 
is  one  greater  than  the  previous  number 


i.e. , 


NEi  =  NEi_1  +  1 


only  the  element  data  line  for  the  first  element  in  the  series  need  be 
given  as  input,  provided. 

(a)  The  end  nodal  point  numbers  are  NI^  =  +  k 

NJ  =  NJ  +  k 

and  the 


(b)  material  property  number 

(c)  element  property  number 

(d)  fixed-end  force  identification  numbers  for  each  element 
load  case 

(e)  eleme.it  release  code 

(f)  orientation  of  local  2-axis  are  the  same  for  each 
element  in  the  series. 


The  value  of  k,  if  left  blank,  is  taken  to  be  one.  The  element  data 
line  for  the  last  beam  element  must  always  be  given. 

(3)  When  successive  beam  elements  have  the  same  stiffness,  orienta¬ 
tion  and  element  loading,  the  program  automatically  skips  recomputation  of 
the  stiffness.  Use  this  feature  when  numbering  the  beams  to  obtain  maximum 
efficiency. 
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PLANE  STRESS  MEMBRANE  ELEMENTS,  ELEMENT  TYPE  3 

Quadrilateral  (and  triangular)  elements  can  be  used  for  plane  stress 
membrane  elements  of  specified  thickness  which  are  oriented  in  an  arbitrary 
plan.  All  elements  have  temperature-dependent  orthotropic  material  proper¬ 
ties.  In  compatible  displacement  modes  can  be  included  at  the  element 
level  in  order  to  improve  the  bending  properties  of  the  elements. 

A  general  quadrilateral  element  is  shown  below: 


Z 


A  local  element  coordinate  system  is  defined  by  a  u-v  system.  The  v-axis 
coincides  with  the  I-J  side  of  the  element.  The  u  axis  is  normal  to  the 
v-axis  and  is  in  the  plane  defined  by  nodal  points  I,  J  and  L.  Node  K  must 


be  in  the  same  plane  if  the  element 
correct.  The  following  sequence  of 
of  TYPE  3  elements. 

A.  Control  Line  (615) 

Notes  Columns  Variables 

1-5 

6-10 

11  -  15 

16  -  20 

(5)  30 


stiffness  calculations  are  to  be 
lines  define  the  input  data  for  a  set 


Entry 

The  number  3 

Total  number  of  plane  stress  elements 

Number  of  material  property  lines 

Maximum  number  of  temperature  points 
for  any  one  material;  see  Section  B 
below. 

Nonzero  numerical  value  will  suppress 
the  introduction  of  incompatible 
displacement  modes. 
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Material  Property  Information:  Orthotropic,  temperature -dependent 
material  properties  are  possible.  For  each  different  material,  the 
following  group  of  lines  must  be  supplied. 


B.  Material  Property  Identification  Line  (2I5,3F10.0) 


Notes 

Columns 

Variables 

Entry 

1-5 

Material  identification  number 

6-10 

Number  of  different  temperatures  for 
which  properties  are  given.  If  this 
field  is  left  blank,  the  number  is 
taken  as  one 

11  -  20 

Veight  density  of  material  (used  to 
calculate  gravity  loads) 

21  -  30 

Mass  density  (used  to  calculate  mass 
matrix) 

u 

1 

* 

o 

Angle  $  in  degrees,  measured 

counter-clockwise  from  the  v-axis  to 
the  n-axis 


K 


The  n-s  axes  are  the  principal  axes  for  the  orthotropic  material. 
Veight  and  mass  densities  need  be  listed  only  if  gravity  and  inertia  loads 
are  to  be  considered. 
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C.  TWo  Temperature  Data  Lines 
Line  1:  (8F10.0) 


Notes 

Columns 

Variables 

Entry 

1-10 

Temperature 

11  -  20 

Modulus  of  Elasticity  - 

21  -  30 

Modulus  of  Elasticity  -  Es 

31  -  40 

Modulus  of  Elasticity  E^ 

41  -  50 

Strain  Ratio  -  v 

ns 

51  -  60 

Strain  Ratio  -  v  . 

nt 

61  -  70 

Strain  Ratio  -  v 

St 

71  -  80 

Shear  Modulus  *  6 

ns 

Line  2: 

(3F10.0) 

Notes 

Columns 

Variables 

Entry 

1-10 

Coefficient  of  thermal  expansion  -  «n 

11  -  20 

Coefficient  of  thermal  expansion  -  a 

8 

21  -  30 

Coefficient  of  thermal  expansion  - 

All  material  constants  must  always  be  specified.  For  plane  stress, 
the  program  modifies  the  constitutive  relations  to  satisfy  the  condition 
that  the  normal  stress  equals  zero. 

D.  Element  Load  Factors  (5F10.0) 

Four  lines  are  used  to  define  the  element  load  cases  A,  B,  C  and  D  as 
fraction  of  the  basic  thermal,  pressure  and  acceleration  loads. 


First  line,  load  case  A:  Second  line,  load  case  B,  etc. 

Notes  Columns  Variables  Entry 

1  -  10  Fraction  of  thermal  load 

1  -  20  Fraction  of  pressure  load 

21  -  30  Fraction  of  gravity  in  X-direction 

31  -  40  Fraction  of  gravity  in  Y-direction 

Fraction  of  gravity  in  Z-direction 


41  -  50 


E.  Eli 


it  Data  Lines  (615, 2F10. 0,215, F10.0) 


One  line  per  element  must  be  supplied  (or  generated)  with  the 
following  information: 


Notes  Columns  Variables  Entry 


1  -  5 


Element  number 


6-10 
11  -  15 
16  -  20 
(4)  21  -  25 

26  -  30 
(3)  31  -  40 

41  -  50 
(2)  51  -  55 

(1)  56  -  60 

61  -  70 


Node  I 
Node  J 
Node  K 

Node  L  (Node  L  must  equal  Node  K  for 
triangular  element) 

Material  identification  number 

Reference  temperature  for  zero  stresses 
within  element 

Normal  pressure  on  I-J  side  of  element 
Stress  evaluation  option  Mn" 

Element  data  generator  nk'' 

Element  thickness 


NOTES: 

(1)  Element  Data  Generation  -  Element  lines  must  be  in  element  number 
sequence.  If  lines  are  omitted,  data  for  the  omitted  elements  will  be 
generated.  The  nodal  numbers  will  be  generated  with  respect  to  the  first 
line  in  the  series  as  follows: 


I 

n 


1 


+  k 


J  « 

n 


n-1 


+  k 


K 


n 


n-1 


+  k 


All  other  element  information  will  be  set  equal  to  the  information  on 
the  last  line  read.  The  data  generation  parameter  k  is  specified  on  that 
line. 
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(2)  Stress  Print  Option  -  The  following  description  of  the  stress 
print  option  applies  to  both  element  types  3  and  4.  The  value  of  the 
stress  print  option  "n"  can  be  given  as  1,  0,  8,  16,  or  20.  Refer  to 
element  type  4  for  more  discussion  regarding  the  option. 

(3)  Thermal  Data  -  Nodal  temperatures  as  specified  nodal  point  data 
lines  are  used  by  element  type  3  and  4  in  the  following  two  ways  (see 
element  type  4  for  more  information): 

(a)  Temperature-dependent  material  properties  are  approximated 
by  interpolating  (or  extrapolating)  the  input  material 
properties  at  the  temperature  T  corresponding  to  the  origin 
of  the  local  s-t  coordinate  system.  The  material  properties 
throughout  the  element  are  assumed  constant  corresponding  to 
this  temperature. 

(b)  For  computation  of  nodal  loads  due  to  thermal  strains  in  the 
element  a  bilinear  interpolation  expansion  for  the 
temperature  change  AT  (s,t)  is  used. 

(4)  Use  of  triangles  -  In  general,  the  elements  are  most  effective 
when  they  are  rectangular,  i.e.,  the  elements  are  not  distorted.  There¬ 
fore,  regular  and  rectangular  element  mesh  layouts  should  be  used  as  much 
as  possible.  In  particular,  the  triangle  used  is  the  constant  strain 
triangle;  and  it  should  be  avoided,  since  its  accuracy  is  not  satisfactory. 

(5)  Use  of  Incompatible  Modes  -  Incompatible  displacement  modes  have 
been  found  to  be  effective  only  when  used  in  rectangular  elements.  They 
should  always  be  employed  with  care.  Since  incompatible  modes  are  used  for 
all  elements  of  a  group,  it  is  recommended  to  use  separate  element  groups 
for  elements  with  incompatible  modes  and  elements  without  incompatible 
modes,  respectively.  (See  note  (2)  for  the  master  control  line). 


25 


TWO-DIMENSIONAL  FINITE  ELEMENTS,  ELEMENT  TYPE  4 


Quadrilateral  (and  triangular)  elements  can  be  used  as: 

_  Axisymmetric  solid  elements  symmetrical  about  the  Z-axis. 

The  radial  direction  is  specified  as  the  Y-axis.  Care 
must  be  exercised  in  combining  this  element  with  other 
types  of  elements 

_  Plane  strain  elements  of  unit  thickness  in  the  Y-Z  plane. 

_  Plane  stress  elements  of  specified  thickness  in  the  Y-Z 
plane. 

All  elements  have  temperature -dependent  orthotropic  material  proper¬ 
ties.  Incompatible  displacement  modes  can  be  included  at  the  element  level 
in  order  to  improve  the  bending  properties  of  the  element. 

A  general  quadrilateral  element  is  shown  below: 


A.  Control  Line  (615) 

Notea  Columns 

Variables 

Entry 

1  -  5 

The  number  4 

6-10 

Total  number  of  elements 

11  -  15 

Number  of  different  materials 

16  -  20 

Maximum  number  of  temperature  lines  for 
any  one  material  -  see  Section  C  below. 
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Notes  Columns  Variables  Entry 

0  for  axisymmetric  analysis 
25  ■  1  for  plane  strain  analysis 

2  for  plane  stress  analysis 

(5)  30  Nonzero  numerical  value  will  suppress 

the  introduction  of  incompatible  dis¬ 
placement  modes.  Incompatible  modes 
cannot  be  used  for  triangular  elements 
and  are  automatically  suppressed. 


Material  Property  Information:  Orthotropic,  temperature -dependent 
material  properties  are  possible.  For  each  different  material  the 
following  group  of  lines  must  be  supplied. 

B.  Material  Property  Identification  Line  (2I5,3F10.0) 

Notes  Columns  Variables  Entry 

I- 5  Material  identification  number 

6-10  Number  of  different  temperature  for 

which  properties  are  given.  If  this 
field  is  left  blank,  the  number  is 
taken  as  one. 

II- 20  Weight  density  of  material  (used  to 

calculate  gravity  loads) 

21  -  30  Mass  density  (used  to  calculate  mass 

matrix) 

31  -  40  Angle  B  in  degrees,  measured  counter¬ 

clockwise  from  the  v-axls  to  the  n-axis. 


PRINCIPAL  MATERIAL  AXES 


The  n-s  axes  are  the  principal  axes  for  the  orthotropic  material. 
Weight  density  is  needed  only  if  gravity  and  inertia  loads  are  to  be 
considered. 

C.  TVo  Temperature  Data  Lines: 


Line  1: 

(8F10.0) 

Notes 

Columns 

Variables 

Entry 

1  -  10 

Temperature 

11  -  20 

Modulus  of  elasticity  -  En 

21  -  30 

Mclrlus  of  elasticity  -  E 

s 

31  -  40 

Modulus  of  elasticity  -  Et 

41  -  50 

Strain  ratio  -  v_ 

ns 

51  -  60 

Strain  ratio  -  v 

61  -  70 

Strain  ratio  - 

71  -  80 

Shear  modulus  -  G _ 

ns 

Line  2: 

(3F10.0) 

Notes 

Columns 

Variables 

Entry 

1-10 

Coefficient  of  thermal  expansion  -  a 

n 

11  -  20 

Coefficient  of  thermal  expansion  - 

21  -  30 

Coefficient  of  thermal  expansion  -  «t 

All 

material 

constants  must 

always  be  specified.  For  plane  stress. 

the  program  modifies  the  constitutive  relations  to  satisfy  the  condition 
that  the  normal  stress  o  equals  zero. 


D.  Element  Load  Factors  (5F10.0) 

Four  lines  are  used  to  define  the  element  load  cases  A,  B,  C,  and  D 
as  fraction  of  the  basic  thermal,  pressure  and  acceleration  loads. 
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First  line,  load  case  A;  Second  line,  load  case  B;  etc. 
Notes  Columns  Variables  Entry 


1-10 

Fraction  of  thermal  load 

11  -  20 

Fraction  of  pressure  load 

21  -  30 

Fraction  of  gravity  in  X-direction 

31  -  40 

Fraction  of  gravity  in  Y-direction 

41  -  50 

Fraction  of  gravity  in  Z-direction 

E. 

Element  Data  Lines  (6I5,2F10.0,2I5,F10.0)  1 

One 

line  per 

element  must  be  supplied  (or  generated)  with  the  | 

following  information:  1 

Notes 

Columns 

Variables  Entry 

1-5 

Element  number 

6-10 

Node  I 

11  -  15 

Node  J 

16  -  20 

Node  K 

(4) 

21  -  25 

Node  L  (  Node  L  must  equal  Node  K  for 
triangular  elements) 

26  -  30 

Material  identification  number 

(3) 

31  -  40 

Reference  temperature  for  zero  stresses 
within  element 

41  -  50 

Normal  pressure  on  I-J  side  of  element 

(2) 

51  -  55 

Stress  evaluation  option  'V 

(1) 

56  -  60 

Element  data  generator  "k" 

NOTES: 

61  -  70 

Element  thickness  (For  plane  strain  set 
equal  to  1.0  by  program) 

(1)  Element  Data  Generation  -  Element  lines  must  be  in  element  number 
sequence.  If  lines  are  omitted  the  omitted  element  data  will  be  generated. 
Then  nodal  numbers  will  be  generated  with  respect  to  the  first  line  in  the 
series  as  follows: 
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I 

n 


1 


+  k 


J 

n 


*  J 


n- 1 


+  k 


Kn  =  K  ,  +  k 
n-1 


L 

n 


u  , 

n-1 


+  k 


All  other  element  information  will  be  set  equal  to  the  information  on  the 
last  line  read.  The  data  generation  parameter  k  is  given  on  that  line. 


(2)  Stress  Print  Option  -  The  following  description  of  the  stress 
print  option  applies  to  both  element  types  3  and  4.  The  value  of  the 
stress  print  option  "n"  can  be  given  as  1,  0,  8,  16,  or  20. 


0  *  origin  of  natural  s-t  coordinates  (Fig.  5-2).  Points  1,  2,  3,  and  4 
are  midpoints  of  sides.  The  points  at  which  stresses  are  output  depend  on 
the  value  of  n  as  described  in  the  following  table. 


m 

Stresses  output  at 

i 

None 

0 

0 

D 

0,  1 

16 

■SSSM 

20 

1 

B 

30 


The  stresses  at  0  are  printed  in  a  local  y-z  coordinate  system. 

For  element  type  3,  side  I-J  defines  the  local  y-z  axes  in  the  plane  of 
the  element.  For  element  type  4  the  local  y-z  axes  are  parallel  to  the 
global  Y-Z  axes. 


Y-Z 
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For  both  element  types  3  and  4  the  stresses  at  each  edge  midpoint  are 
output  in  a  rectangular  n-p  coordinate  system  defined  by  the  outward  normal 
to  the  edge  (n  axis)  and  the  edge  (p  axis).  The  positive  p  axis  for  points 
1,  2,  3  and  4  is  from  L  to  I,  J  to  K,  I  to  J  and  K  to  L  respectively 
(positive  direction  is  counterclockwise  about  element). 


The  stresses  for  an  element  are  output  under  the  following  headings: 
Sll,  S22,  S12,  S33,  S-MAX,  S-MIN,  ANGLE.  The  normal  stresses:  Sll  and  S22 
and  the  shear  stress  S12  are  as  described  above.  S-MAX  and  S-MIN  are  the 
principle  stresses  in  the  plane  of  .he  element  and  S33  is  the  third  princi¬ 
pal  stress  acting  on  the  plane  of  the  element.  ANGLE  is  the  angle  in 
degrees  from  (1)  the  local  y  axis  at  point  0,  or  (2)  the  n  axis  at  the 
midpoints,  to  the  axis  of  the  algebraically  largest  principal  stress. 

For  triangular  elements  the  stress  print  option  is  as  described  above 
except  that  n=20  is. not  valid.  If  n=20  is  input,  n  will  be  set  to  16  by 
the  program. 

(3)  Thermal  Data  -  Nodal  temperatures  as  specified  nodal  point  data 
lines  are  used  by  element  type  3  and  4  in  the  following  two  ways: 

(a)  Temperature -dependent  material  properties  are 

approximated  by  interpolating  (or  extrapolating)  the 
input  material  properties  at  the  temperature  T 
corresponding  to  the  origin  of  the  local  s-t  coordinate 
system  (see  Figure  below  for  description  of  local  element 
coordinates).  The  material  properties  throughout  the 
element  are  assumed  constant  corresponding  to  this 
temperature. 


(b)  For  computation  of  nodal  loads  due  to  thermal  strains  in 
the  element  a  bilinear  interpolation  expansion  for  the 
temperature  change  AT  (s,t)  is  used. 

4 

AT  (s.t)  =  I  h  (s,t)  T  -  T 
i=l  1  r 

where  T.  are  the  nodal  temperatures  specified  on  the  joint  data  lines, 
Ty  is  tne  reference  stress  free  temperature  and  h.  (s,t)  are  the 
interpolation  functions. 


33 


(4)  Use  of  Triangles  -  In  general,  the  elements  are  most  effective 
when  they  are  rectangular,  i.e.  the  elements  are  not  distorted.  Therefore, 
regular,  and  rectangular  element  mesh  layouts  should  be  used  much  as 
possible.  In  particular,  the  triangle  used  is  the  constant  strain 
triangle;  and  it  should  be  avoided,  since  its  accuracy  is  not  satisfactory. 

(5)  Use  of  Incompatible  Modes  -  Incompatible  displacement  modes  have 
been  found  to  be  effective  only  when  used  in  rectangular  elements.  They 
should  always  be  employed  with  care.  Since  incompatible  modes  are  used  for 
all  elements  of  a  group  it  is  recommended  to  use  separate  element  groups 
for  elements  with  incompatible  modes  and  elements  without  incompatible 
modes,  respectively.  (See  note  (2)  for  the  master  control  line). 


34 


THREE-DIMENSIONAL  SOLID  ELEMENTS,  ELEMENT  TYPE  5 


General  three-dimensional,  eight-node,  isoparametric  elements  with 
three  translational  degrees  of  freedom  per  node  are  identified  by  the 
number  5.  Isotropic  material  properties  are  assumed.  The  element  load 
cases  (A,  B,  C  and  D)  are  defined  as  a  combination  of  surface  pressure, 
hydrostatic  loads,  inertia  loads  in  three  directions  and  thermal  loads. 
The  six  components  of  stress  and  three  principal  stresses  are  computed  at 
the  center  of  each  element.  Also,  surface  stresses  are  evaluated.  Nine 
incompatible  displacement  modes  are  assumed  in  the  formation  of  element 
stiffness  matrices.  For  8-node  elements  without  incompatible  modes  use 
element  type  8. 


A.  Control  Line  (415) 


Notes 


B. 


Columns 

Variables 

Entry 

1  -  5 

The  number  5 

6-10 

Number  of  8-node  solid  elements 

11  -  15 

Number  of  different  materials 

16  -  20 

Number  of  element  distributed  load  sets 

Material  Property  Lines  (I5,4F10.0)  One  line  for  each  different 
material 


Notes 


Columns 

Variables 

Entry 

1-5 

Material  identification  number 

6  -  15 

Modulus  electability  (only  elastic, 
isotropic  materials  are  considered) 

16  -  25 

Poisson's  ratio 

26  -  35 

Weight  density  of  material  (for 
calculation  of  gravity  loads  of  mass 
matrix) 

36  -  45 

Coefficient  of  thermal  expansion 

C.  Distributed  Surface  Loads  (2I5,2F10.2,I5)  One  line  is  required 
for  each  unique  set  of  uniformly  distributed  surface  loads  and 
for  each  reference  fluid  levex  for  hydrostatically  varying 
pressure  loads.  See  notes  (4)  and  (5)  for  sign  convention. 
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Notes 

Columns 

1-5 

Variables 

Entrv 

Load  set  identification  number 

(5) 

6-10 

LT  (load  type) 

LT  =  1  if  this  line  specifies  a 
uniformly  distributed  load. 

LT  =  2  if  this  line  specifies  a 

hydrostatically  varying  pressure. 

(5)  11-20  P 

If  LT  =  1,  P  is  the  magnitude  of  the 
uniformly  distributed  load 
If  LT  =  2,  P  is  the  weight  density  of 
the  fluid  causing  the 
hydrostatic  pressure 

(5)  21-30  Y 

If  LT  =  1,  leave  blank 

If  LT  =  2,  Y  is  the  global  Y  coordinate 
of  the  surface  of  fluid  causing 
hydrostatic  pressure  loading 

(4)  31  -  35  Element  face  number  on  which  surface 

(5)  load  acts.  Face  numbers  are  from  1  to 

6  as  described  in  note  (5)  for 
uniformly  distributed  loads  and  can  be 
only  faces  2,  4  or  6  for  hydrostati¬ 
cally  varying  pressures. 

D.  Acceleration  Dae  to  Gravity  (F10.2) 


Notes 

Columns 

Variables 

Entrv 

1-10 

Acceleration  due  to  gravity  (for 
calculation  of  mass  matrix) 

E.  Element  Load  Case  Multipliers  (5  lines  of  4F10.2) 

Multipliers  on  the  element  load  cases  are  scaling  factors  in 
order  to  provide  flexibility  in  modifying  applied  loads. 


Line  1: 

Notes 

Columns 

Variables 

Entrv 

(7) 

1-10 

PA-i 

11  -  20 

PB  -  Pressure  load 

21  -  30 

PC  multipliers 

31  -  40 

PIV 
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Notes 

(6) 


Notes 


Notes 


Notes 


PA  is  a  factor  used  to  scale  the  complete  set  of  distributed 
surface  loads.  This  scale  set  of  loads  is  assigned  to  element 
load  case  A.  Note  that  zero  is  a  valid  multiplier.  PC,  PC  and 
PD  are  similar  to  PA  except  that  scaled  loads  are  assigned  to 
element  load  cases  B,  C  and  D  respectively.  For  the  majority  of 
applications  these  factors  should  be  1.0 

Line  2: 


Columns 

Variables 

Entry 

1-10 

TA-i 

11  -  20 

TB 

Thermal  load 

21  -  30 

TC 

multipliers 

31  -  40 

TD-' 

TA  is  a  factor  used  to  scale  the  complete  set  of  thermal  loads. 
The  scaled  set  of  loads  are  then  assigned  to  element  load  case  A. 
TB,  TC  and  TD  are  similar  and  refer  to  element  load  cases  B,  C 
and  D  respectively. 


Line  3: 

Columns 

Variables 

Entry 

1  -  10 

GXA-i 

11  -  20 

GXB 

- 

Gravity  load 

21  -  30 

GXC 

multipliers  for  +  X 

31  -  40 

gxdJ 

global  direction 

Line  4: 

Columns 

Variables 

Entry 

1-10 

GYA-i 

11  -  20 

GYB 

- 

Gravity  load 

21  -  30 

GYC 

Multipliers  for  +  Y 

31  -  40 

gydJ 

global  direction 

Line  5: 

Columns 

Variables 

Entry 

1-10 

GZAi 

11  -  20 

GZB 

■ 

Gravity  load 

21  -  30 

GZC 

multipliers  for  +  Z 

31  -  40 

gzdj 

global  direction 
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Gravity  loads  are  computed  from  the  weight  density  of  the  material  and 
from  the  geometry  of  the  element.  GXA  is  a  multiplier  which  reflects  the 
location  of  the  gravity  axis  and  any  load  factors  used.  The  program  com¬ 
putes  the  weight  of  the  element,  multiplies  it  by  GXA  and  assigns  the 
resulting  loads  to  the  +  X  direction  of  element  load  case  A.  Consequently 
GXA  is  the  product  of  the  component  of  gravity  along  the  +  X  global  axis 
(from  -  1.0  to  1.0)  and  any  desired  load  factor.  GXB,  GXC  and  GXD  are 
similar  to  GXA  and  refer  to  element  load  cases  B,  C  and  0  respectively. 

GYA  and  GZA  refer  to  the  global  Y  and  Z  directions  respectively. 

F.  Element  Data  Lines  (1215,412, 211, F10. 2) 

Notes  Columns  Variables  Entry 


1  -  5  Element  number 

6  -  10-j  rl 

11  -  15  Global  node  point  2 

16  -  20  numbers  corresponding  3 

(3)  21  -  25  -  to  element  nodes  ■  4 

26-30  5 

31-25  6 

36-40  7 

41  -  45-*  L8 

(2)  46  -  50  Integration  Order 

51-55  Material  Number 


(1)  56  -  60  Generation  Parameter  (INC) 

61-62  LSA  -i  LSA  is  the  distributed 

63  -  64  LSB  -  surface  load  set  identi- 

65  -  66  LSC  fication  number  of  the 

67  -  68  LSD  -*  distributed  load  acting 

on  this  element  to  be 
assigned  to  element  load 
case  A.  LSB,  LSC  and  LSD 
refer  to  element  load  cases 
B,  C,  and  D,  respectively 

(8)  69  -  70  Face  number  for  stress  output 

71  -  80  Stress-free  element  temperature 

NOTES: 

(1)  Element  Generation 

(a)  Element  lines  must  be  in  ascending  order 

(b)  If  a  series  of  element  lines  are  omitted,  generation  is 
possible  as  follows: 


*  Nodal  point  numbers  are  generated  by  adding  INC  to  those 
of  the  preceding  element.  (If  omitted,  INC  is  set 
equal  to  1. ) 

*  Same  material  properties  are  used  as  for  the  preceding 
element. 

*  Same  temperature  is  used  for  succeeding  elements. 
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•  If  on  first  line  for  the  series  the  integration  order  is: 

>  0  Same  value  is  used  for  succeeding  elements. 

=  0  A  new  element  stiffness  is  not  formed. 

Element  stiffness  is  assumed  to  be  identical 
to  that  of  the  preceding  element. 

<  0  Absolute  value  is  used  for  the  first  element 

of  the  series,  and  the  same  element  stiffness 
is  used  for  succeeding  element. 

*  If  on  first  line  for  the  series,  the  distributed 
load  number  (for  any  load  case)  is: 

>  0  Same  load  is  applied  to  succeeding  elements. 

<  0  The  load  case  is  applied  to  this  element  but 

not  to  succeeding  elements  in  the  series. 

(c)  Element  line  for  the  last  element  must  be  supplied. 

(2)  Integration  Order 

Computation  time  (for  element  stiffness)  increases  with  the  third 
power  of  the  integration  order.  Therefore,  the  smallest  satisfactory  order 
should  be  used.  This  is  found  to  be: 

2  for  rectangular  element 

3  for  skewed  element 

4  may  be  used  if  element  is  extremely  distorted  in  shape, 
but  not  recommended. 

Hesh  should  be  selected  to  give  "rectangular"  elements  as  far  as 
possible. 

(3)  Element  coordinate  System 

Local  element  coordinate  systems  is  a  natural  system  for  this  element 
in  which  the  element  maps  onto  a  cube.  Local  element  numbering  is  shown  in 
the  diagram  below: 
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(4)  Identification  of  Eleaent  Faces 


Element  faces  are  numbered  as  follows: 


Face  1 
2 

3 

4 

5 

6 
0 


corresponds  to 
corresponds  to 
corresponds  to 
corresponds  to 
corresponds  to 
corresponds  to 
corresponds  to 


+  a  direction  - 

-  a  direction 
+  b  direction 

-  b  direction 
+  c  direction 

-  c  direction  ■ 
the  center  of  the 


Faces  1,3,5  are 
positive  faces 

Faces  2,4,6  are 
negative  faces 

eleaent 


(5)  Distributed  Surface  Loads 


TVo  types  of  surface  loadings  nay  be  specified;  load  type  1  (LT  =  1), 
uniformly  distributed  surface  load  and  load  type  2  (LT  *  2),  hydrostatically 
varying  surface  pressure  (but  not  surface  tension).  Both  loading  types  are 
for  loads  nomal  to  the  surface  and  do  not  include  surface  shears.  Surface 
loadings  that  do  not  fall  into  these  categories  aust  be  input  as  nodal 
loads  on  the  concentrated  load  data  lines  (see  Section  V). 

(1)  LT  *  1:  A  positive  surface  load  acts  in  the  direction 
of  the  outward  nomal  of  a  positive  eleaent  face  and 
along  the  inward  nomal  of  a  negative  eleaent  face  as 
shown  in  the  following  diagram. 


If  the  uniformly  distributed  surface  loading  P  is  input  as  a  positive 
quantity  then  it  describes  pressure  loading  on  faces  2,  4  or  6  and  tensile 
loading  on  faces  1,  3  or  5.  If  P  is  input  as  a  negative  quantity  then  it 
describes  tensile  loading  on  faces  2,  4  or  6  and  pressure  on  faces  1,  3  or  5. 
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(2)  LT  =  2:  A  hydrostatically  varying  surface  pressure 
on  element  faces  2,  4  or  6  can  be  specified  by  a 
reference  fluid  surface  and  a  fluid  weight  density  Y  as 
input.  Only  one  hydrostatic  surface  pressure  line  need 
be  input  in  order  to  specify  a  hydrostatic  loading  on 
the  complete  structure.  The  consistent  nodal  loads  are 
calculated  by  the  program  as  follows.  At  each  numerical 
integration  point  "i"  on  an  element  surface  the  pressure 
is  calculated  from 


Pi  =  Y  (Yi 


Yref> 


Where  is  the  global  Y  coordinate  of  the  point  in  question  and 
Y  ^  specifies  the  fluid  surface  assuming  gravity  acts  along  the  -Y  axis 


If  P^  >  0,  corresponding  to  surface  tension,  the  contribution  is 
ignored,  tf  an  element  face  is  such  that  Y^,  >  Yf  _  for  all  i  (16 
integration  points  are  used  by  program)  then  no  nodal  loads  will  be  applied 
to  the  element.  If  some  P^  >  0  and  some  P,  <  0  for  a  particular  face,  then 
approximate  nodal  loads  are  obtained  for  tne  partially  loaded  surface. 

(6)  Thermal  Loads 

Thermal  loads  are  computed  assuming  a  constant  temperature  increase  AT 
throughout  the  element. 


T 

avg 


=  the  average  of  the  8  nodal  point  temperatures 
specified  on  nodal  point  data  lines 


T  =  stress  free  element  temperature  specified  on  the 
element  line 
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(7)  Element  Load  Cases 

Element  load  case  A  consists  of  all  the  contributions  from  distributed 
loadings,  thermal  loadings  and  gravity  loading  for  all  the  elements  taken 
collectively. 

Load  case  A  s  Z  (PA  x  pressure  loading 
+  TA  x  thermal  loading 
+  GXA  x  gravity  X  loading 
+  GYA  x  gravity  Y  loading 
+  GZA  x  gravity  Z  loading) 

Element  load  case  A  for  the  set  of  three  dimensional  solid  elements  is 
added  to  element  load  case  A  for  the  other  element  types  in  the  analysis. 

The  treatment  of  element  load  cases  B,  C  and  D  is  analogous  to  that  of  ele¬ 
ment  load  case  A.  The  loading  cases  for  the  structure  are  obtained  by  add* 
ing  linear  combinations  of  element  load  cases  A,  B,  C  and  D  to  the  nodal 
loads  specified  on  the  joint  data  lines. 

(8)  Output  of  Element  Stresses 

(a)  At  the  centroid  of  the  element,  stresses  are  referred  to  the 
global  axes.  Three  principal  stresses  are  also  presented. 

(b)  At  the  center  of  an  element  face,  stresses  are  referred  to  a 
set  of  local  axes  (x,y,z).  These  local  axes  are  individually 
defined  for  each  face  as  follows:  Let  nodal  points  I, 

J,  K  and  L  be  the  four  corners  of  the  eleaient  face.  Then: 

x  is  specified  by  LI  -  JK,  where  LI  and  JK  are 
midpoints  of  sides  L-I  and  J-K 

z  is  normal  to  x  and  to  the  line  Joining  midpoints  IJ 
and  KL 

y  is  normal  to  x  and  z,  to  coaplete  the  right-handed 
system 
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The  corresponding  nodal  points  1,  J,  K  and  L  in  each  face  are  given  in 
the  table. 


FACE 

NODAL  POINTS 

I 

J 

K 

L 

1 

1 

2 

6 

5 

2 

4 

3 

7 

8 

3 

3 

7 

6 

2 

4 

4 

8 

5 

1 

5 

8 

5 

6 

7 

6 

4 

1 

2 

3 

IWo  surface  principal  stresses  and  the  angle  between  the  algebraically 
largest  principal  stress  and  the  local  x  axis  are  printed  with  the  output. 
It  is  optional  to  choose  one  or  two  locations  of  an  element  where  stresses 
are  to  be  computed.  In  the  output,  "face  zero”  designates  the  centroid  of 
the  element. 
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PLATE  AND  SHELL  QUADRILATERAL  EIEMENT,  ELEMENT  TYPE  6 

Quadrilateral  or  triangular  shaped  thin  plate/shell  elements  are 
identified  as  eleaent  type  "6".  For  the  quadrilateral,  four  nodes  I,  J,  X, 
and  L  way  be  located  arbitrarily  in  the  global  (X,Y,Z)  system.  A  local 
eleaient  coordinate  system  (x,y,z)  is  defined  as  follows: 


LOCAL  COORDINATE  SYSTEM  FOR  THE  3-DIMENSIONAL  THIN  PLATE/SHELL  ELEMENT 

1.  The  x-axis  la  specified  by  the  line  LI-JK,  where  LI  and  JK  are  the 

midpoints  of  sides  L-I  and  J-K,  respectively. 

2.  The  z-axis  is  normal  to  the  x-direction  and  to  the  line  joining 
midpoints  IJ  and  KL. 

3.  The  y-axis  is  normal  to  both  the  z  and  x-directions  to  complete 
the  remainder  of  the  right-handed  system. 

This  local  system  is  used  to  express  all  physical  and  kinematic  shell 
properties  (stresses,  strains,  material,  law,  etc.),  except  that  gravity 
forces  act  in  the  global  (X,Y,Z)  coordinate  system. 

For  the  analyses  of  smooth  shells,  rotational  constraints  normal  to 
the  surface  may  be  imposed  by  the  addition  of  boundary  elements  at  the 
nodes  (i.e.,  element  type  "7"). 

Plate/shell  elements  are  described  by  the  following  sequence  of  data 
lines. 
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A. 

Notes 


B. 


Notes 


Notes 


Control  Line  (315) 


Columns 

Variables 

Entrv 

1-5 

The  number  6 

6-10 

Number  of  shell  elements 

11  -  15 

Number  of  different  materials 

Material  Property  Infon 

lation 

Anisotropic  material  properties  are  permitted.  For  each 
different  material,  two  lines  must  be  supplied. 

Line  1: 

(I10,20X,4F10. 

°) 

Columns 

Variables 

Entrv 

1-10 

Material  identification  number 

31  -  40 

Mass  density 

41  -  50 

Thermal  expansion  coefficient  e^ 

51  -  60 

Thermal  evened™  coefficient  „y 

61  -  70 

Thermal  expansion  coefficient 

Line  2: 

(6F10.0) 

Columns 

Variables 

Entrv 

1-10 

11  -  20 

21  -  30 

31  -  40 

41  -  50 

51  -  60 

Elasticity  eleamnt  C 

Elasticity  element  C** 

Elasticity  element  C*^ 

Elasticity  element  C*8 

Elasticity  elearant  C^ 

Elasticity  element  Gys 

*y 

Elements  in  plane  stress 
material  matrix  [C] 

0 

c  c  c 

c 

o“ 

a 

c*x  cxy  cxs 

e** 

Tyy 

c*y  cyy  Gy® 

Yyy 

xsj 

xs  ys  xy 

[  *yJ 
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C.  Element  Load  Multipliers  (5  lines) 


Line  1: 

(4F10.0) 

Notes 

Columns 

Variables 

Entry 

1-10 

Distributed 
load  case  A 

lateral  load  multiplier  for 

11  -  20 

Distributed 
load  case  B 

lateral  load  multiplier  for 

21  -  30 

Distributed 
load  case  C 

lateral  load  multiplier  for 

31  -  40 

Distributed 
load  case  D 

lateral  load  multiplier  for 

Line  2: 

(4F10.0) 

Notes 

Columns 

Variables 

Entry 

1-10 

Temperature 

multiplier  for  load  case  A 

11  -  20 

Temperature 

multiplier  for  load  case  B 

21  -  30 

Temperature 

multiplier  for  load  case  C 

31  -  40 

Tesiperature 

multiplier  for  load  case  D 

Line  3: 

(4F10.0) 

Notes 

Columns 

Variables 

Entry 

1-10 

X-direction 
case  A 

acceleration  for  load 

11  -  20 

X-direction 
case  B 

acceleration  for  load 

21  -  30 

X-direction 

case  C 

acceleration  for  load 

31  -  40 

X-directlon 
case  D 

acceleration  for  load 

Line  4: 

(4F10.0)  Same 

as  Line  3  for  Y-direction 

Line  5: 

(4F10.0)  Some 

as  Line  3  for  Z-direction 

D.  Kl— ant  Data  Lines  (8I5.F10.0) 
One  line  for  each  elenent 


Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

Element  number 

(2) 

6-10 

Node  I 

11  -  15 

Node  J 

16  -  20 

Node  K 

21  -  25 

Node  L 

(3) 

26  -  30 

Node  0 
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Notes 

Columns  Variables 

Entrv 

31 

-  35 

Material  identification  (if  left  blank, 
taken  as  one) 

Cl) 

36 

-  40 

Element  data  generator  K 

Element  thickness 

41 

-  50 

51 

-  60 

Distributed  lateral  load  (pressure) 

61 

-  70 

Mear  temperature  variation  T  from  the 
reference  level  in  undeformed  position 

NOTES: 

71 

-  80 

Mean  temperature  gradient  9T/3z  across 
the  shell  thickness  (a  positive 
temperature  gradient  produces  a 
negative  curvature) 

CD 

Element  Data  Generation 

Element  lines  must  be  in  element  number  sequence.  If  element  data 
lines  are  omitted,  the  program  automatically  generates  the  omitted 
information  as  follows: 

The  increment  for  element  number  is  one 

i.e. :  NEi+1  =  NE^  +  1 

The  corresponding  increment  for  nodal  number  is  Kn 


ni.a. 

=  NI , 

+ 

K 

i+1 

i 

n 

nj,a. 

=  NJ. 

+ 

K 

i+1 

i 

n 

NK . . . 

=  NK1 

+ 

K 

i+1 

n 

NL, . , 

=  NL, 

+ 

K 

i+1 

i 

n 

Material  identification,  element  thickness,  distributed  lateral  load, 
temperature  and  temperature  gradient  for  generated  elements  are  the  same. 
Always  include  the  complete  last  element  line. 

(2)  Nodal  Points  and  Coordinate  Systems 

The  nodal  point  numbers  I,  J,  K  and  L  are  in  sequence  in  a 
counter-clockwise  direction  around  the  element.  The  local  element 
coordinate  system  (x,  y,  z)  is  defined  as  follows: 

x  Specified  by  LI  -  JK,  where  LI  and  JK  are  midpoints  of 
sides  L-I  and  J-K. 

z  Normal  to  x  and  to  the  line  joining  midpoints  IJ  and  KL. 

y  Normal  to  x  and  z  to  complete  the  right-handed  system. 
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This  system  is  used  to  express  all  physical  and  kinematic  shell 
properties  (stresses,  stains,  material  law,  etc.),  except  that  the 
body  force  density  is  referred  to  the  global  coordinate  system  (X,Y,Z). 


2 


J 


For  analysis  of  flat  plates,  the  stiffness  associated  with  the 
rotation  normal  to  the  plate  surface  is  not  defined,  therefore,  the 
rotation  normal  degree  of  freedom  must  not  be  included  in  the  analysis. 

For  the  analyses  of  shallow  shells,  rotational  constraints  normal  to 
the  surface  may  be  imposed  by  the  addition  of  boundary  elements  at  the 
nodes  (element  type  #7). 

(3)  Node  0 

When  columns  26  -  30  are  left  blank,  mid-node  properties  are  computed 
by  averaging  the  four  nodes.  The  thin  plate  and  shell  element  is  a  quadri¬ 
lateral  of  arbitrary  geometry  formed  from  four  compatible  triangles.  The 
central  node  is  located  at  the  average  of  the  coordinate  of  the  four  corner 
nodes  if  a  node  k  is  not  specified. 

(A)  Element  Stress  Calculations 


Output  are  moments  per  unit  length  and  membrane  stresses. 


MAXIMUM 

STRESS 

ELEMENT 

NUMBER  OF 

COMPONENT 

OUTPUT 

TYPE 

COMPONENTS 

NUMBER 

SYMBOL 

D  E  S  C  R 

I  P  T  I  0 

6 

6 

(1) 

(XX-S/R) 

XX-Stress 

Resultant 

(2) 

(YY-S/R) 

YY-Stress 

Resultant 

(3) 

(XY-S/R) 

XY-Stress 

Resultant 

(A) 

(XX-M/R) 

XX -Moment 

Resultant 

(5) 

(YY-M/R) 

YY-Moment 

Resultant 

(6) 

(XY-M/R) 

XY-Moment 

Resultant 
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BOUNDARY  ELEMENTS,  ELEMENT  TYPE  7 


This  element  is  used  to  constrain  nodal  displacement  to  specified 
values,  to  compute  support  reactions  and  to  provide  linear  elastic  supports 
to  nodes.  If  the  boundary  condition  code  for  a  particular  degree  of  free¬ 
dom  is  specified  as  1  on  the  structure  nodal  point  data  lines,  the  dis¬ 
placement  corresponding  to  that  degree  of  freedom  is  zero  and  no  support 
reactions  are  obtained  with  the  printout.  Alternatively,  a  boundary 
element  can  be  used  to  accomplish  the  same  effect  except  that  support 
reactions  are  obtained  since  they  are  equal  to  member  end  forces  of  the 
boundary  elements  which  are  printed.  In  addition  the  boundary  element  can 
be  used  to  specify  nonzero  nodal  displacements  in  any  direction  which  is 
not  possible  using  the  nodal  point  data  lines. 

The  boundary  element  is  defined  by  a  single  directed  axis  though  a 
specified  nodal  pint,  by  a  linear  extensional  stiffness  along  the  axis  or 
by  a  linear  rotational  stiffness  about  the  axis.  The  boundary  element  is 
essentially  a  spring  which  can  have  axial  displacement  stiffness  and  axial 
rotational  stiffness.  There  is  no  limit  to  the  number  of  boundary  elements 
which  can  be  applied  to  any  joint  to  produce  the  desired  effects.  Boundary 
elements  have  no  effect  on  the  size  of  the  stiffness  matrix. 

A.  Control  Line  (215) 

Notes  Columns  Variables  Entry 

(4)  1-5  The  number  7. 

6-10  Total  number  of  boundary  elements . 

B.  Element  Load  Multipliers  (4F10.0) 

Notes  Columns  Variables  Entry 

1-10  Multiplier  for  load  case  A 

11  -  20  Multiplier  for  load  case  B 

21  -  30  Multiplier  for  load  case  C 

31  -  40  Multiplier  for  load  case  D 

C.  Element  Lines  (8I5,3F10.0) 

One  line  per  element  (in  ascending  nodal  point  order)  except 

where  automatic  element  generation  is  used. 
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Notes  Columns  Variables  Entry 

1-5  Node  N,  at  which  the  element  is  placed 

(1)  6-10 
11  -  15 

16  -  20 

21  -  25 

(2)  26  -  30 

31  -  35 

(3)  36  -  40 

41  -  50 

51  -  60 

61  -  70 


NOTES: 

(1)  Direction  of  boundary  element 

The  direction  of  the  boundary  element  at  node  is  specified  in  one  of 
two  ways. 

(a)  A  second  nodal  point  I  defines  the  direction  of  the 
element  from  node  N  to  node  I. 

(b)  Four  nodal  points  I,  J,  K  and  L  specify  the  direction  of  the 
element  as  the  normal  to  the  plane  defined  by  two 
intersecting  straight  lines  (vectors  a  and  b,  see  Fig. 
below).  The  direction  vector  n  is  given  by  the  vector  cross 
product  n  =  a  x  b  and  defines  the  direction  of  the  boundary 
element. 


a 


K 

n«axb 

rotational  constraint 

IN  THIN  SHELL  ANALYSIS 

The  four  points  I,  J,  K  and  L  need  not  be  unique.  A  useful 
application  for  the  analysis  of  shallow  thin  shells  employs  the  boundary 
element  to  approximate  rotational  constraint  about  the  surface  normal  as 
shown  above. 


Node  In 

Node  J  -  Leave  columns  11-25  blank 
Node  K  if  only  node  I  is  needed. 

Node  IP 

Code  for  displacement 
Code  for  rotation 
Data  generator  K 

Specified  displacement  along  element  axis 
Specified  rotation  about  element  axis 
Spring  stiffness  (set  to  10  if  left 
blank)  for  both  extension  and  rotation. 
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Note  that  node  I  in  case  (a)  and  nodes  I,  J,  K  and  L  in  case  (b)  are  used 
only  to  define  the  direction  of  the  element  and  if  convenient  may  be  any 
nodes  used  to  define  other  elements.  However  'artificial  nodes'  may  be 
created  to  define  directions  of  boundary  elements.  These  'artificial  nodes' 
are  input  on  the  nodal  point  data  lines  with  their  coordinates  and  with  all 
the  boundary  condition  codes  specified  as  fixed. 

It  should  be  noted  that  node  N  is  the  structure  node  to  which  the 
boundary  element  is  attached.  In  case  (a),  a  positive  displacement  moves 
node  N  towards  node. I.  Correspondingly,  a  positive  force  in  the  element 
means  compression  in  the  tiement.  In  case  (b),  a  positive  displacement 
moves  node  N  in  the  positive  direction  of  n  (see  Fig.). 

(2)  Displacement  and  rotation  codes 

Displacement  code  =  1:  When  this  code  is  used,  the  displacement  6, 
specified  in  columns  41-50,  and  the  spring  stiffness  k,  specified  in 
columns  61-70,  are  used  by  the  program  in  the  following  way.  The  load  P, 
evaluated  from  P  =  kfi,  is  applied  to  node  N  in  the  direction  node  N  to  node 
I  in  case  (a)  and  in  the  direction  n  in  case  (b),  if  6  is  positive.  If  k 
is  much  greater  than  the  stiffness  of  the  structure  at  node  N  without  the 
boundary  element,  then  the  net  effect  is  to  produce  a  displacement  very 
nearly  equal  to  6  at  node  N.  If  6  =  0,  then  P  =  0  and  the  stiff  spring 
approximates  a  rigid  support.  Note  that  the  load  P  will  contribute  to  the 
support  reaction  for  nonzero  6.  The  boundary  condition  codes  specified  on 
the  structure  nodal  point  data  lines  must  be  consistent  with  the  fact  that 
a  load  P  is  being  applied  to  node  N  to  effect  the  desired  displacement 
(even  when  this  displacement  is  zero). 

Rotation  code  =  1:  This  case  is  analogous  to  the  situation  described 
above.  A  torque  T,  evaluated  from  T  =  k  8,  is  applied  to  node  N  about  the 
axis  (direction)  of  the  element.  The  rotation  0  is  specified  in  columns 
51-60. 

(3)  Data  generator  Kr 

When  a  series  of  nodes  are  such  that: 

(a)  All  have  identical  boundary  elements  attached 

(b)  All  boundary  elements  have  same  direction 

(c)  All  specified  displacements  and  rotations  are  identical 

(d)  The  nodal  sequence  form  an  arithmetic  sequence,  i.e..,  N, 

N  +  K  ,  N  +  2K  etc. , 
n  n 

then  only  the  first  and  last  node  in  the  sequence  need  be  input.  The 
increment  is  input  in  columns  36-40  of  the  first  line. 
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(4)  Element  load  multipliers 


Each  of  the  four  possible  element  load  cases  A,  B,  C  and  D  associated 
with  the  boundary  elements  consists  of  the  complete  set  of  displacements  as 
specified  on  the  boundary  element  lines  multiplied  by  the  element  load 
multiplier  for  the  corresponding  load  case.  As  an  example,  suppos^that 
displacement  of  node  N  is  specified  as  1.0,  spring  stiffness  as  10  and  no 
other  boundary  element  displacements  are  specified.  Let  case  A  multiplier 
be  0.0  and  case  B  multiplier  be  2.0  For  element  load  case  A  the  specified 
displacement  is  0.0  x  1.0  =  0.0  while  that  for  B  is  2.0  x  1.0  =  2.0. 

Linear  combinations  of  element  load  cases  A,  B,  C  and  D  for  all  types  of 
elements  collectively  for  a  particular  problem  are  specified  on  the  struc¬ 
ture  element  load  multiplier  lines.  As  far  as  the  boundary  element  is 
concerned,  this  device  is  useful  when  a  particular  node  has  a  support 
displacement  in  one  load  case  but  is  fixed  in  others. 


(5)  Recommendations  for  use  of  boundary  elements 


If  a  boundary  element  is  aligned  with  a  global  displacement  direction, 
only  the  corresponding  diagonal  element  in  the  stiffness  matrix  is  modified. 
Therefore,  no  stiffness  matrix  ill-conditioning  results.  However,  when  the 
boundary  element  couples  degrees  of  freedom,  large  off-diagonal  elements 
introduce  ill-conditioning  into  the  stiffness  matrix  which  can  cause  solu¬ 
tion  difficulties. 


In  the  analysis  of  shallow  shells  boundary  elements  with  stiffness  a 
fraction  of  the  element  bending  stiffness  should  be  used  (say  less  than  or 
about  10  percent). 

In  dynamic  analysis  "artificially  stiff"  boundary  elements  should  not 
be  used  (see  note  (8)  in  the  Mode  Shapes  and  Frequencies  section). 
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THICK  SHELL  AND  THREE-DIMENSIONAL  ELEMENTS,  ELEMENT  TYPE  8 


A  minimum  of  8  and  a  maximum  of  21  nodes  are  used  to  describe  a  general 
three  dimensional  isoparametric  element;  the  element  is  used  to  represent 
orthotropic,  elastic  media.  The  element  type  is  identified  by  the  number 
eight  (8).  Three  translational  degrees  of  freedom  are  assigned  to  each 
node,  and  at  least  the  eight  corner  nodes  must  be  input  to  define  a  hexahe¬ 
dron.  Input  of  nodes  9  to  21  is  optional;  the  figures  below  illustrate 
some  of  the  most  commonly  used  node  combinations. 

Element  load  cases  (A,B,C,...)  are  formed  from  combinations  of  applied 
surface  pressure,  hydrostatic  loads,  inertia  loads  in  the  three  directions 
X,Y,Z  and  thermal  loads.  Six  global  stresses  are  output  at  up  to  seven  (7) 
locations  within  the  element;  these  output  locations  are  selected  by  means 
of  appropriate  data  entries. 

Node  temperatures  input  with  the  nodal  point  data  are  used  to  form  an 
average  element  temperature,  which  is  the  basis  of  material  property 
selection  for  the  element.  If  thermal  loads  are  applied,  node  temperatures 
are  used  to  establish  the  temperature  field  within  the  element,  and  the 
temperature  interpolation  functions  are  the  same  as  those  assumed  to 


represent 

element  displacements. 

A. 

Control  Line  (1015) 

Notes 

Columns 

1  -  5 

Variables 

Entry 

Enter  the  number  "8" 

6-10 

NS0L21 

Number  of  solid  elements 

11  -  15 

NUMMAT 

Number  of  different  materials 

(1) 

16  -  20 

MAXTP 

Maximum  number  of  temperature  points 
used  in  the  table  for  any  material; 

EQ.  0  the  default  is  set  to  "l" 

(2) 

21  -  25 

N0RTH0 

Number  of  different  sets  of  martial  axis 
orientation  data; 

EQ.  0;  all  properties  are  defined  in 
the  X,Y,Z,  system 

(3) 

26  -  30 

NDLS 

Number  of  different  distributed  load 
(i.e.,  pressure)  sets 

(4) 

31  -  35 

MAXN0D 

Maximum  number  of  nodes  used  to  describe 
any  one  element; 

8  S  MAXN0D  <  21 

EQ.  0  the  default  is  set  to  "21" 

(5) 

36  -  40 

N0PSET 

Number  of  sets  of  data  requesting  stress 

output  at  various  element  locations; 
EQ.  0;  centroid  output  only 


53 


A.  Control  Line  (1015)  (continued) 

Notes  Columns  Variables  Entry 

(6)  41  -  45  INTRS  Standard  Integration  order  for  the  natural 

(r,s)  directions; 

2  £  INTRS  S  4 

EQ.  0  the  default  is  set  to  "2" 

(6)  46-50  . INTT  Standard  integration  order  for  the  natural 

(t)-  direction; 

2  £  INTT  £  4 

EQ.  0  the  default  is  set  to  "2" 

NOTES: 

(1)  The  variable  MAXTP  limits  the  number  of  temperature  points  that 
can  be  input  for  any  one  of  the  NUMMAT  material  sets;  i.e.,  the  variable 
NTP  on  the  material  property  line  (section  B  below)  cannot  exceed  the  value 
of  MAXTP. 

(2)  N0RTH0  specifies  the  number  of  material  axes  orientation  lines  to 
be  read  and  if  omitted,  all  orthotropic  material  axes  are  assumed  to 
coincide  with  the  global  cartesian  axes  X,Y,Z. 

(3)  NDLS  specifies  the  number  of  line  pairs  of  distributed  surface 
load  data  to  be  read.  NDLS  must  be  a  positive  integer  if  any  pressure 
loads  are  to  be  applied  to  solid  element  faces. 

(4)  MAXNOD  specifies  the  maximum  number  of  nonzero  node  numbers 
assigned  to  any  one  of  the  NS0L21  elements  input  in  Section  F  Element 
Lines.  Locations  of  the  element's  21  possible  nodes  are  shown  in  the 
figure  below  in  which  the  element  is  shown  mapped  into  its  natural  r,s,t 
coordinate  system.  The  eight  corner  nodes  must  be  input  for  every  element, 
and  nodes  9  to  21  are  input  optionally.  If  MAXNOD  is  9  or  greater,  all  21 
node  entries  are  read  for  each  element  (Lines  2  and  3,  Section  7),  but  only 
the  first  MAXNOD  nonzero  entries  encountered  when  reading  in  sequence  from 
1  to  21  will  be  used  for  element  description.  As  an  example,  for  the 
16-17-  and  20-node  elements  MASNOD  has  values  of  16,  17,  20  respectively. 

(5)  As  a  means  of  controlling  the  amount  of  solution  output,  stress 
output  location  sets  are  defined  in  Section  D,  and  the  total  number  of 
these  output  requests  is  specified  by  the  variable  NOPSET.  For  the  case  of 
NOPSET  equals  0,  no  data  is  input  in  Section  D,  and  the  only  stress  output 
produced  by  the  program  is  at  the  element  centroid.  Otherwise,  stress 
output  can  be  requested  at  up  to  seven  (7)  locations  (selected  from  a  table 
of  27  possible  locations)  by  means  of  the  data  entries  given  in  Section  D. 
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(6)  The  entries  INTRS  and  INTT  control  the  number  of  integration 
points  to  be  used  in  numerical  evaluation  of  integrals  over  volumes  in  the 
(r,s)  and  (t)-coordinate  directions,  respectively.  When  solid  elements  are 
used  to  represent  shell  structures,  the  through-the-thickness  integration 
(i.e.,  in  the  natural  t-axis  direction)  can  be  evaluated  less  accurately 
than  those  in-plane  (i.e.,  in  the  r,s  plane).  For  this  case  INTRS  might  be 
3  and  INTT  would  be  chosen  typically  as  2.  The  entries  INTRS  and  INTT  are 
standard  or  reference  values  and  are  used  if  the  integration  order  entries 
on  the  element  lines  (Line  1,  Section  F)  are  omitted.  Nonzero  entries  for 
integration  order(s)  given  on  the  element  lines  over-ride  the  standard 
values  posted  on  this  line. 

B.  Material  Property  Lines 

Orthotropic,  temperature  dependent  material  properties 
are  allowed.  For  each  different  material  that  is  requested 
on  the  Control  Line,  (Section  A)  the  following  data  set  must  be 
supplied  (i.e.,  NUMMAT  sets  the  total  number  of  data  sets): 

Material  Identification  Line  (2I5,2F10.0,6A6) 


Notes 

Columns 

Variables 

Entrv 

(1) 

1  -  5 

M 

Material  identification  number; 
liHi  NUMMAT 

6-10 

NTP 

Number  of  different  temperatures  at 
which  properties  are  given; 

LE.MAXTP 

EQ.  0  the  default  is  set  to  "l" 

CM 

11  -  20 

WIDEN 

Weight  density  of  the  material  used  to 
computed  static  gravity  loads 

21  -  30 

MASSDN 

Mass  density  of  the  material  used  to 
compute  the  mass  matrix  in  a  dynamic 

analysis; 

EQ.  0;  default  set  to  "WIDEN/ 386. 4” 

31  -  66 

Material  description  used  to  label  the 
output . 

NOTES: 

(1) 

Material 

numbers  (M)  must 

be  input  in  ascending  sequence 

beginning  with  1 

and  ending  with 

"NUMMAT";  omissions  or  repetitions  are 

not  permitted. 

(2)  Weight  density  is  used  to  compute  static  node  forces  due  to 
applied  gravity  loads;  mass  density  is  used  to  calculate  element  mass 
matrices  for  use  in  connection  with  a  dynamic  analysis. 

Material  Lines  (7F10.0,6F10.0) 

NTP  pairs  of  lines  are  input  in  order  of  algebraically  increasing 
value  of  temperature. 
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Material  Line  1 


Notes  Columns  Variables 


(1)  1  -  10 
(2)  11  -  20 

21  -  30 
31  -  40 
41  -  50 
51-60 
61  -  70 

Material  Line  2 

Notes  Columns  Variables 


1-10 
11  -  20 
21  -  30 
31  -  40 
41  -  50 
51  -  60 

NOTES: 

(1)  The  12  entries  following  the  temperature  value  T  are  physical 
properties  known  at  T  .  When  two  or  more  temperature  points  describe  a 
material,  interpolation  based  on  average  element  temperature  is  performed 
to  establish  a  property  set  for  the  element.  Hence,  the  range  of  tempera¬ 
ture  points  for  a  material  table  must  span  the  expected  range  of  average 
element  temperatures  for  all  elements  associated  with  the  material. 

(2)  The  12  constants  (E. . ,E._, . . . . ,a.)  are  defined  with  respect  to  a 
set  of  axes  (Xj,X„,X3)  which  are  the  principal  material  directions  for  an 
orthotropic,  elastic  medium.  The  stress-strain  relations  with  respect  to 
the  (X^XgjX^)  system  is  written  as  follows: 


'*11 

1/En 

"V12/E22 

"  V13/E33  ° 

0 

0 

O 

I-* 

*22 

‘v21/E11 

-1/e22 

-  V23/E33  0 

0 

0 

°22 

*33 

-v31/EU 

’V32/E22 

-vl/E33  0 

0 

0 

°33 

Y 

12 

0 

0 

°  1/G12 

0 

0 

*12 

Y23 

0 

0 

0  0 

!/G23 

0 

T23 

Y31 

0 

0 

0  0 

0 

I/G13 

. 

-  {ATOj  AT«2  AT<*3  0  0  0JT 


Entry 


Temperature,  T 


E11 

E22 

In 

vl3 

V23 


at 

at 

at 

at 

at 

at 


T 

T" 

T° 

T" 

T11 

T* 

n 


n 


Entry 


G12 

G13 

G23 

o3 


at 

at 

at 

at 

at 

at 


T 

T° 

T“ 

t* 

T" 

n 


59 


where  e,  and  o,.  are  normal  strains  and  stresses  in  the  X,  directions;  Y 

and  t..  are  shear  strains  and  stresses  on  the  principal  material  planes; 

are  tne  coefficients  of  thermal  expansion,  and  AT  is  the  increase  in  1 

temperature  from  stress  free  (reference  temperature  (T  )  distributed  over 

the  element  volume  specified  on  the  element  data  line  (i.e.,  AT  -  T  -  T  ). 

no 

Material  Axes  Orientation  Sets  (415) 


If  NORTHO  is  zero  on  the  Control  Line,  skin  this  data 
section,  and  all  material  axes  (X^X^tX^)  will  be  assumed  to 
coincide  with  the  global  cartesian  system  (X,Y,Z).  Otherwise, 
NORTHO  lines  must  be  input  as  follows: 


Notes 

Columns 

Variables 

Entry 

(1) 

1  - 

5 

M 

Identification  number; 
GE.  1  and  LE. NORTHO 

(2) 

6  - 

10 

NI 

Node  number  for  point 

11  - 

15 

NJ 

Node  number  for  point 

16  - 

20 

NX 

Node  number  for  point 

NOTES: 

(1)  Identification  number  (M)  must  be  input  in  increasing  sequence 
beginning  with  "l"  and  ending  with  NORTHO”. 


(2)  Orthotropic  material  axes  orientations  are  specified  by  means  of 
the  three  node  numbers  NI,NJ,NK.  For  the  special  case  where  orthotropic 
material  axes  coincide  with  the  global  axes  (X,Y,Z),  it  is  not  necessary  to 
input  data  in  this  section;  see  Section  F,  note  (4).  Let  f j, f_ , be  the 
three  orthogonal  vectors  which  define  the  axes  of  material  ortfiotropy,  then 
their  directions  are  as  shown  below: 


Node  numbers  NI,NJ,NK  are  only  used  to  locate  point  i,j,k, 
respectively,  and  any  convenient  nodes  may  be  used. 
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C.  Distributed  Surface  Load  Data 


NDLS  pairs  of  lines  are  to  be  input  in  this  section  in  order  of 
increasing  set  number  (N).  These  data  describe  surface  loads  acting 
on  element  faces  and  may  be  prescribed  directly  in  terms  of  face 
corner  node  pressures  or  indirectly  by  means  of  a  hydrostatic  pressure 
field. 


Surface  Load  Control  Line  (315) 


Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

N 

Load  set  identification  number; 

1  S  N  <.  NDLS 

(2) 

6-10 

NFACE 

Element  face  number  on  which  this 
distributed  load  is  acting; 

1  S  NFACE  £  6 

(3) 

11  -  15 

LT 

Load  type  code; 

EQ.  1;  prescribed  normal  pressure 
intensities 

EQ.  2;  hydrostatically  varying 
pressure  field 
EQ.  0;  default  set  tc  "l" 

NOTES: 

(1)  The  surface  load  data  sets  established  in  this  section  are 
assigned  to  the  elements  in  Section  F. 

(2)  Hexahedra  have  six  quadrilateral  faces  each  uniquely  described  by 
four  node  numbers  at  the  corners  of  the  face.  The  face  number  convention 
established  for  elements  is  given  in  the  Table  below. 


Table  8.1.  Corner  Node  Numbers  for  the  Solid  Element  Faces 


FACE  NATURAL  CORNER  NODE  NUMBERS 

NUMBER  COORDINATES  Nj  N2  N3  N4 


1  (+1,  s,  t) 

2  (-1,  8,  t) 

3  (  r,+l,  t) 

4  (  r, -1,  t) 

5  (  r,  s,+l) 

6  (  r,  s,-l) 


14  8  5 

2  3  7  6 

15  6  2 

4  8  7  3 

12  3  4 

5  6  7  8 


(3)  Two  types  of  surface  pressure  loads  may  be  applied  to  faces  of 
the  elements.  If  LT  equals  0  (or  1),  a  normal  pressure  distribution  is 
prescribed  directly  by  means  of  pressure  intensities  at  the  face  corner 
nodes.  If  LT  equals  2,  the  face  is  exposed  to  hydrostatic  pressure  due  to 
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Noj 

ml  Pressure  Data 

(4F10.0)  (LT  equals  1,  only) 

Notes 

Columns 

Variables 

Entry 

(1) 

1 

-  10 

PI 

Pressure  at  face  node  N^ 

(2) 

11 

-  20 

P2 

Pressure  at  face  node  N. 
EQ.  0;  default  set  to  "Pi" 

21 

-  30 

P3 

Pressure  at  face  node  N,; 
EQ.  0;  default  set  to  "Pi" 

31 

o 

■ 

P4 

Pressure  at  face  node  N. ; 

EQ.  0;  default  set  to  "Pi" 

NOTES: 

(1)  The  pressure  distribution  acting  on  an  element  face  is  defined  by 
specifying  intensities  P1,P2,P3,P4  at  the  face  corner  nodes  as  shown  below: 


The  tec.*  corner  node  numbers  are  given  on  the  element  data  line 
(Section  F)  and  positive  pressure  tends  to  compress  the  volume  of  the 
element. 


The  variation  of  pressure  over  the  element  face,  p(a,b),  is  given 

as: 


p(a,b) 


where 


i  «  Plxhj 

+  P2xh2  + 

P3xh3  + 

’  (1/4) 

(1+a) 

(1+b) 

*  (1/4) 

(1-a) 

( 1+b) 

«  (1/4) 

(1-a) 

(1-b) 

*  (1/4) 

(1+a) 

(1-b) 

P4xh. 


in  quadrilateral  natural  face  coordinates  (a,b) 


(2)  If  any  of  the  entries  P2,P3,P4  are  omitted,  these  values  are  reset 
to  the  value  of  PI;  i.e.,  for  a  uniformly  distributed  pressure  (p),  we  have 
PI  equals  p  and  columns  11*40  blank.  If  t2  is  zero  specify  a  small  number. 
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Notes 

Columns 

Variables 

Entrv 

(1) 

1  -  10 

GAMMA 

Height  density  of  the  fluid,  Y;  GT.  0 

(2) 

11  -  20 

XS 

X-ordinate  of  point  s  in  the  free 
surface  of  the  fluid 

21  -  30 

YS 

Y-ordinate  of  point  s  in  the  free 
surface  of  the  fluid 

31  -  40 

ZS 

Z -ordinate  of  point  s  in  the  free 
surface  of  the  fluid 

41-50 

XN 

X-ordinated  of  a  point  n  on  the  normal 
to  the  fluid  surface 

51  -  60 

YN 

Y-ordinate  of  a  point  n  on  the  normal 
to  the  fluid  surface 

61  -  70 

ZN 

Z -ordinate  of  a  point  n  on  the  normal 
to  the  fluid  surface 

NOTES: 

(1)  GAMMA  is  the  weight  density  (i.e.,  units  of  force  per  unit  of 
fluid  volume)  of  the  fluid  in  contact  with  element  face  number  NFACE. 

(2)  Point  "s"  is  any  point  in  the  free  surface  of  the  fluid,  and 
point  "n"  is  located  such  that  the  direction  from  s  to  n  is  normal  to  the 
free  surface  and  is  positive  with  increasing  depth. 


Hydrostatic  pressure  in  contact  with  an  element  face  causes  element 
compression;  i.e.,  pressure  resultant  acts  toward  the  element  centroid. 
Nodes  located  above  the  fluid  surface  are  automatically  assigned  zero  pres¬ 
sure  intensities  if  an  element  face  is  not  (or  only  partially)  submerged  in 
the  fluid. 
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D.  Stress  Ontput  Request  Location  Sets  (715) 

If  NOPSET  is  zero  on  the  Control  Line,  skip  this  section,  and 
global  stresses  will  be  computed  and  output  at  the  elements 


centroid  only. 

Otherwise, 

NOPSET  lines 

must  be  input  as 

follows : 

Notes  Columns 

Variables 

Entrv 

(1)  1 

- 

5 

LOCI 

Location 

number 

of 

output 

point 

1 

6 

- 

10 

L0C2 

Location 

number 

of 

output 

point 

2 

11 

- 

15 

L0C3 

Location 

number 

of 

output 

point 

3 

16 

- 

20 

L0C4 

Location 

number 

of 

output 

point 

4 

21 

- 

25 

L0C5 

Location 

number 

of 

output 

point 

5 

26 

- 

30 

L0C6 

Location 

number 

of 

output 

point 

6 

31 

- 

35 

L0C7 

Location 

number 

of 

output 

point 

7 

• 

:} 

Up  to  27 

points 

NOTES: 

(1)  27  element  locations  are  assigned  numbers  as  shown  in  the  Figure 

below.  Locations  1  to  21  correspond  to  node  numbers  1  to  21,  respectively. 
Locations  22  to  27  are  element  face  centroids.  The  first  zero  (or  blank) 
entry  on  a  location  line  terminates  reading  of  location  numbers  for  the 
output  set;  hence,  fewer  than  seven  locations  can  be  requested  in  an  output 
set.  Location  numbers  must  be  input  in  order  of  increasing  magnitude; 
i.e.,  L0C2  is  greater  than  LOCI,  L0C3  is  greater  than  L0C2,  etc.  In 
dynamic  analysis  FACE  1, 

FACE  2, . . . ,  FACE  6  correspond  to  output  locations 
22, 23,..., 27  respectively.  (See  Table  8.1). 

E.  Element  Load  Case  Multipliers 

Five  (5)  lines  must  be  input  in  this  section  specifying  the 
fraction  of  gravity  (X,Y,Z),  the  fraction  of  thermal  loads 
and  the  fraction  of  pressure  loads  to  be  added  to  each  of  the 
element  loading  combination  (A,B, — ).  Load  a  case 
multiplier  data  affect  static  analysis  calculations  only. 


X -Direct ion 

Gravity, 

Line  1  (4F10.0) 

Notes 

Columns 

Variable 

Entrv 

(1) 

1  -  10 

GXA 

Fraction  of  X-direction  gravity  to  be 
applied  in  element  load  case  A 

31  -  40 

GXD 

Fraction  of  X-direction  gravity  to  be 
applied  in  element  load  case  D 

64 
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Y -Direct ion  Gravity,  Line  2  (4F10.0) 
Z -Direction  Gravity,  Line  3  (4F10.0) 
Thermal  loads.  Line  4  (4F10.0) 


Notes 

Columns 

Variables 

Entry 

(2) 

1 

- 

10 

TA 

Fraction  of  thermal 

loads 

to  be 

applied 

in  element  load  case 

A 

11 

- 

20 

TB 

Fraction  of  thermal 

loads 

to  be 

applied 

in  element  load  case 

B 

21 

- 

20 

TC 

Fraction  of  thermal 

loads 

to  be 

applied 

in  element  load  case 

C 

31 

- 

40 

TD 

Fraction  of  thermal 

loads 

to  be 

applied 

in  element  load  case 

D 

Pressure 

loads.  Line  5 

(4F10.0) 

Notes 

Columns 

Variables 

Entry 

(3) 

1 

- 

10 

PA 

Fraction  of  pressure 

loads 

to  be 

applied 

in  element  load  case 

A 

11 

- 

20 

PB 

Fraction  of  pressure 

loads 

to  be 

applied 

in  element  load  case 

B 

21 

- 

20 

PC 

Fraction  of  pressure 

loads 

to  be 

applied 

in  element  load  case 

C 

31 

- 

40 

PD 

Fraction  of  pressure 

loads 

to  be 

applied 

in  element  load  cose  D 

NOTES: 

(1)  Gravity  loads  on  the  structure  due  to  static  body  forces  are 
computed  from  the  weight  density  of  element  materials  and  the  element 
geometry.  These  loads  are  assigned  to  the  element  load  combinations  by 
aieans  of  the  entries  on  Lines  1,2  and  3  for  forces  in  the  X,Y,Z  directions, 
respectively. 

(2)  Thermal  loads  are  computed  knowing  the  node  temperatures  input  in 
the  Nodal  Point  Data  Section,  the  stress  free  reference  temperature  (T  ) 
input  in  Section  F  and  the  element's  material  properties  and  node 
coordinates.  The  temperature  distribution  within  the  element  is  described 
using  the  same  Interpolation  functions  which  describe  the  variation  of 
displacements  within  the  element. 

(3)  Pressure  loads  are  first  assigned  to  element  load  cases  (A,B,...) 
by  means  of  the  entries  (scale  factors)  on  Line  5,  and  the  distributed  load 
sets  which  were  input  in  Section  C  are  then  applied  to  the  elements 

individually  for  cases  (A,B, _ )  by  means  of  load  set  references  given  in 

Section  F. 

F.  Element  Lines 

Two  line  (if  MAXNOD  equals  8)  or  three  lines  (if  MAXNOD.GT.  8) 
must  be  prepared  for  each  element  that  appears  in  the  input, 
and  the  format  for  these  lines  is  as  follows: 
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Eleaent  Line  1  (6I5,F10. ,415,412) 


Notes 

Columns 

Variable 

Entry 

(1) 

1 

-  5 

M 

Element  number; 

1  £  M  £  NS0L21 

(2) 

6 

-  10 

NDIS 

Number  of  nodes  to  be  used  in  describing 
the  element's  displacement  field; 

EQ.  0;  default  set  to  "MAXN0D" 

(3) 

11 

-  15 

NXYZ 

Number  of  nodes  to  be  used  in  the  descrip¬ 

tion  of  element  geometry; 

EQ.  0;  default  set  to  "NDIS" 

EQ.  NDIS  isoparametric  element 

LT.NDIS  subparametric  element 

16 

-  20 

NMAT 

Material  identification  number 

GE.  1  AND  LE.NUMMAT 

(4) 

21 

-  35 

MAXES 

Identification  number  of  the  material 
axis  orientation  set; 

1  £  MAXES  £  N0RTH0 

EQ.  0;  material  axes  default  to  the 

global  X,  Y,  Z  system 

(5) 

26 

-  30 

I0P 

Identification  number  of  the  stress  output 
location  set; 

1  £  I0P  £  N0PSET 

EQ.  0;  centroid  output  only 

31 

-  40 

TZ 

Stress  free  reference  temperature,  T 

Node  number  increment  for  element  da?a 

(6) 

41 

-  45 

KG 

generation; 

EQ.  0;  default  set  to  "l” 

46 

-  50 

NRSINT 

Integration  order  for  natural  coordinate 
(r,s)  direction; 

EQ.  0;  default  set  to  "INTRS" 

51 

-  55 

NTINT 

Integration  order  for  natural  coordinate 
(t)  direction: 

EQ.  0;  default  set  to  "INTT" 

(7) 

56 

-  60 

IREUSE 

Flag  indicating  that  the  stiffness  and 
mass  matrices  for  this  element  are  the 
same  as  those  for  the  preceding  element; 
EQ.  0;  no 

EQ.  1;  yes 

(8) 

61 

-  62 

LSA 

Pressure  set  for  element  load  case  A 

63 

-  64 

LSB 

Pressure  set  for  element  load  case  B 

65 

-  66 

LSC 

Pressure  set  for  element  load  case  C 

67 

-  68 

LSD 

Pressure  set  for  element  load  case  D 

LSD  £  NDLS  defined  in  Section  A 
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Notes 


(9) 


(10) 


Note 


Element  Line  2  (1615) 


Columns 

Variables 

Entry 

1  - 

5 

Node 

1 

number 

6  - 

10 

Node 

2 

number 

11  - 

15 

Node 

3 

number 

16  - 

20 

Node 

4 

number 

21  - 

25 

Node 

5 

number 

26  - 

30 

Node 

6 

number 

31  - 

35 

Node 

7 

number 

36  - 

40 

Node 

8 

number 

41  - 

45 

Node 

9 

number 

46  - 

50 

Node 

10 

number 

51  - 

55 

Node 

11 

number 

56  - 

60 

Node 

12 

number 

61  - 

65 

Node 

13 

number 

66  - 

70 

Node 

14 

number 

71  - 

75 

Node 

15 

number 

76  - 

80 

node 

16 

number 

Slew 

snt 

Line  3  (515) 

(required 

if 

MAXNOD. GT.  8) 

Columns 

Variables 

Entry 

1  - 

5 

Node 

17 

number 

6  - 

10 

Node 

18 

number 

11  - 

15 

Node 

19 

number 

16  - 

20 

Node 

20 

number 

21  - 

25 

Node 

21 

number 

NOTES: 

(1)  Element  lines  must  be  input  in  ascending  element  number  order 
beginning  with  "l"  and  ending  with  "NS0L21".  Repetition  of  element  numbers 
is  not  permitted,  but  element  lines  may  be  omitted,  and  missing  element 
data  are  generated  according  to  the  procedure  described  in  note  (7). 


(2)  NDIS  is  a  count  of  the  node  numbers  actually  posted  on  Lines  2 
and  3  which  must  immediately  follow  Line  1.  NDIS  must  be  at  least  eight 
(8),  but  must  be  less  than  or  equal  to  the  limit  (MAXNOD)  which  was  given 
on  the  Control  Line,  Section  A.  Element  displacements  are  assigned  at  the 
NDIS  nonzero  nodes,  and  thus,  the  order  of  the  element  matrices  is  three 
(i.e.,  translations  X,Y,Z)  times  NDIS.  The  eight  corner  nodes  of  the 
hexahedron  must  be  input,  but  nodes  9  to  21  are  optional,  and  any  or  all  of 
these  optional  nodes  may  be  used  to  describe  the  element ' s  displacement 
field. 
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(3)  When  element  edges  are  straight  it  is  unnecessary  computationally 
to  include  side  nodes  in  the  numerical  evaluation  of  coordinate  derivatives, 
the  Jacobian  matrix,  etc.,  and  since  regular  element  shapes  are  common,  an 
option  has  been  included  to  use  fewer  nodes  in  these  geometric  calculations 
than  are  used  to  describe  element  displacements.  The  first  NXYZ  nonzero 
nodes  posted  on  Lines  2  and  3  are  used  to  evaluate  those  parameters  which 
pertain  to  element  geometry  only.  NXYZ  must  be  at  least  eight  (8),  and  if 
omitted  is  reset  to  NDIS.  A  common  application  might  be  a  20  node  element 
(i.e. ,  NDIS  equals  20)  with  straight  edges  in  which  case  NXYZ  would  be 
entered  as  "8". 

(4)  AXES  (unless  omitted)  refers  to  one  of  the  material  axes  set 
defined  in  Section  B.  If  omitted,  the  material  (NMAT)  orientation  is  such 
that  the  (X^.X^jX^)  axes  coincide  with  the  (X,Y,Z)  axes,  respectively. 

(5)  IOP  (unless  omitted)  refers  to  one  of  the  output  location  sets 
given  in  Section  D.  If  IOP  equals  0,  stress  output  is  quoted  at  the 
element  centroid  only.  Stress  output  at  a  point  consists  of  three  normal 
and  three  shear  components  referenced  to  the  global  (X,Y,Z)  axes. 

(6)  When  element  lines  are  omitted,  element  data  are  generated 
automatically  as  follows: 

(a)  all  data  on  Line  1  for  generated  elements  is  taken  to  be  the 
same  as  that  given  on  the  first  element  line  in  the  sequence; 

(b)  nonzero  node  numbers  (given  on  Lines  2  and  3  for  the  first 
element)  are  incremented  by  the  value  "KG"  (which  is  given 
on  Line  1  of  the  first  element)  as  element  generation 
progresses;  zero  (or  blank)  node  number  entries  are 
generated  as  zeroes . 

(c)  The  last  element  cannot  be  generated. 

(7)  The  flag  IREUSE  allows  the  program  to  bypass  stiffness  and  mass 
matrix  calculations  providing  the  current  element  is  identical  to  the  pre¬ 
ceding  element;  i.e.,  the  preceding  and  current  elements  are  identical 
except  for  the  rigid  body  translation.  If  IREUSE  equals  0,  new  matrices 
are  computed  for  the  current  element.  If  IREUSE  equals  1  it  is  also 
assumed  that  the  node  temperatures  of  the  element  (for  calculation  of 
thermal  loads)  are  the  same  as  those  of  the  preceding  element. 

(8)  Pressure  loads  are  assigned  (i.e.,  applied)  to  the  element  by 
means  of  load  set  references  in  columns  61-62  for  combination  A,  columns 
63-64  for  B,  etc.  A  zero  entry  means  that  no  pressure  acts  on  the  element 
for  that  particular  element  load  combination. 

(9)  The  first  eight  node  numbers  establish  the  corners  or  vertices  of 
a  general  hexahedron  and  must  be  all  nonzero,  (see  the  Figures  associated 
with  Section  A).  Node  numbers  must  be  input  in  the  sequence  indicated 
otherwise  volume  and  surface  area  integrations  will  be  indefinite. 
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(10)  The  number  of  lines  required  as  input  for  each  element  depends 
on  the  variable  MAXNOD.  For  the  case  of  MAXNOD  equals  8,  only  Line  2  is 
required.  If  MAXNOD  .GT.  8,  Lines  2  and  3  are  required  for  all  elements. 

Nodes  9  to  21  are  optional,  and  only  those  nodes  actually  used  to 
describe  the  element  are  input.  The  program  will  read  all  21  entries  if 
MAXNOD  was  given  as  9  or  greater,  but  only  NDIS  nonzero  values  are  expected 
to  be  read  on  Lines  2  and  3.  If  for  example  one  element  is  described  by  10 
nodes,  then  columns  1-40  on  Line  2  would  be  the  eight  corner  node  numbers, 
and  the  remaining  two  node  numbers  would  be  posted  somewhere  on  Lines  2  and 
3. 
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THREE-DIMENSIONAL  STRAIGHT  OR  CURVED  PIPE  ELEMENTS,  ELEMENT  TYPE  12 


Cylindrical  pipe  elements  are  identified  by  the  number  twelve  (12). 
Axial  and  shear  forces,  torque  and  bending  moments  are  calculated  for  each 
member.  Gravity  loading  in  the  global  (X,Y,Z)  directions,  uniform 
temperature  changes  (computed  from  input  nodal  temperatures),  and 
extensional  effects  due  to  internal  pressure  form  the  basic  member  loading 
conditions.  Pipe  element  input  is  described  by  the  following  sequence  of 
lines: 


A. 

Control  Line 

(1415) 

Notes 

Columns 

Variables 

Entry 

4-5 

Enter  the  number  "12" 

(1) 

6-10 

NPIPE 

Number  of  pipe  elements 

11  -  15 

NUMMAT 

Number  of  material  sets 

16  -  20 

MAXTP 

Maximum  number  of  temperature  points 
used  in  the  table  for  any  material 

GE.  1;  at  least  one  point 

21  -  25 

NSECT 

Number  of  section  property  sets;  GE.  1 

(2) 

26  -  30 

NBRP 

Number  of  branch  point  nodes  at  which 

output  is  required; 

EQ.  0;  no  branch  point  output  is 


produced 


31  -  35 

MAXTAN 

Maximum  number  of  tangent  elements 
common  to  any  one  branch  point  node; 

EQ.  0;  default  set  to  "4" 

36  -  40 

NPAR  (8) 

Blank 

41  -  45 

NPAR  (9) 

Tangent  stiffness  load  matrix  dump  flag 
EQ.  1;  Print 

EQ.  0;  Suppress  printing 

46  -  50 

NPAR(10) 

Bend  stiffness  load  matrix  dump  flag 

EQ.  1;  Print 

EQ.  0;  Suppress  printing 

51  -  55 

NPAR(ll) 

Element  parameters  dump  flag 

EQ.  1;  Print 

EQ.  0;  Suppress  printing 

NOTES: 

(1)  The  number  of  pipe  elements  ("NPIPE")  counts  both  tangent  and 
bend  geometries,  and  both  the  material  and  section  property  tables  can 
reference  either  the  bend  or  tangent  element  types. 

(2)  A  branch  point  is  defined  as  a  nodal  locations  where  at  least 
three  tangent  pipe  element  connect.  The  two  input  parameters  "NBRP"  and 
"MAXTAN"  reserve  storage  for  an  index  array  created  during  the  processing 
of  pipe  element  data;  posting  a  larger  number  of  maximum  common  tangents 
than  actually  exist  is  not  considered  a  fatal  error.  Branch  point  data  is 
read  if  requested,  but  not  currently  used;  i.e.  to  be  used  in  future 
program  versions. 
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B.  Material  Property  Lines 


Temperature -dependent  Young's  modules  (E),  Poisson's  ratio 
(v)  and  thermal  expansion  coefficient  (a)  are  allowed.  If 
more  than  one  temperature  point  is  input  for  a  material 
table,  then  the  program  selects  properties  using  linear 
interpolation  between  input  temperature  values.  The 
temperature  used  for  property  selection  is  the  average 
element  temperature  which  is  denoted  as  Tfl: 

T  =  (T.  +  T . )/2 
a  i  j 

where  T  and  T,  are  input  nodal  temperatures  for  ends  "i"  and  "j"  of  the 
pipe.  For  each  different  material,  the  following  set  of  lines  must  be 
input : 


Material  Identification  Line  (215, 6A6) 


Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

M 

Material  identification  number; 

1  £  M  £  NUMMAT 

6-10 

11  -  46 

NT 

Number  of  different  temperatures  at 
which  properties  are  given; 

EQ.  0;  one  temperature  point  is 
assumed  to  be  input 

Material  description  used  to  label 
the  output  for  this  material 

NOTES: 

CD 

Material 

ident  i  f  icat  ion 

number  must  be  input  between  one  and  the 

total  number  of  materials  specified  (NUMMAT) 

Material  Lines  (4F10.0) 

Notes 

Columns 

Variables 

Entry 

(1) 

1-10 

T(N) 

Temperature,  T 

11  -  20 

E(N) 

Young's  modulus,  E 

21  -  30 

XNU(N) 

Poisson's  ratio,  v 

31  -  40 

ALP(N) 

Thermal  expansion  coefficient,  aQ 

NOTES: 

(1) 

Supply  one  line  for  each  temperature  point  in  the  material  table; 

at  least  one  line  is  required.  Temperatures  must  be  input  in  increasing 
(algebraic)  order.  If  two  or  more  points  are  used,  care  must  be  taken  to 
Insure  that  the  table  covers  the  expected  range  of  average  temperatures 
existing  in  the  elements  to  which  the  material  table  is  assigned. 
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Section  Property  Lines  (I5,5F10.0,3A6) 


Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

N 

Section  property  identification  number; 
MNS  NSECT 

(2) 

6-15 

16  -  25 

26  -  35 

Outside  diameter  of  the  pipe,  dQ 

Pipe  wall  thickness,  t 

Shape  factor  for  shear  distortion,  a 

(3) 

36-45 

Weight  per  unit  length  of  section, 

(4) 

46  -  55 

56  -  73 

Mass  per  unit  length  of  section, 

Section  description  (used  to  label  the 
output ) 

NOTES: 

(1)  Section  property  identification  number  must  be  input  in  an  ascend 
ing  sequence  beginning  with  one  and  ending  with  the  total  number  of  section 
specified  (NSECT). 

(2)  Assuming  that  (y,z)  are  the  section  axes  and  that  the  x-axis  is 
normal  to  the  section,  the  properties  for  the  section  are  computed  from  the 
input  parameters  [ d^,  t  and  o^]  as  follows: 

(a)  inner  and  outer  pipe  radii; 


(b)  cross -sectional  area  (axial  deformations); 


(c)  principal  moments  of  inertia  (bending); 

Iy  *  (»/4)  (rQ  '  **) 


(d)  polar  moment  of  inertia  (torsion); 


(e)  effective  shear  areas  (shear  distortions); 
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Note  that  the  shape  factor  for  shear  distortion  (a  )  may  be  input 
directly.  If  the  entry  is  omitted,  the  shape  factor  is  computed  using  the 
equation: 

«,  *  W3)  (r l  -  r*)/|(r*  +  r*)  (r<)  -  r^] 

=  2.0 

An  input  value  for  a  greater  than  one  hundred  (100.)  causes  the 
program  to  neglect  shear  distortions  entirely.  If  used,  the  same  shape 
factor  is  applied  to  both  in  and  out-of -plane  shear  distortions. 

(3)  The  weight  per  unit  length  of  section  (Y^)  is  used  to  compute 
gravity  loadings  on  the  elements.  Fixed  end  shears,  moments,  torques,  etc. 
are  computed  automatically  and  applied  as  equivalent  nodal  loads.  These 
forces  will  not  act  on  the  structure  unless  first  assigned  to  one  of  the 
element  load  cases  (A,B,C,D)  in  Section  D  below. 

(4)  The  mass  per  unit  length  is  only  used  to  form  the  lumped  mass 
matrix  for  a  dynamic  analysis  case.  If  no  entry  is  input,  then  the  program 
will  redefine  the  mass  density  from  the  weight  density  using: 

Pl  =  Yj/386.4 

Either  a  nonzero  weight  density  or  mass  density  will  cause  the  program 
to  assign  masses  to  all  pipe  element  nodes. 

C.  Branch  Point  Node  Numbers 


If  the  number  output  branch  point  nodes  has  been  omitted  from  the 
control  line  (i.e.,  columns  26-30  are  blank),  skip  this  section  of 
input,  and  no  branch  point  data  will  be  read.  Otherwise,  supply  node 
numbers  for  a  total  number  of  branch  points  requested  on  control  line, 
ten  nodes  per  line: 


First  Line  (1015) 

Notes  Columns  Variables  Entry 


(1) 


1-5 

6-10 


Node  number  at  branch  point  1 
Node  number  at  branch  point  2 


45  -  50  Node  number  at  branch  point  10 

Second  Line  (1015)  --  if  required 
Notes  Columns  Variable  Entry 

1-5  Node  number  at  branch  point  11 

NOTES: 

(1)  A  node  does  not  define  a  branch  point  unless  at  least  three 
tangent  elements  are  common  to  the  node.  Branch  point  output  is  only 
produced  for  static  analysis  cases. 
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D.  Element  Load  Case  Multipliers 

Five  lines  must  be  input  in  this  section  specifying  the  fraction  of 
gravity  (in  each  of  the  X,Y,Z  coordinate  directions),  the  fraction  of 
thermal  loading  and  the  fraction  of  internal  pipe  pressure  loading  the 
be  added  to  each  of  four  possible  element  loading  combinations 
(A,B,C,D). 


X-Direction  Gravity,  Line  1  (4F10.0) 


Notes 

Columns 

Variables  Entry 

(1) 

1  -  10 

Fraction  of  X-direction  gravity  to  be 
applied  in  element  lo*d  case  A 

11  -  20 

Fraction  of  X-direction  gravity  to  be 
applied  in  element  load  case  B 

21  -  30 

Fraction  of  X-direction  gravity  to  be 
applied  in  element  load  case  C 

31  -  40 

Fraction  of  X-direction  gravity  to  be 
applied  in  element  load  case  D 

Y -Direction  Gravity,  Line  2  (4F10.0) 

Z -Direction  Gravity,  Line  3  (4F10.0) 

Thermal 

loads.  Line  4  (4F10.0) 

Notes 

Columns 

Variables  Entry 

(2) 

1-10 

Fraction  of  thermal  loading  to  be 
applied  in  element  load  case  A 

11  -  20 

Fraction  of  thermal  loading  to  be 
applied  in  element  load  case  B 

21  -  30 

Fraction  of  thermal  loading  to  be 
applied  in  element  load  case  C 

31  -  40 

Fraction  of  thermal  loading  to  be 
applied  in  element  load  case  D 

Internal 

Pressure,  Line  5  (4F10.0) 

Notes 

Columns 

Variables  Entry 

(3) 

1-10 

Fraction  of  pres sure -induced  loading 
applied  in  element  load  case  A 

11  -  20 

Fraction  of  pres sure -induced  loading 
applied  in  element  load  case  B 

21  -  30 

Fraction  of  pressure- induced  loading 
applied  in  element  load  case  C 

31  -  40 

Fraction  of  pressure- induced  loading 

applied  in  element  load  case  D 
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VERTICAL  TANGENT 


LOCAL  COORDINATE  SYSTEMS  FOR  PIPE  ELEMENTS 
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NOTES: 

(1)  No  gravity  loads  will  be  produced  if  the  weight  per  unit  length 
was  input  as  zero  on  all  section  property  lines.  Otherwise,  a  multiplier 
of  1.0  input  for  an  element  load  case  means  that  100  percent  of  deadweight 
will  be  assigned  to  that  load  combination. 

(2)  No  thermal  loading  will  result  if  the  coefficient  of  thermal 

expansion  has  been  omitted  from  all  the  material  lines.  Otherwise,  thermal 
loads  are  computed  for  each  element  using  the  AT  between  the  average  element 
temperature  (T  )  and  the  stress -free  temperature  (T  )  given  with  each  pipe 
element  line  (Section  E,  below).  ° 

(3)  Element  distortions  are  computed  for  each  element  due  to  internal 
pressure,  and  these  loads  are  combined  into  element  load  cases  by  means  of 
appropriate  nonzero  entries  in  Line  5. 

Gravity,  thermal  or  pressure  induced  loads  cannot  act  on  the  structure 
unless  first  combined  in  one  or  more  of  the  element  load  sets  (A,B,C,D). 

Once  defined,  element  loads  are  assigned  (via  scale  factors)  to  the 
structure  load  cases  by  means  of  Element  Load  Multipliers  given  in  Section 
D.  An  element  load  case  combination  may  be  used  multiple  times  when 
defining  the  various  structure  loading  conditions. 


E. 

Pipe  Element  Lines 

Line  Type 

1 

Notes 

Columns 

Variables 

Entry 

(1) 

1 

- 

4 

N 

Pipe  element  number; 

5 

GE.  1  and  LE.NPIPE 

Geometric  type  code: 

"T"  (or  blank);  tangent  section 
"B"  ;  bend  (circular)  section 

6 

- 

10 

I 

Node  I  number 

11 

- 

15 

J 

Node  J  number 

16 

— 

20 

MAT 

Material  identification  number; 

GE.  1  and  LE.NUMMAT 

21 

— 

25 

I  SECT 

Section  property  identification  number 
GE.  1  and  LE.NSECT 

(2) 

26 

- 

35 

Stress-free  temperature,  T 

(3) 

36 

- 

45 

Internal  pressure,  p 

(4) 

46 

' 

55 

Positive  projection  of  a  local  y- 
vector  on  the  global  X-axis;  A(yX) 

56 

65 

Positive  projection  of  a  local  y- 
vector  on  the  global  Y-axis;  A(yY) 

66 

75 

Positive  projection  of  a  local  y- 
vector  on  the  global  Z-axis  A  (yZ) 

(5) 

76 

- 

80 

KG 

Node  number  increment  for  tangent 

element  generation; 

EQ.  0;  default  set  to  ”l" 
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NOTES: 

(1)  Line  type  1  is  used  for  both  tangent  and  bend  elements;  a  second 
line  (line  type  2,  below)  must  be  input  immediately  following  line  type  1 
if  the  pipe  element  is  a  bend  (i.e.,  "B"  in  column  5).  Note  that  element 
lines  must  be  input  in  ascending  sequence  beginning  with  one  and  ending 
with  the  total  number  of  pipe  elements.  If  tangent  elements  are  omitted, 
generation  of  the  intermediate  elements  will  occur;  the  generation 
algorithm  is  described  below.  An  attempt  to  generate  bend  type  elements  is 
considered  an  error. 

(2)  The  stress-free  temperature,  T  ,  is  subtracted  from  the  average 
element  temperature,  T  ,  to  compute  the  uniform  temperature  difference 
acting  on  the  element: 


AT  =  T  -  T 
a  o 

The  entire  element  is  assumed  to  be  at  this  uniform  value  of 
temperature  difference. 

(3)  The  value  of  pressure  is  used  to  compute  a  set  of  self-equilibrating 
joint  forces  arising  from  member  distortions  due  to  pressurization;  i.e., 
the  mechanical  equivalent  of  thermal  loads.  For  bend  elements,  the  pressure 
is  also  used  to  compute  the  bend  flexibility  factor,  k  .  The  curved  pipe 
subjected  to  bending  is  more  flexible  than  elementary  Beam  theory  would 
predict.  The  ratio  of  "actual"  flexibility  to  that  predicted  by  beam  theory 

is  denoted  by  k  ,  where 
P 

kp  -  ( 1. 65/h)/[ 1  +  (6p/Eh)  (R/t)4/3J  *  1 

in  which 

h  =  tR/r2 
r  -  (dQ  -  t)/2 

and 


t  -  pipe  wall  thickness 
R  =  radius  of  the  circular  bend 
r  -  mean  radius  of  the  pipe  cross  section 
d  =  outside  diameter  of  the  pipe 
£  =  Young's  modulus 
p  =  internal  pressure 

The  flexibility  factor  is  computed  and  applied  to  all  bend 
elements;  pressure  stiffening  is  neglected  if  the  entry  for 
internal  pressure  ("p")  is  omitted. 
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(4)  The  global  projections  of  the  local  y-axis  for  a  tangent  member 
may  be  omitted  (columns  46-75  are  blank);  for  this  case,  the  following 
convention  for  the  local  system  is  assumed: 

(a)  tangents  parallel  to  the  global  Y-axis  (vertical  axis) 
have  their  local  y-axes  directed  parallel  to  and  in  the 
same  directions  as  the  global  Z-axis 

(b)  tangents  not  parallel  to  the  global  Y-axis  have  their 
local  Y-axes  contained  in  a  vertical  (global)  plane  such 
that  local  y  projects  positively  on  the  positive  global 
Y-axis. 

For  bend  elements,  the  global  projections  of  the  local  y-axis  are 
not  used;  instead,  the  local  axis  convention  is  defined  as  follows: 

(a)  the  local  y-axis  is  directed  positively  toward  and 
intersects  the  center  of  curvature  of  the  bend  (i.e., 
radius  vector); 

(b)  the  local  x-axis  is  tangent  to  the  arc  of  the  bend  and  is 
directed  positively  from  node  I  to  node  J. 

Note  that  for  all  elements,  the  local  x,  y,  z  system  is  a  right-handed 
set  (see  Figure  on  page  76). 

(5)  If  a  tangent  element  sequence  exists  such  that  each  element 
number  (NE^)  is  one  greater  than  the  previous  number  (NE^_^);  i.e., 

NBi  =  NE  i_i  +  i 

only  the  element  line  for  the  first  tangent  in  the  series  need  be  input. 

The  node  numbers  for  the  missing  tangents  are  computed  using  the  formulae: 

NIa  =  NI1_1  +  KG 

NJi  =  NJ11  +  KG 

where  "KG"  is  the  node  number  increment  input  in  columns  76-80  for  the 
first  element  in  the  series,  and  the 

(a)  material  identification  number 

(b)  section  property  identification  number 

(c)  stress-free  temperature 

(d)  internal  pressure 

(e)  Y-axis  global  projections 

for  each  tangent  in  the  generation  sequence  are  taken  to  be  the  same  as 
those  input  on  the  first  line  in  the  series.  The  node  number  increment 
("KG")  is  reset  to  one  if  left  blank  on  the  first  line  in  the  series.  The 
last  (highest)  element  cannot  be  generated;  i.e.,  it  must  be  input. 
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Bend  element  data  cannot  be  generated  because  two  Input  lines  are 
required  for  each  bend.  Also,  the  element  just  prior  to  a  bend  element 
must  appear  on  an  input  line.  Several  bends  may  be  input  in  a  sequence, 
but  each  bend  must  appear  (on  two  lines)  in  the  input  stream. 


Line  Type  2  (F10. 0, . 3X.A2.4F10. 0) 
Notes  Columns  Variables  Entry 


(1)  1  -  10  R 

(2)  14  -  15 


16  -  25 
26  -  35 
36  -  45 
46  -  55 


Radius  of  the  bend  element,  R 
Third  point  type  code:  "TI" 

(or  blank);  third  point  is  the 
tangent  intersection  point  "CC" ; 
third  point  is  the  center  of 
curvature 

X-ordinate  of  the  third  point,  X_ 

Y -ordinate  of  the  third  point, 
Z-ordinate  of  the  third  point,  Z^ 
Fraction  of  wall  thickness  to  be  used 
for  dimensional  tolerance  tests; 

EQ.  0;  default  set  to  "0.1" 


NOTES: 

(1)  The  radius  of  the  bend  ("R")  must  be  input  regardless  of  the 
method  ("Tl"  or  "CC")  used  to  define  the  third  point  for  the  bend. 


(2)  If  the  tangent  intersection  point  is  used,  the  program  computes  a 
radius  for  the  bend  and  compares  the  computed  value  with  the  input  radius. 
An  error  is  declared  if  the  two  radii  are  different  by  more  than  the 
specified  fraction  (or  multiple)  of  the  section  wall  thickness.  The 
lengths  of  the  two  tangent  lines  (I  to  TI  and  J  to  TI)  are  compared  for 
equality,  and  an  error  will  be  flagged  if  the  two  values  are  discrepant  by 
more  than  the  dimensional  tolerance. 


If  the  center  of  curvature  is  input,  the  distances  from  the  third 
point  to  nodes  I  and  J  are  compared  to  the  input  radius;  discrepancies 
larger  than  the  user  defined  tolerance  are  noted  as  errors. 

This  second  element  line  is  only  to  be  input  for  the  bend  type 
element. 


F.  Element  Stress  Output 

Stress  output  for  pipe  elements  consists  of  forces  and  moments  acting 
in  the  member  cross  sections  at  the  ends  of  each  member  and  at  the 
midpoints  of  the  arcs  in  bend  elements.  Output  quantities  act  on  the 
element  segment  connecting  the  particular  output  station  and  end  i; 
i.e.,  j  to  i,  center  to  i,  or  AX  to  i  (where  AX  ■*  0).  Positive 
force/moment  vectors  are  directed  into  the  positive  local  (x,y,z) 
directions,  as  shown  in  the  accompanying  figure. 
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FORCE  SIGN  CONVENTION  FOR  PIPE  ELEMENT  OUTPUT 
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CONCENTRATED  LOAD/MASS  DATA  (2I5.6F10.4) 


Notes 

Columns 

Variables 

Entry 

(1) 

1  - 

5 

N 

Nodal  point  number 

(2) 

6  - 

10 

L 

Structure  load  case  number; 

6E.  1;  static  analysis 

EQ.  0;  dynamic  analysis 

11  - 

20 

FX(N,L) 

X-direction  force  (or  translational 
mass  coefficient) 

21  - 

30 

FY(N,L) 

Y-direction  force  (or  translational 
mass  coefficient) 

31  - 

40 

FZ(N,L) 

Z-direction  force  (or  translational 
mass  coefficient) 

41  - 

50 

MX(N,L) 

X-axis  moment  (or  rotational  inertia) 

51  - 

60 

MY(N,L) 

Y-axis  moment  (or  rotational  inertia) 

NOTES: 

61  - 

70 

MZ(M,L) 

Z-axis  moment  (or  rotational  inertia) 

(1) 

For  i 

a  static  analysis 

case  (NDYN  =  0),  one  line  is  required  for 

each  nodal  point  ("N")  having  applied  (nonzero)  concentrated  forces  or 
moments.  All  structure  load  cases  must  be  grouped  together  for  the  node 
("N")  before  data  is  entered  for  the  next  (higher)  node  at  which  loads  are 
applied.  Only  the  structure  load  cases  for  which  node  N  is  loaded  need  be 
given,  but  the  structure  load  case  number  ("L")  which  are  referenced  must 
be  supplied  in  ascending  order.  Node  loadings  must  be  defined  (input)  in 
increasing  node  number  order,  but  again,  only  those  nodes  actually  loaded 
are  required  as  input.  The  static  loads  defined  in  this  section  act  on  the 
structure  exactly  as  input  and  are  not  scaled,  factored,  etc.  by  the  element 
load  case  (A,B,C,D)  multipliers  (refer  to  the  following  Element  Load  Multi- 
pliers  section).  Nodal  forces  arising  from  element  loadings  are  combined 
(additively)  with  any  concentrated  loads  given  in  this  section.  Applied 
force/moment  vectors  act  on  the  structure,  positive  in  the  positive  global 
directions.  Only  one  line  is  allowed  per  node  per  load  case. 

For  a  dynamic  analysis  case  (NDYN  =  1,  2,  3  or  4),  structure  load 
cases  have  no  meaning,  but  the  program  expects  to  read  data  in  this  section 
nonetheless.  In  place  of  concentrated  loads,  lumped  mass  coefficients  for 
the  nodal  degrees  of  freedom  may  by  input  for  any  (or  all)  nodes.  The  mass 
matrix  is  automatically  constructed  by  the  program  from  element  geometry 
and  associated  material  densities;  the  mass  coefficients  read  in  the  this 
section  are  combined  (additively)  with  the  existing  element-based  lumped 
mass  matrix.  For  mass  input,  a  node  may  only  be  specified  once,  and  the 
load  case  number  ("L")  must  be  zero  (or  blank). 
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The  program  terminates  reading  loads  (or  mass)  data  when  a  zero  node 
number  "(N")  is  encountered.  For  the  special  case  of  a  static  analysis 
with  no  concentrated  loads  applied,  input  only  one  line  in  this  section. 
Similarly,  a  dynamic  analysis  in  which  the  mass  matrix  is  not  to  be  aug¬ 
mented  by  any  entries  in  this  section  requires  only  one  line  as  input. 

(2)  For  a  static  analysis,  structure  load  case  numbers  range  from  "l" 
to  the  total  number  of  load  cases  requested  on  the  Master  Control  Line 
("LL");  thus,  1  £  L  5  LL  for  NDYN  >  0.  For  a  dynamic  analysis,  only  zero 
is  allowed  (L  =  0  for  NDYN  equal  to  1,  2,  3,  or  4). 
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ELEMENT  LOAD  MULTIPLIERS  (4F10.0) 


Notes 

Columns 

Variable 

Entrv 

(1,2) 

1-10 

EM  (1) 

Multipliers 

for 

element 

load 

case 

A 

11  -  20 

EM  (2) 

Multipliers 

for 

element 

load 

case 

B 

21  -  30 

EM  (3) 

Multipliers 

for 

element 

load 

case 

C 

31  -  40 

EM  (4) 

Multipliers 

for 

element 

load 

case 

D 

END  OF  DATA  SET  FOR  STATIC  ANALYSIS  (NDYN  =  0) 


NOTES: 

(1)  One  line  must  be  included  for  each  static  (NDYN  equals  0)  structure 
load  case  requested  on  the  Master  Control  Line  ("LL").  Each  line 
references  the  respective  load  case  number  in  ascending  order.  The  four 
element  load  sets  (A,B,C,D),  if  created  in  any  of  the  element  data  sets  are 
combined  with  any  concentrated  loads  specified  in  the  Concentrated 
Load/Mass  section  for  the  structure  loads  cases.  For  example,  suppose  an 
analysis  case  calls  for  seven  static  structure  loading  conditions  (i.e.,  LL 
—7),  then  the  program  expects  to  read  seven  lines  in  this  section. 

Further,  suppose  line  number  three  (i.e.,  structure  load  case  three)  in 

this  section  contains  the  entries: 

[EM(1),EM(2),EM(3)EM(4)J  =  { -3. 0,0. 0,2. 0,0.0] 

Structure  load  case  three  will  then  be  constructed  using  100  percent 
of  any  concentrated  loads  specified  in  the  Concentrated  Load/Mass  section 
minus  ( -)  300  percent  of  the  loads  in  element  set  A  plus  (+)  200  percent  of 
the  loads  in  element  set  C.  Load  sets  B  and  D  will  not  be  applied  in 
structure  load  case  3.  Element  load  sets  may  be  referenced  any  number  of 
tines  in  order  to  construct  different  structure  loading  conditions. 
Element-based  loads  (gravity,  thermal,  etc.)  can  only  be  applied  to  the 
structure  by  means  of  the  data  entries  in  this  section. 

Static  analysis  input  is  complete  with  this  section.  If  a  new  problem 
is  not  desired,  provide  two  additional  lines:  one  representing  the  heading 
line  with  "END"  entered  beginning  in  Column  1,  the  second  representing  the 
master  control  card  with  the  number  of  nodes  (NUMNP)  equal  to  zero.  Begin 
a  new  problem  data  set  with  a  new  heading  line. 

(2)  If  this  case  calls  for  one  of  the  dynamic  analysis  options, 
supply  only  one  blank  line  in  this  section.  If  the  job  is  a  dynamic 
restart  case  (NDYN  equals  -2  or  -3),  skip  this  section. 
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DYNAMIC  ANALYSES 


Four  types  of  dynamic  analysis  can  be  performed  by  the  program. 

The  type  of  analysis  is  indicated  by  the  number  "NDYN"  specified  in  columns 
21-25  of  the  Master  Control  Line.  If: 

NDYN  .EQ.  1;  Determination  of  system  mode  shapes  and  frequencies  only 
(complete  Mode  Shape  and  Frequency  section  only) 

NDYN  .EQ.  2;  Dynamic  Response  Analysis  for  arbitrary  time  dependent 

loads  using  mode  superposition  (complete  both  Mode  Shape 
and  Frequency,  and  Response  History  sections) 

NDYN  .EQ.  3;  Response  Spectrum  Analysis  (complete  both  Mode  Shape  and 
Frequency,  and  Response  Spectrum  sections) 

NDYN  .EQ.  4;  Dynamic  Response  Analysis  for  arbitrary  time  dependent 
loads  using  step-by-step  direct  integration  (complete 
Response  History  section  only) 

In  any  given  dynamic  analysis  case  only  one  value  of  NDYN  will  be 
considered.  However,  if  NDYN  equals  2  or  3,  the  program  oiust  first  solve 
the  eigenvalues  problem  for  structure  modes  and  frequencies.  These 
eigenvalues/vectors  are  then  used  as  input  to  either  the  Forces  Response 
Analysis  (NDYN  equals  2)  or  to  the  Response  Spectrum  Analysis  (NDYN  equals 
3).  Hence,  options  1,  2  or  3  all  require  that  the  control  parameters  for 
eigenvalue  extraction  be  supplied. 

In  case  of  a  direct  step-by-step  integration  analysis  (NDYN  equals  4) 
do  not  provide  the  eigenvalue  solution  control  line. 

For  the  special  case  of  dynamic  analysis  restart  (NDYN  equals  -2  or 
-3), data  input  consists  of  the  Heading  Line,  the  Master  Control  Line,  and 
either  the  Response  History  section  (NDYN  =  -2)  or  the  Response  Spectrum 
section  (NDYN  =  -3),  below.  Restarting  is  possible  only  if  a  previous 
eigenvalue  solution  using  the  same  model  was  performed  with  NDYN  equals  1, 
and  the  results  from  this  eigenvalue  solution  were  saved  on  the  restart 
file  (refer  to  dynamic  analysis  restart  instructions). 

Up  to  this  section  the  program  processes  (i.e.,  expects  to  read) 

essentially  the  same  blocks  of  data  for  either  the  static  or  dynamic  analy¬ 

sis  cases;  certain  of  these  preceding  data  lines,  however,  are  read  by  the 
program  but  are  not  used  in  the  dynamic  analysis  phase.  In  general,  the 
purpose  of  the  preceding  data  sections  is  to  provide  information  leading  to 
the  formation  of  the  systems  stiffness  and  mass  matrices  (appropriately 
modified  for  displacement  boundary  conditions).  For  example,  element  load 
sets  (A,  B,  C,  D)  may  be  constructed  as  through  a  static  case  were  to  be 
considered,  but  these  data  are  not  used  in  a  dynamic  analysis;  i.e.,  the 

same  data  file,  through  the  Element  Library  section,  can  be  used  for  either 

type  of  analysis.  The  concept  of  structure  loading  conditions  is  not 
applicable  for  dynamic  analysis,  and  input  for  Concentrated  Load/Mass  Data 
and  Element  Load  Multipliers  must  be  prepared  specially. 
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MODE  SHAPES  AND  FREQUENCIES  (3I5,2F10.0) 


This  section  is  required  for  NDYN  equals  1,  2  or  3. 

A  diagonal  (lumped)  mass  matrix  is  formed  automatically  using  element 
geometry  and  assigned  material  density  or  densities.  The  mass  matrix  so 
defined  contain  only  translational  mass  coefficients  calculated  from 
tributary  element  volumes  common  to  each  node.  Known  rotational  inertias 
must  be  input  for  the  individual  nodal  degrees  of  freedom  in  the 
Concentrated  Load/Mass  section. 

Nonzero  impressed  displacement  (or  rotations)  input  by  means  of  the 
BOUNDARY  element  (type  "7")  are  ignored;  instead  the  component  is  restrained 
against  motion  during  dynamic  motion  of  the  structure. 

The  program  does  not  change  the  order  of  the  system  by  performing  a 
condensation  of  those  nodal  degrees  of  freedom  having  no  (zero)  mass  coef¬ 
ficients  (i.e.,  a  zero  mass  reduction  is  not  performed).  No  distinction  is 
made  between  static  and  dynamic  degrees  of  freedom  (i. e. ,  they  are  identical 
in  sequence,  type  and  total  number) . 


Notes 

Columns 

Variables 

Entry 

(1) 

1  -  5 

IFPF 

Flag  for  printing  intermediate 
matrices,  norms,  etc.  calculated 
during  the  eigenvalue  solution; 

EQ.  0;  do  not  print 

EQ.  1;  print 

(2) 

6-10 

IFSS 

Flag  for  performing  the  STURM  SEQUENCE 
check; 

EQ.  0;  check  to  see  if  eigenvalues 
were  missed 

EQ.  1;  pass  on  the  check 

(3) 

11  -  15 

NITEM 

Maximum  number  of  iterations  allowed 
to  reach  the  convergence  tolerance; 

EQ.  0;  default  set  to  "16" 

(4) 

16  -  25 

RTOL 

Convergence  tolerance  (accuracy)  for 
the  highest  ("NF”)  requested  eigenvalue; 
EQ.  0;  default  set  to  "l.OE-5" 

(5) 

26  -  35 

COFQ 

Cut-off  frequency  (cycles/unit  time) 

EQ.  0;  NF  eigenvalues  will  be  extracted 
GT.  0;  extract  only  those  values 
below  COFQ 

(6) 

36  -  40 

NFO 

Number  of  starting  Iteration  vectors 

to  be  read  from  TAPE10 
END  OF  DATA  SET  FOR  NDYN  =  1 

NOTES: 

(1)  Extra  output  produced  by  the  eigenvalue  solutions  can  be  requested; 
output  produced  by  this  option  can  be  quite  voluminous.  Normal  output  pro¬ 
duced  by  the  program  consists  of  an  ordered  list  of  eigenvalues  followed  by 
the  eigenvectors  for  each  mode.  The  number  of  modes  found  and  printed  is 
specified  by  the  variable  "NF"  given  in  column  16-20  of  the  Master  Control 
Line. 
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(2)  The  program  performs  the  solution  for  eigenvalues/vectors 
using  either  of  two  (2)  district  algorithms; 

(a)  the  DETERMINANT  SEARCH  algorithm  requires  that  the  upper 
triangular  band  of  the  system  stiffness  matrix  fit  into 
memory  (in-core);  i.e.,  one  equation  "block". 

(b)  the  SUBSPACE  ITERATION  algorithm  is  used  if  only  portions 
(fractions)  of  the  system  matrix  can  be  retained  in  memory; 
i.e.,  the  matrix  (even  though  in  band  form)  must  be 
manipulated  in  blocks. 

The  program  will  automatically  select  the  SUBSPACE  ITERATION  procedure 
for  eigenvalue  solution  if  the  model  is  too  large  for  the  in-core 
algorithm. 

The  entries  "IFSS",  "NITEM"  and  "RTOL"  ignored  if  the  program  can  use 
the  DETERMINANT  SEARCH  to  find  eigenvalues.  Whether  or  not  a  model  is  too 
large  for  the  DETERMINANT  SEARCH  depends  on  the  amount  of  core  allocated 
(by  the  programmer  and  not  the  user)  for  array  storage.  The  program 
variable  "MTOT"  equals  the  amount  of  working  storage  available. 

Define: 


MBAND  =  maximum  equation  bandwidth  (coefficients) 

*  (maximum  element  node  number  difference) 

X  (average  number  of  degree  of  freedom  per 
node) 

NEQ  =  total  number  of  degrees  of  freedom  in  the 
model 

=  (6)  x  (total  number  of  nodes)  -  [number  of 
fixed  (deleted)  degrees  of  freedom] 

NEQB  =  number  of  equations  per  block  of  storage 

*  MTOT/  MBAND/  2  (for  large  systems) 

If  NEQB  is  less  than  NEQ,  the  model  is  too  large  for  the  DETERMINANT 
SEARCH  algorithm,  and  the  SUBSPACE  ITERATION  procedure  will  be  used. 

If  the  SUBSPACE  ITERATION  algorithm  is  used  the  user  may  request  that 
the  STURM  SEQUENCE  check  be  oer formed.  By  experience  the  algorithm  has 
always  produced  the  lowest  NF  eigenvalues,  but  there  is  no  formal  mathe¬ 
matical  proof  that  the  calculated  NF  eigenvalues  will  always  be  the  lowest 
ones.  The  STURM  SEQUENCE  check  can  be  used  to  verify  that  the  lowest  NF 
eigenvalues  have  been  obtained.  It  should  be  noted  that  the  computational 
effort  expended  in  performing  the  STURM  SEQUENCE  check  is  not  trivial.  A 
factorization  of  the  complete  system  matrix  is  performed  at  a  shift  just  to 
the  right  of  the  NFth  eigenvalue. 

If  during  the  SUBSPACE  ITERATIONS  the  NFth  eigenvalue  fails  to  converge 
to  a  tolerance  of  "RTOL"  (normally  1.0E-5,  or  5  significant  figures)  within 
"NITEM"  (normally  "16")  iterations,  then  the  STURM  SEQUENCE  flag  ("IFSS") 
is  ignored. 
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(3)  The  maximum  number  of  iterations  to  reach  convergence  ("NITEM") 
applies  only  to  the  SUBSPACE  ITERATIONS  algorithm.  If  columns  11-15  are 
left  blank,  a  default  valued  of  "16"  for  NITEM  is  assumed. 

(4)  The  convergence  tolerance  ("RTOL")  is  applicable  only  if  the 
SUBSPACE  ITERATION  algorithm  is  used.  This  tolerance  test  applies  to  the 
NFth  eigenvalue,  and  all  eigenvalues  lower  than  the  NFth  one  will  be  more 
accurate  than  RTOL.  The  lowest  mode  is  found  most  accurately  with 
precision  decreasing  with  increasing  mode  number  until  the  highest 
requested  mode  ("NF")  is  accurate  to  a  tolerance  of  RTOL.  Iteration  is 
terminated  after  cycle  number  (k+1)  if  the  NFth  eigenvalue  (X,  say) 
satisfies  the  inequality: 

l|X(k+l)  -  X(k) | /X(k) ]  <  RTOL 

If  the  determinant  search  algorithm  is  used,  the  eigenpairs  are 
obtained  to  a  high  precision,  which  is  indicated  by  the  "physical  error 
bounds" 

-  HrJlj  /|  |K4±|  l2 

where 

ri  =  (K  -  Wl2M)  #i, 

2 

and  (w  i  )  are  the  i'th  eigenvalue  and  eigenvector  obtained  in  the 
solution. 

(5)  The  cut-off  frequency  ("COFQ")  is  used  by  both  eigenvalue 
algorithms  to  terminate  computations  if  all  eigenvalues  below  the  specified 
frequency  have  been  found. 

The  DETERMINANT  SEARCH  algorithm  computes  eigenvalues  in  order  from 
"l"  to  "NF".  If  the  Nth  eigenvalue  (IS  NF)  has  a  frequency  greater  than 
"COFQ",  the  remaining  (NF-N)  eigenvalues  are  not  computed. 

The  SUBSPACE  ITERATION  algorithm  terminates  calculation  when  the  Nth 
eigenvalue  is  accurate  (i.e.,  does  not  change  with  iteration)  to  a  tolerance 
of  RTOL.  As  before,  the  Nth  eigenvalue  is  the  nearest  eigenvalue  higher 
than  COFQ.  If  the  SUBSPACE  ITERATION  solution  determines  N  eigenvalues 
less  than  COFQ  (where,  N<NF),  the  STURM  SEQUENCE  check  (if  requested)  is 
performed  using  the  Nth  (rather  than  the  NFth)  eigenvalue  as  a  shift. 

Only  those  modes  whose  frequencies  are  less  than  COFQ  will  be  used  in 
the  RESPONSE  HISTORY  or  RESPONSE  SPECTRUM  analyses. 

(6)  The  starting  iteration  vectois,  together  with  control  informa¬ 
tion,  must  be  written  onto  TAPE10  before  the  program  execution  is  started. 
The  section  Restart  Instruction  for  Obtaining  Additional  Eigenvalues  and 
Eigenvectors  describes  the  creation  of  TAPE10  and  gives  the  required 
control  lines. 


88 


(7)  The  program  does  not  calculate  rigid  body  modes,  i.e,  the  system 
must  be  restrained  so  that  no  rigid  body  modes  are  present.  In  exact 
arithmetic  the  element  d  of  the  matrix  D  in  the  triangular  factorization 
of  the  stiffness  matrix  f?.e.  K  =  LDL1)  is  zero  if  a  rigid  body  mode  is 
present.  In  computer  arithmetic  the  element  d  is  small  when  compared 
with  the  other  elements  of  the  matrix  D.  If  tR?s  condition  occurs  the 
program  stops  with  an  ill-condition  matrix  message. 

(8)  If  many  "artificially"  stiff  boundary  elements  are  used,  the 
average  of  the  elements  of  D  will  be  artificially  large.  Consequently,  d 
may  be  small  in  comparison,  and  although  no  rigid  body  modes  may  be  present, 
the  program  will  stop.  In  a  dynamic  analysis  it  is  recommended  not  to  use 
very  stiff  boundary  elements. 
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RESPONSE  BISTORT  ANALYSIS  (NDYN  equals  2  or  NDYN  equals  4) 

The  NDYN  equals  2  option  uses  the  ("NF")  mode  shapes  and  frequencies 
computed  in  the  proceeding  section  to  perform  a  mode  superposition  solution 
for  forced  response.  The  NDYN  equals  4  option  initiates  a  direct 
step-by-step  integration  of  the  coupled  system  equation,  i.e.,  no 
eigenvalue  solution  has  been  performed  and  no  transformation  to  the 
eigenvector  basis  is  carried  out.  The  data  input  is  identical  to  the  case 
NDYN  equals  2  except  for  the  definition  of  damping.  Dynamic  response  can 
be  produced  by  two  general  types  of  forcing  function: 

(1)  ground  acceleration  input  in  any  (or  all)  or  the  three  (3)  global 
(X,Y,Z)  directions; 

and/or 

(2)  time  varying  loads  (forces/moments)  applied  in  any  (or  all)  nodal 
degrees  of  freedom  (except  -  "slave"  degrees  of  freedom) 

Time  dependent  forcing  functions  (whether  loads  or  ground  acceleration 
components)  are  described  in  two  steps.  First,  a  number  (1  or  more  are 
possible)  of  nondimens ional  time  functions  are  specified  tabularly  by  a  set 
of  discrete  points:  [f(t  ),t.],  where  i  =  l,2,...,k.  Each  different  time 
function  may  have  a  different  number  of  definition  points  (k).  A  particular 
forcing  foundation  applied  at  some  point  on  the  structure  is  then  defined 
by  a  scalar  multiplier  ("IS",  say)  and  reference  to  one  of  the  input  time 
functions  ("f(t),",  say).  The  actual  force  (or  acceleration)  at  any  time 
(’V'.say)  equals  3  X  f(t);  f (t )  is  found  by  linear  interpolation  between 
two  of  the  input  time  points  ft^ti+lj,  where  t^  £  t  it^+^. 

Assuming  that  the  solution  begins  at  time  zero  (0),  an  independent 
arrival  time  (tg,  where  t  k  0)  may  be  assigned  to  each  forcing  function. 

The  forcing  function  is  not  applied  to  the  system  until  the  solution  time 
("t",  say)  equals  the  arrival  time,  t  .  Interpolation  for  function  values 
is  based  on  relative  time  within  the  function  table;  i.e.,  g(i)  *  f(t  -  t  ). 

The  structure  is  assumed  to  be  at  rest  at  time  zero;  i.e.,  zero 
initial  displacements  and  velocities  are  assumed  at  time  of  solution  start. 

The  following  data  are  required  for  a  Forced  Dynamic  Response 
History  Analysis: 


A.  Response  History  Control  Line  (5I5,2F10.0) 


Notes 

Columns 

Variables 

Entry 

(1) 

1  -  5 

NFN 

Number 
6E.  1 

of  different  time  functions; 

(2) 

6-10 

NGM 

Ground 
EQ.  0; 
EQ.  1; 

motion  indicator; 
no  ground  motion  is  input 
read  ground  motion  control 
line  (Section  C  below) 
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Notes 


Columns 


Variables 


Entry 


(3) 

11 

-  15 

NAT 

Number  of  different  arrival  times 

for  the  forcing  functions; 

EQ.  0;  all  arrival  times  are  zero 

(A) 

16 

-  20 

NT 

Total  number  of  solution  time 
steps;  .GE.  1 

(5) 

21 

-  25 

NOT 

Output  print  interval  for 
stresses,  displacements,  etc. 

1  S  NOT  £  NT 

(4) 

26 

-  35 

DT 

Solution  time  step.  At; 

GT.  0 

(6) 

36 

-  45 

DAMP 

Damping  factor  to  be  applied  to  all 
NF  modes  (fraction  of  critical) 

GE.  0 

In 

case  of 

direct  integration  (NDYN  equals  4)  use 

(6) 

36 

-  45 

ALPHA 

Damping  factor  a 

(7) 

46 

-  55 

BETA 

Damping  factor  3 

NOTES: 

(1) 

At 

least  one  time  function  must  be  input. 

(2) 

If 

no  ground  acceleration  acts  on  the  structure,  set  "NGfT  to 

zero  and  skip  the  ground  motion  control  line,  (Section  C  below).  Both 
ground  acceleration  and  nodal  force  input  are  allowed. 

(3)  If  no  arrival  time  values  are  input,  all  forcing  functions  begin 
acting  on  the  structure  at  time  zero.  The  same  arrival  time  value  may  be 
reference  by  different  forcing  function.  "NAT"  determines  the  number  of 
nonzero  entries  that  the  program  expects  to  read  in  the  arrival  time  line, 
(Section  D  below). 

(4)  The  program  performs  a  step-by-step  integration  of  the  equations 
of  motion  using  a  scheme  which  is  unconditionally  stable  with  respect  to 
time  step  size.  At.  In  case  NDYN  equals  2  the  modal  uncoupled  equations  of 
motion  are  integrated.  In  case  NDYN  equals  4  the  coupled  system  equations 
are  integrated.  If  "T"  is  the  period  of  the  highest  numbered  mode  (normally 
the  NFth  mode)  that  is  to  be  included  in  the  response  calculation.  At  should 
be  chosen  such  that  At/T  <0.1.  A  larger  time  step  (i.e..  At  >  0.1T)  will 
not  cause  failure  (instability),  but  participation  of  the  higher  modes  is 
"filtered"  from  the  predicted  response.  In  general,  with  increasing  time 
step  size  the  solution  is  capable  of  capturing  less  of  the  higher  frequency 
participation. 

(5)  The  program  computes  system  displacement  at  every  solution  time 
step,  but  printing  of  displacements  and  recovery  of  element  stresses  is 
only  performed  at  solution  step  intervals  of  "NOT”.  NOT  must  be  at  least 
"l"  and  is  normally  selected  In  the  range  of  10  to  100. 
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(6)  The  damping  factor  ("DAMP")  is  applied  to  all  NF  modes.  The 
admissible  range  for  DAMP  is  between  0.0  (no  damping)  and  1.0  (100  percent 
of  critical  viscous  damping). 


(7)  For  direct  integration  (NDYN  equals  4)  the  damping  matrix  used  is 
C  ■  oM  +  BK,  where  o  and  B  are  defined  in  columns  36  -  55. 


B.  Time-Varying  Load  Lines  (4I5,F10.0) 


Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

NP 

Nodal  point  number  where  the  load 
component  (force  or  moment)  is  applied; 
1  S  NP  S  NUMNP 

EQ.  0  last  line  only 

(2) 

10 

IC 

Degree  of  freedom  number; 

1  S  IC  S  6 

( 6X=1 , 6Y=2 , 6Z=3 , *X=4 , *Y5 , *Z=6 ) 

(3) 

11  -  15 

IFN 

Time  function  number; 

1  S  IFN  S  NTFN 

(4) 

16  -  20 

IAT 

Arrival  time  number; 

EQ.  0;  load  applied  at  solution  start 
6E.  1;  nonzero  arrival  time 

(5) 

.1 

NOTES: 

21  -  30 

P 

Scalar  multiplier  for  the  time  function 
EQ.  0;  no  load  applied 

(1)  One  line  is  required  for  each  nodal  degree  of  freedom  having 
applied  time  varying  loads.  Lines  must  be  input  in  ascending  node  point 
order.  This  sequence  of  lines  must  be  terminated  with  a  blank  line.  A 
blank  line  must  be  supplied  even  if  no  loads  are  applied  to  the  systems. 

(2)  The  same  node  may  have  more  than  one  degree  of  freedom  loaded; 
arrange  degrees  of  freedom  reference  ("IC")  in  ascending  sequence  at  any 
given  node. 

(3)  A  nonzero  time  function  number  ("IFN")  must  be  given  for  each 
forcing  function.  IFN  must  be  between  1  and  NFN.  The  time  functions  are 
input  tabularly  in  the  time  function  definition  section,  below.  Function 
values  at  times  between  input  time  point  are  computed  with  linear  inter¬ 
polation. 

(4)  If  "IAT"  is  zero  (or  blank),  the  forcing  function  is  assumed  to 
act  on  the  system  beginning  at  time  zero.  If  IAT  is  input  as  a  positive 
integer  between  1  and  NAT,  the  IATth  arrival  time  (defined  in  the  arrival 
time  section,  below)  is  used  to  delay  the  application  of  the  forcing 
function;  i.e.,  the  forcing  function  beings  acting  on  the  structure  when 
the  solution  reaches  the  IATth  arrival  time  value. 

(5)  The  actual  magnitude  of  force  (or  moment)  acting  on  the  model  at 
time,  t,  equals  the  product:  ("P")  x  (value  of  function  number  "IFN"  at 
time,  t). 
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C.  Ground  Motion  Control  Line  (615) 


Notes 

Columns 

Variables 

Entrv 

(1) 

1-5 

NFNX 

Time  function  number  describing  the 
ground  acceleration  in  the  X-direction 

6-10 

NFNY 

Time  function  number  describing  the 
ground  acceleration  in  the  Y-direction 

11  -  15 

NFNZ 

Time  function  number  describing  the 
ground  acceleration  in  the  Z-direction 

(2) 

16  -  20 

NATX 

Arrival  time  number,  X-direction 

21  -  25 

NATY 

Arrival  time  number,  Y-direction 

26  -  30 

NATZ 

Arrival  number,  Z-direction 

NOTES: 

(1)  This  line  must  be  input  only  in  the  ground  motion  indicator 
("NGM")  was  set  equal  to  one  (1)  on  the  Response  History  Control  Line, 
above.  A  zero  time  function  number  indicates  that  no  ground  motion  is 
applied  for  that  particular  direction. 

(2)  Zero  arrival  time  references  mean  that  the  ground  acceleration 
(if  applied)  begins  acting  on  the  structure  at  time  zero).  Nonzero 
references  must  be  integers  in  the  range  1  to  NAT. 


D. 

Arrival  Time  Lines 

a.  Line 

One  (8F10.0) 

Notes 

Columns 

Variables 

Entrv 

(1) 

1-10 

AT(1) 

Arrival 

time 

number  1  4 

11  -  20 

AT(2) 

Arrival 

time 

number  2  w. 

71  -  80 

AT(8) 

Arrival 

time 

number  8 

b.  Line 

TWo  (8F10.0)  - 

(required 

if  NAT.GT.  8) 

Notes 

Columns 

Variables 

Entrv 

1  -10 

AT(9) 

Arrival 

time  number  9 

etc. 

etc. 

NOTES: 

(1) 

The  entry  ("NAT")  given  in  column 

11-15  on  the  Response  History 

Control  Mne  above  specifies  the  total  number  of  arrival  time  entries  to  be 
read  in  this  section.  Input  as  many  lines  as  are  required  to  define  "NAT" 
different  arrival  times,  eight  entries  per  line.  If  no  arrival  times  were 
requested  (NAT  equals  0)  supply  one  blank  line  in  this  section. 
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E.  Time  Function  Definition  Lines 

Supply  one  set  (line  1  and  card(s)  2)  of  input  for  each  of  the  " NFN " 
time  functions  requested  in  column  1-5  of  the  Response  History  Control  Line 
above.  At  least  one  set  of  time  function  lines  is  expected  in  this  section. 
The  line  sets  are  input  in  ascending  function  number  order. 

Line  1  (I5,F10.0, 12A5) 


Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

NLP 

Number  of  function  definition  points; 
GE.  2 

(2) 

6  -  15 

SFTR 

Scale  factor  to  be  applied  to  f(t) 
values; 

EQ.  0;  default  set  to  "l.O" 

16  -  75 

HED(12) 

Label  information  (to  be  printed 
with  output)  describing  this 
function  table 

NOTES: 

(1)  At  least  two  points  (i.e.,  2  pairs:  f(t  ),t^)  must  be  specified 
for  each  time  function.  Less  than  two  points  would  preclude  linear 
interpolation  in  the  table  for  f(t). 

(2)  The  scale  factor  "SFTR"  is  used  to  multiply  function  values  only; 
i.e.,  input  time  values  are  not  changed.  If  the  scale  factor  is  omitted, 
SFTR  is  re-set  by  the  program  to  ”1.0"  thereby  leaving  input  function 
values  unchanged. 


Line' 2  (12F6.0) 

Notes 

Columns 

Variables 

Entry 

(1) 

1-6 

T(l) 

Time  values  at  point  1,  t^ 

7-12 

F(l) 

Function  value  at  point  1, 

f(t.) 

13  -  18 

T(2) 

Time  value  at  point  2,  t_ 

19  -  24 

F(2) 

Function  value  at  point  2, 

f(t.) 

etc. 

etc. 

4 

NOTES: 

(1)  Input  as  many  line(s)  2  as  are  required  to  define  the  time 
function.  The  "NLP"  pairs  must  be  input  in  order  of  ascending  time  value. 
Time  at  point  one  must  be  zero,  an£  care  must  be  taken  to  ensure  that  the 
highest  (last)  input  time  values  (NLP  )  is  at  least  equal  to  the  value  of 
time  at  the  end  of  solution;  i.e.,  the  time  span  for  all  functions  must 
cover  the  solution  time  period  otherwise  the  interpolation  for  function 
values  will  fail.  For  the  case  of  nonzero  arrival  times  associated  with  a 


particular  function,  the  shortest  grrlval  time  reference  ("t.",  say) 
plus  (+)  the  last  function  time  £"  NLP")  must  at  least  equalAthe  time  at 
the  end  of  the  solution  period  (  END),  say);  i.e.,  t^  +  t^p  2  t^^. 
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F.  Output  Definition  Lines 


To  minimize  the  amount  of  the  output  which  would  be  produced  by  the 
program  if  all  displacement,  stresses,  etc.  were  printed,  output  requests 
for  specific  components  must  be  given  in  this  section.  Time  histories  for 
selected  components  appear  in  tables;  the  solution  step  output  printing 
interval  is  specified  as  "NOT"  which  is  given  in  column  21-25  of  the 
Response  History  Control  Line  above. 


Displacement  Output  Requests.  Control  Line  (215) 


Notes 

Columns 

Variables 

Entrv 

(1) 

1-5 

KKK 

Output  type  indicator; 

EQ.  1;  print  histories  and  maxima 

EQ.  2;  printer  plot  histories  and 
recovery  of  maxima 

EQ.  3;  recovery  maxima  only 

(2) 

6-10 

ISP 

Printer  plot  spacing  indicator 

NOTES: 

(1)  The  type  of  output  to  be  produced  by  the  program  applies  to  all 
displacement  requests.  KKK  equals  0  is  not  permitted. 

(2)  "ISP"  controls  the  vertical  (down  the  page)  spacing  for  printer 
plots.  Output  points  are  printed  on  every  (ISP+l)th  line.  The  horizontal 
(across  the  page)  width  of  printer  plots  is  a  constant  ten  (10>  inches  (100 
print  positions).  ISP  is  used  only  if  KKK  equals  2. 


Node  Displacement  Request  Lines  (715) 

Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

NP 

Node  number 

GE.  1  and  LE.NUMNP 

EQ.  0  last  line  only 

(2) 

6-10 

IC(1) 

Displacement  component. 

request  1 

11  -  15 

IC(2) 

Displacement  component. 

request  2 

16  -  20 

IC(3) 

Displacement  component. 

request  3 

21  -  25 

IC(4) 

Displacement  component. 

request  4 

26  -  30 

IC(5) 

Displacement  component, 

request  5 

31  -  36 

IC(6) 

Displacement  component. 

request  6 

GE .  1  and  LE . 6 

EQ.  0  terminates  requests  for  the  node 

NOTES: 

(1)  Only  those  nodes  at  which  output  is  to  be  produced  (or  at  which 
maxima  are  to  be  determined)  are  entered  in  this  section.  Lines  must  be 
input  in  ascending  node  number  order.  Node  numbers  may  not  be  repeated. 
This  section  must  be  terminated  with  a  blank  line. 


(2)  Displacement  component  request  ("IC")  range  from  1  to  6,  where 
1=6X,  2=6 Y ,  3=6 Z ,  4=#X,  5=fY,  6=#Z.  The  first  zero  (or  blank)  encountered 
while  reading  IC(1),  IC(2), . . . ,IC(6)  terminates  information  for  the  line. 
Displacement  components  at  a  node  may  be  requested  in  any  order.  For 
example,  suppose  that  6Y,4X  and  iZ  are  to  be  output  at  node  34;  the  line 
could  be  written  as  /34,2,4,6,0/,  or  /34,6,4,2,0/,  etc.  but  only  four  (4) 
fields  would  have  nonzero  entries. 


Element  Stress  Component  Output  Requests.  Control  Line  (215) 


Notes 

Columns 

Variables 

Entry 

(1) 

1  -  5 

KKK 

Output  type  indicator; 

EQ.  1;  print  histories  and  maxima 
EQ.  2;  printer  plot  of  histories 
and  recovery  of  maxima 

EQ.  3;  recover  maxima  only 

(2) 

6-10 

ISP 

Plot  spacing  indicator 

NOTES: 

(1)  The  type  of  output  to  be  produced  by  the  program  applies  to  all 
displacement  requests.  KKK  equals  0  is  not  permitted. 

(2)  "ISP"  controls  the  vertical  (down  the  page)  spacing  for  printer 
plots.  Output  points  are  printed  on  every  (ISP+l)th  line.  The  horizontal 
(across  the  page)  width  of  printer  plots  is  a  constant  ten  (10)  inches  (100 
print  positions).  ISP  is  used  only  if  KKK  equals  2. 

Element  Stress  Component  Request  Lines  (1315) 

Requests  are  grouped  by  element  type;  "NELTYP"  groups  must  be  input. 

A  group  consists  of  a  series  of  element  stress  component  request  lines 
terminated  by  a  blank  line.  Element  number  references  within  an  element 
type  (e.g.  TRUSS)  grouping  must  be  in  ascending  order.  Element  number 
references  may  be  omitted  but  not  repeated.  The  program  processes  element 
groups  in  the  same  order  originally  input  in  the  Element  Data  Section.  If 
no  output  is  to  be  produced  for  an  element  type,  then  input  one  blank  line 
for  its  group. 


Notes 

Columns 

Variables 

Entry 

(1) 

1  -  5 

NEL 

Element  number 
GE.  1 

EQ.  0;  last  line 

in  the 

group  only 

(2) 

6-10 

IS(1) 

Stress  component 
request  1 

number 

for  output. 

11  -  15 

IS(2) 

Stress  component 

number 

for  output. 

request  2 

61  -  65 

IS( 12) 

Stress  component 
request  12 

number 

for  output, 

END  OF  DATA  SET  FOR  NDYN  =  2  OR  NDYN  =  4 
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NOTES: 

(1)  Terminate  each  different  element  output  group  (type)  with  a  blank 
line.  Elements  within  a  group  must  be  in  element  number  order  (ascending); 
element  number  repetitions  are  not  permitted. 

(2)  The  first  zero  (or  blank)  request  encountered  while  reading 
IS(1),  IS(2),...,  IS(12)  terminated  information  for  the  line.  No  more 
than  twelve  (12)  different  components  may  be  output  for  any  one  of  the 
elements.  Appendix  A  lists  the  stress  component  numbers  and  corresponding 
descriptions  for  the  various  elements  types.  Some  element  types  (TRUSS, 
for  example)  have  fewer  than  12  components  defined;  only  the  stress 
component  number  listed  in  Appendix  A  are  acceptable  references. 
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RESPONSE  SPECTRUM  ANALYSIS  (NDYN  equals  3) 

This  option  combines  all  (NF)  mode  shapes  and  frequencies  computed 
during  the  eigenvalue  solution  to  calculate  root  mean  square  (R.M.S.) 
stresses/deflections  due  to  an  input  displacement  (or  acceleration) 
spectrum.  The  input  spectrum  is  applied  in  varying  proportions  in  the 
global  X,Y,Z  directions.  For  the  case  of  a  nonzero  cut-off  frequency 
"COFQ"  only  those  modes  whose  frequencies  are  less  than  COFQ  will  be 
combined  in  the  R.M.S  analysis. 

A.  Control  Line  (3F10.0.I5) 


Notes 

Columns 

Variables 

Entry 

(1) 

1-10 

FX 

Factor  for  X-direction  input 

11  -  20 

FY 

Factor  for  Y-direction  input 

21  -  30 

FZ 

Factor  for  Z-direction  input 

EQ.  0;  not  acting 

(2) 

31  -  35 

1ST 

Input  spectrum  type; 

EQ.  0;  displacement  vs.  period 
EQ.  1;  acceleration  vs.  period 


NOTES: 

(1)  All  three  direction  factors  may  be  nonzero  in  which  case  the 
entries  represent  the  X,Y,Z  components  of  the  input  direction  vector. 

(2)  1ST  defines  the  type  of  spectrum  table  to  be  input  immediately 

following.  The  spectral  displacements  ("S.")  and  accelerations  (S  )  are 
assumed  to  be  related  as  follows:  S  =  (4f*f*)  (S.).  a 


B.  Spectrum  Lines 

a.  Heading  Line  (12A6) 

Notes  Columns  Variables  Entry 

1-72  HED(12)  Heading  information  used  to  label 

the  spectrum  table 

b.  Control  Line  (I5,F10.0) 


Notes 

Columns 

Variable 

1  -  5 

NPTS 

6-15 

SFTR 

Entry 

Number  of  definition  points  in  the 
spectrum  table; 

6E.  2 

Scale  factor  used  to  adjust  the 
displacement  (or  acceleration) 
ordinates  in  the  spectrum  table 
EQ.  1.0;  no  adjustment 
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c.  Specti 

m  Data  (2F10.0) 

Notes 

Columns 

Variables 

Entry 

Cl) 

1-10 

T 

Period  (reciprocal  of  frequency) 

(2) 

11  -  20 

S 

Value  of  displacement  (or  acceleration 
if  1ST  equals  1) 

END  OF 

DATA  SET  FOR 

NDYN  =  3 

NOTES: 

(1)  Input  one  definition  point  per  line;  "NPTS"  lines  are  required  in 
this  section.  Lines  oust  be  arranged  in  ascending  value  of  period. 

(2)  S  is  interpreted  to  be  a  displacement  quantity  if  "1ST"  was  input 
as  zero.  For  1ST  equals  1,  "S"  is  an  acceleration  value. 
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EXECUTION  INSTRUCTIONS 


MS-DOS  COMPUTERS 

Execution  of  SAP  IV  on  a  microcomputer  can  only  be  done  on  a  fixed 
disk  system  with  a  math  coprocessor  chip.  The  primary  memory  requirement 
is  640K,  and  the  fixed  disk  must  have  at  least  10  Megabytes  available  dur¬ 
ing  a  large  run.  Some  smaller  runs  can  have  smaller  disk  space  require¬ 
ments.  The  program  reads  data  from  the  "SAP4.DAT"  file  and  if  the  run  is  a 
restart,  the  appropriate  restart  files  are  used. 

The  master  diskette  contains  several  files.  Tie  three  SAP  IV  diskettes 
are  360K  double  sided,  double  density  disks  so  most  machines  can  read  them. 
There  are  three  parts  of  the  SAP4  executable  program:  SAP4.EXE,  SAP4C0M.0VL, 
and  SAP4DYN.OVL.  The  RUNSAP4.BAT  file  is  set  up  so  the  user  can  name  the 
data  file  as  desired  and  the  output  file  will  be  named  with  the  same  name, 
however  the  input  data  file  must  use  the  extension  .DAT  and  the  output  file 
will  be  assigned  the  extension  .OUT.  The  RUNSAP4.BAT  file  also  erases  all 
of  the  TAPE  files  that  are  required  for  restarts  and  post -processing.  This 
convention  can  be  changed  to  suit  user  restart  needs. 

The  current  version  of  SAP  IV,  supplied  on  the  distribution  diskette, 
has  been  set  up  so  it  will  solve  medium  sized  problems.  The  problem  size 
is  determined  by  the  number  of  nodes,  resulting  in  the  following  size 
limitations: 

Maximum  number  of  nodes  is  1,500 

Maximum  storage  allocation  (MTOT)  is  15,000 

Executable  program  size  is  400K  bytes 

Number  of  disk  files  simultaneously  open  is  13 

Recommended  number  of  file  transfer  buffers  is  13  (7K) 

There  are  a  number  of  options  available  to  modify  the  program  so  it  can 
address  larger  problems.  NCEL  can  assist  users  who  have  requirements  that 
exceed  the  current  capacity  of  the  program.  Options,  such  as  extended 
memory  and  random-access  memory  (RAM)  disk  are  available. 


OTHER  COMPUTERS 

This  program  can  be  used  on  other  computers.  Contact  NCEL  Code  L51 
for  instructions. 
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DATA  FILES  CREATED  DURING  EXECUTION 

SAP  IV  creates  eleven  data  files  during  execution.  These  files  are 
used  for  various  reasons  depending  upon  the  type  of  solution  being  conducted. 
Their  size  is  also  problem  dependent.  Most  problems  will  run  if  10  mega¬ 
bytes  of  disk  space  are  available.  The  files  are  unformatted,  meaning  they 
can  not  be  examined  unless  a  special  program  is  used.  However,  the 
TAPE11,  TAPE 12,  SAP4.DAT,  and  the  SAP4.0UT  files  are  formatted  so  they  can 
be  examined.  The  table  below  provides  some  information  about  each  file. 

File  Contents 

TAPE1  Restart:  Element  description  data 

TAPE2  Restart:  Stiffness  matrix,  displacement 

TAPE3  Scratch:  Mass  matrix 

TAPE4  Scratch:  Nodal  displacement 

SAP4.DAT  Initial  input  data 

SAP4.0UT  Printer  output  data 

TAPE7  Restart:  Subspace  iteration 

TAPE8  Restart:  Dynamic  analysis 

Postprocessing:  Geometry  and  element  data 
TAPE9  Restart:  Mass  matrix 

TAPE10  Postprocessing:  Element  stress  results  4 
TAPE11  Postprocessing:  Time  history  results 
TAPE12  Postprocessing:  Problem  size  data 
TAPE13  Scratch:  Stress -displacement  transformation  data 
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DYNAMIC  ANALYSIS  RESTART  INSTRUCTIONS 


The  purpose  of  this  section  is  to  describe  the  procedure  required  for 
program  restart  following  an  eigenvalue/ eigenvector  analysis.  The  restart 
option  has  been  include  in  the  program  in  order  to  make  a  repeated  forced 
response  or  spectrum  analysis  possible  without  solving  each  time  for  the 
required  eigensystem.  For  medium-to- large  size  models,  eigenvalue  solution 
is  computation  intense  when  compared  to  the  forced  response  calculations; 
hence,  excessive  computation  time  may  be  incurred  if  the  entire  job  has  to 
be  rerun  due  to  improper  specification  of  forcing  functions  or  input 
spectra,  inadequate  requests,  etc.  For  small  models  (less  than  100  nodes) 
the  extra  effort  required  for  restart  is  normally  not  justified. 

A  complete  dynamic  analysis  utilizing  the  restart  feature  requires 
that  the  analysis  be  conducted  in  two  steps: 

Step  1:  Eigenvalue  extraction  solution  only,  after  which  program  files 

TAPE1,  TAPE2,  TAPE7,  TAPES,  and  TAPE9  are  saved. 

Step  2:  Response  Analysis  (NDYN  equals  -2)  or  a  Response  Spectrum  Analysis 

(NDYN  equals  -3)  using  previously  saved  files  TAPE1,  TAPE2,  TAPE7, 
TAPE 8,  and  TAPE9. 

For  a  given  model,  the  first  step  creating  the  restart  tapes  is 
executed  once.  The  restart  tapes  then  contsin  all  the  initial  information 
required  by  the  program  at  the  beginning  of  a  forced  response  analysis. 

Step  2.  More  than  one  step  2  may  be  executed  using  the  restart  tapes. 

Step  1  PSEUDO  Problem  Data  File: 

HEADING  LINE 
MASTER  CONTROL  LINE  with 
LL  ■  0 
NF  2  1 
NDYN  «  0 
MODEX  *=  0 
JOINT  DATA 
ELEMENT  DATA 
CONCENTRATED  MASS  DATA 
ELEMENT  LOAD  MULTIPLIERS 

DYNAMIC  ANALYSIS:  Mode  shapes  and  Frequencies 
Blank  line 
Blank  line 
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Step  2  PSEUDO  Problen  Data  File: 

HEADING  LINE 

MASTER  CONTROL  LINE  with 
LL  «  0 
NF  *  1 

NDYN  =  -2  or  -3 
MODEX  =  0 
DYNAMIC  ANALYSIS 

Dynaaic  Response  Analysis  (NDYN  equals  -2) 
or 

Response  Spectrum  Analysis  (NDYN  equals  -3) 

Blank  line 
Blank  line 

Normally,  the  number  of  frequencies  ("NF")  entered  on  the  MASTER  CONTROL 
LINE  for  a  restart  case  has  the  same  value  as  was  specified  earlier  when 
the  eigenvalue  problem  was  solved  in  Step  1.  If  a  value  for  the  cut-off 
frequency  ("COFQ")  was  entered  on  the  "Mode  Shapes  and  Frequencies"  control 
line  (Step  I)  and  the  program  extracted  fewer  than  "NF"  frequencies 
(eigenvalues),  then  only  the  actual  number  of  eigenvalues  computed  by  the 
program  in  (Step  1)  is  specified  for  "NF"  in  this  restart  run. 
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RESTART  INSTRUCTIONS  FOR  OBTAINING  ADDITIONAL 
EIGENVALUES  AND  EIGENVECTORS 


In  the  dynamic  analysis  of  large-order  systems,  the  solution  of  the 
required  eigensystem  is  normally  the  most  computation  intensive  phase.  The 
option  described  in  this  section  demonstrates  how  it  is  possible  to  use  NFO 
previously  calculated  eigenvalues  and  vectors  when  the  solution  for  NF  2; 

NFO  eigenvalues  and  eigenvectors  is  required. 

Assume  that  in  Step  1,  the  solution  for  NFO  eigenvalues  and 
eigenvectors  was  performed.  At  the  end  of  this  step,  TAPE2  and  TAPE7  must 
have  been  saved.  Assuming  that  in  Step  2  the  solution  of  NF  additional 
eigenvalues  and  eigenvectors  is  required,  then  prior  to  the  execution  of 
this  job,  the  previously  saved  TAPE2  and  TAPE7  files  are  combined  into  one 
file,  TAPE10,  and  that  file  is  rewound. 
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PRE-  AND  POST-PROCESSING 


The  key  to  using  the  finite  element  method  is  computer  applications. 

A  program  like  SAP  IV  requires  a  large  amount  of  coordinate,  connectivity, 
element  property,  and  load  information  to  solve  common  structural  problems. 
The  results  of  the  analysis  are  often  very  voluminous  especially  for 
dynamic  analysis  problems.  This  section  describes  various  ways  the  user 
can  simplify  the  task  of  preparing  the  input  data  and  of  examining  the 
output  data. 

PRE-PROCESSING 

Pre-processing  involves  computer  applications  to  generate  the  modal 
coordinates,  modal  connectivity,  element  properties,  and  load  information. 

A  number  of  options  are  available  to  simplify  data  preparation: 

Manual  Preparation 

The  SAP  IV  program  uses  fixed  format  line  (line)  image  method  of 
input.  The  lines  of  data  are  collected  into  an  input  data  file.  This 
method  of  input  is  very  laborious  for  large  problems  and  therefore  is 
not  recommended.  However,  small  problems  can  often  be  quickly  posed 
using  this  method  with  the  aid  of  a  screen  editor  program.  The  screen 
editor  program  tab  features  for  right  and  decimal  justification  will 
simplify  preparing  the  input  file.  The  screen  editor  program  must  be 
capable  of  preparing  an  alpha-numeric  file  for  SAP  IV  with  the  10  numbers 
and  26  capital  characters  and  no  special  control  characters.  Experience 
has  shown  that  this  method  of  input  can  be  very  expedient  and  easy  to 
modify  for  the  smaller  class  problems.  Additionally,  the  SAP  IV  coordinate 
and  element  generation  features,  while  somewhat  limited,  can  be  used  to 
simplify  preparing  the  data. 

Automatic  Generators 

There  are  a  few  automatic  input  data  generators  that  have  been 
successfully  used  to  generate  SAP  IV  data. 

GEN3D  is  Naval  Civil  Engineering  Laboratory  program  that  generates 
the  node  and  element  connectivity  data.  This  program  is  technically  one 
of  the  most  advanced  mech  generators  that  is  available.  The  primary 
limitation  is  the  user  interface  is  less  not  desirable,  since  it  does 
not  support  immediate  user  interaction. 

Finite  Element  Modeling  (FEM)  software  is  the  best  way  available 
to  prepare  input  data  for  SAP  IV.  Such  a  package  encourages  the  user  to 
develop  a  three-dimensional  model  of  the  structure  and  then  provides 
comprehensive  commands  to  superimpose  a  finite  element  mesh  on  top  of 
the  model.  The  mesh  contains  the  coordinate  and  element  connectivity 
information.  The  package  can  also  be  used  to  generate  the  element 
properties  and  load  cases.  The  features  include  ways  to  check  the 
quality  of  the  mesh,  element,  and  load  information  before  the  analysis 
program  is  executed. 
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The  Naval  Civil  Engineering  Laboratory  has  developed  some  special 
purpose  generation  software  to  prepare  data  information  for  specific  types 
of  problems.  A  generator  for  overhead  electric  traveling  bridge  cranes  has 
been  prepared.  A  two-dimensional  drydock  model  generator  has  also  been 
prepared.  The  primary  strength  of  the  SAP  IV  program  is  having  the 
flexibility  of  developing  problem  specific  generation  software. 

Automatic  load  generation  has  been  accomplished.  Special  software 
has  been  developed  to  generate  hydrostatic  pressures  on  complex  curved 
surfaced.  Time  series  analysis  software  has  been  used  to  prepare  dynamic 
test  data  so  it  can  be  used  as  input  data.  The  dynamic  loads  include 
ocean  waves,  ship  motion,  wind,  and  blast  pressures. 

POST-PROCESSING 

Post -processing  involves  reducing  the  large  quantities  of  analysis 
results  into  meaningful  presentations  that  can  be  quickly  assimilated  and 
evaluated.  SAP  IV  prepares  a  number  files  (TAPE1,  TAPE8,  TAPE10,  TAPE11, 
TAPE12)  that  contain  the  model  undeformed  geometry,  the  displacement,  and 
the  various  stress,  force,  moment,  velocity,  and  accelerations  quantities. 
Programs  have  been  prepared  that  read  these  files  so  that  information  can 
be  further  processed  to  meet  specific  analysis  objectives. 

The  undeformed  geometry  can  be  plotted  using  the  SAPLOT  program. 

This  is  helpful  to  verify  the  correctness  of  the  coordinates  and  the 
element  connectivity. 

The  deformed  geometry  can  be  plotted  using  the  SAPPLOTD  program. 

This  program  will  plot  the  deformed  structure  with  or  without  the  unde¬ 
formed  structure  in  the  background. 

Various  programs  have  been  written  by  the  Naval  Civil  Engineering 
Laboratory  to  process  the  force,  moment,  and  stress  results.  A  beam 
stress  program  has  been  used  to  evaluate  the  beam  element  results  with 
respect  to  bending  and  shear  stress.  Moment  and  shear  diagrams  have 
also  been  produced.  The  two  and  three  dimensional  element  results  have 
been  studied  using  MOVIE. BYU  which  prepares  contour  and  continuous  tone 
displays  of  requested  information.  An  Interface  has  been  prepared  to 
aid  studying  time  history  data  using  time  series  analysis  techniques. 
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SAP  IV  -  VERSION  1.0 


FEEDBACK  REPORT 


The  Naval  Civil  Engineering  Laboratory  is  fully  dedicated  to  support¬ 
ing  GEMS  users.  A  primary  requirement  for  this  task  is  to  establish  a 
priority  listing  of  user  requirement.  It  would  be  of  great  value  to  the 
development  of  new  software  if  you,  the  user,  would  complete  the  feedback 
questions  below.  Since  each  individual  user  may  have  specific  requirement, 
please  reproduce  this  page  as  many  times  as  necessary. 


Please  circle  the  number  that  best  applies  in  questions  1  through 
4,  complete  the  other  questions,  fold  at  tic  marks,  and  mail  to  NCEL 

with  franked  label  on  reverse  side  or  to  address  at  bottom  of  page. 

1.  Was  the  software  beneficial  (productive)? 

No  benefit  0123456789  10  Very  beneficial 

2.  Was  it  easy  to  use  (user  friendly)? 

Difficult  0123456789  10  Very  easy 


3.  Does  this  software  make  decisions  more  reliable? 

No  0  1  2  3  4  5  6  7  8  9  10  Yes 

4.  Does  it  better  document  the  design? 

No  0123456789  10  Yes 

5.  Did  it  save  time? 

Yes _  No  Estimated  percent  saved _ 

6.  What  would  make  future  software  more  user  friendly? 


7.  What  further  support  would  you  like  to  have  on  the  GEMS  systems? 


8.  What  other  comments  or  remarks  would  you  like  to  add? 


Activity _ 

Telephone _ 

Mail  address  is: 

NAVFAC  GEMS  Support  Group 

Naval  Civil  Engineering  Laboratory 

Code  L54 

Port  Hueneme,  CA  93043-4328 
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APPENDIX  A 

OUTPUT  STRESS  COMPONENTS 


ELEMENT  TYPE 

1.  TRUSS 

MAXIMUM 
NUMBER  OF 
COMPONENTS 

(2) 

STRESS 

COMPONENT 

NUMBER 

(1) 

(2) 

OUTPUT 

SYMBOL 

(P/A  ) 
(P  ) 

DESCRIPT 

AXIAL  STRESS 
AXIAL  FORCE 

ION 

2. 

BEAM 

(12) 

(  1) 

(PI(I)  ) 

1 -Force 

at 

end 

I 

(  2) 

(V2(I)  ) 

2 -Shear 

at 

end 

I 

(  3) 

(V3(I)  ) 

3 -Shear 

at 

end 

I 

(  4) 

(TI(I)  ) 

1 -Torque 

at 

end 

I 

(  5) 

(M2(I)  ) 

2 -Moment 

at 

end 

I 

(  6) 

(M3(I)  ) 

3 -Moment 

at 

end 

I 

(  7) 

(P1(J)  ) 

1 -Force 

at 

end 

J 

(  8) 

(V2(J)  ) 

2 -Shear 

at 

end 

J 

(  9) 

(V3(J)  ) 

3 -Shear 

at 

end 

J 

(10) 

(TI( J)  ) 

1 -Torque 

at 

end 

J 

(ID 

(M2( J)  ) 

2 -Moment 

at 

end 

J 

(12) 

(M3(J)  ) 

3 -Moment 

at 

end 

J 

PLANE- 

(20) 

(  1) 

(11-SO 

) 

V-  Stress 

at 

point  0 

STRESS/ 

(  2) 

(22-SO 

) 

U-  Stress 

at 

point  0 

PLANE- 

(  3) 

(33-SO 

) 

T-  Stress 

at 

point  0 

STRAIN 

(  4) 

( 12-SO 

) 

UV-Stress 

at 

point  0 

(  5) 

(11-SI 

) 

V-  Stress 

at 

point  1 

(  6) 

(22-SI 

) 

U-  Stress 

at 

point  1 

(  7) 

(33-SI 

) 

T-  Stress 

at 

point  1 

(  8) 

( 12-SI 

) 

UV-Stress 

at 

point  1 

(  9) 

( 11-S2 

) 

V-  Stress 

at 

point  2 

(10) 

(22-S2 

) 

U-  Stress 

at 

point  2 

(ID 

(33-S2 

) 

T-  Stress 

at 

point  2 

(12) 

(12-S2 

) 

UV-Stress 

at 

point  2 

(13) 

( 11-S3 

) 

V-  Stress 

at 

point  3 

(14) 

(22-S3 

) 

U-  Stress 

at 

point  3 

(15) 

(33-S3 

) 

T-  Stress 

at 

point  3 

(16) 

(12-S3 

) 

UV-Stress 

at 

point  3 

(17) 

(V  -S4 

) 

V-  Stress 

at 

point  4 

(18) 

(U  -S4 

) 

U-  Stress 

at 

point  4 

(19) 

(T  -S4 

) 

T-  Stress 

at 

point  4 

(20) 

(UV-S4 

) 

UV-Stress 

at 

point  4 

A- 1 


(20) 


4.  TWO- 

DIMENSIONAL 

ELEMENTS 


5.  EIGHT  (12) 

NODE 
BRICK 


(  1) 

(11-SO 

) 

(  2) 

(22-SO 

) 

(  3) 

(33-SO 

) 

(  4) 

( 12-S0 

) 

(  5) 

(11-SI 

) 

(  6) 

(22-SI 

) 

(  7) 

( 33-SI 

) 

(  8) 

( 12-SI 

) 

(  9) 

( 11-S2 

) 

(10) 

(22-S2 

) 

(ID 

(33-S2 

) 

(12) 

( 12-S2 

) 

(13) 

( 11-S3 

) 

(14) 

(22-S3 

) 

(15) 

(33-S3 

) 

(16) 

( 12-S3 

) 

(17) 

(V  -S4 

) 

(18) 

(U  -S4 

) 

(19) 

(T  -S4 

) 

(20) 

(UV-S4 

) 

(  1) 

(XX-SL1) 

(  2) 

(YY-SL1) 

(  3) 

(ZZ-SL1) 

(  4) 

(XX-SL1) 

(  5) 

(YY-SL1) 

(  6) 

(ZZ-SL1) 

(  7) 

(XX-SL2) 

(  8) 

(YY-SL2) 

(  9) 

(ZZ-SL2) 

(10) 

(XX-SL2) 

(ID 

(YY-SL2) 

(12) 

(ZX-SL2) 

V-  Stress  at  point  0 
U-  Stress  at  point  0 
T-  Stress  at  point  0 
UV-Stress  at  point  0 

V-  Stress  at  point  1 
U-  Stress  at  point  1 
T-  Stress  at  point  1 
UV-Stress  at  point  1 

V-  Stress  at  point  2 
U-  Stress  at  point  2 
T-  Stress  at  point  2 
UV-Stress  at  point  2 
V-  Stress  at  point  3 
U-  Stress  at  point  3 
T-  Stress  at  point  3 
UV-Stress  at  point  3 

V-  Stress  at  point  4 
U-  Stress  at  point  4 
T-  Stress  at  point  4 
UV-Stress  at  point  4 


XX-Stress  at  location  1 
YY-Stress  at  location  1 
ZZ-Stress  at  location  1 
XY-Stress  at  location  1 
YZ-Stress  at  location  1 
ZX-Stress  at  location  1 

XX-Stress  at  location  2 
YY-Stress  at  location  2 
ZZ-Stress  at  location  2 
XY-Stress  at  location  2 
YZ-Stress  at  location  2 
ZX-Stress  at  location  2 


PLATE/ 

(  6) 

(  1) 

(XX-S/R) 

XX-Stress  Resultant 

SHELL 

(  2) 

(YY-S/R) 

YY-Stress  Resultant 

(  3) 

(XY-S/R) 

XY-Stress  Resultant 

(  4) 

(XX-M/R) 

XX-Moment  Resultant 

(  5) 

(YY-M/R) 

YY- Moment  Resultant 

(  6) 

(XY-M/R) 

XY-Moment  Resultant 

A-2 


7 .  BOUNDARY 


C  1) 
(  2) 


(BDRY-F)  BOUNDARY  Force 
(BDRY-M)  BOUNDARY  Moment 


(  2) 


8.  THICK  (  1) 

SHELL  (42)  (  2) 

AND  (3) 

3-DIAM  (  4) 

(  5) 

(  6) 


(  8) 
(  9) 
(10) 
(ID 
(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 
(AO) 

(41) 

(42) 


(SXX(O)) 

XX-Stress 

at 

centroid  (0) 

(SYY(O) ) 

XY-Stress 

at 

centroid  (0) 

(SZZ(O)) 

ZZ-Stress 

at 

centroid  (0) 

(SXY(O)) 

XY-Stress 

at 

centroid  (0) 

(SYZ(O) ) 

YZ-Stress 

at 

centroid  (0) 

(SZX(O)) 

ZX-Stress 

at 

centroid  (0) 

(SXX(l)) 

XX-Stress 

at 

center 

face 

1 

(SYY(l)) 

XY-Stress 

at 

center 

face 

1 

(SZZ(l)) 

ZZ-Stress 

at 

center 

face 

1 

(SYZ(l)) 

XY-Stress 

at 

center 

face 

1 

(SZX(l)) 

YZ-Stress 

at 

center 

face 

1 

(SXX(2)) 

XX-Stress 

at 

center 

face 

2 

(SXX(2)) 

YY-Stress 

at 

center 

face 

2 

(SYY(2)) 

ZZ-Stress 

at 

center 

face 

2 

(SZZ(2) ) 

XY-Stress 

at 

center 

face 

2 

(SYZ(2)) 

YZ-Stress 

at 

center 

face 

2 

(SZX(2)) 

ZX-Stress 

at 

center 

face 

2 

(SXX(3)) 

XX-Stress 

at 

center 

face 

3 

(SYY(3)) 

YY-Stress 

at 

center 

face 

3 

(SZZ(3)) 

ZZ-Stress 

at 

center 

face 

3 

(SXY(3)) 

YZ-Stress 

at 

center 

face 

3 

(SZX(3)) 

ZX-Stress 

at 

center 

face 

3 

(SZX(3)) 

ZX-Stress 

at 

center 

face 

3 

(SXX(4)) 

XX-Stress 

at 

center 

face 

4 

(SYY(4)) 

YY-Stress 

at 

center 

face 

4 

(SZZ(4)) 

ZZ-Stress 

at 

center 

face 

4 

(SXY(4)) 

XY-Stress 

at 

center 

face 

4 

(SYZ(4)) 

YZ-Stress 

at 

center 

face 

4 

(SZX(4)) 

ZX-Stress 

at 

center 

face 

4 

(SXX(5)) 

XX-Stress 

at 

center 

face 

5 

(SYY(5)) 

YY-Stress 

at 

center 

face 

5 

(SZZ(5)) 

ZZ-Stress 

at 

center 

face 

5 

(SXY(5)) 

XY-Stress 

at 

center 

face 

5 

(SYZ(5)) 

YZ-Stress 

at 

center 

face 

5 

(SZX(5)) 

ZX-Stress 

at 

center 

face 

5 

(SXX(6)) 

XX-Stress 

at 

center 

face 

6 

(SYY(6)) 

YY-Stress 

at 

center 

face 

6 

(SZZ(6)) 

ZZ-Stress 

at 

center 

face 

6 

(SXY(6)) 

XY-Stress 

at 

center 

face 

6 

(SYZ(6)) 

YZ-Stress 

at 

center 

face 

6 

(SZX(6)) 

ZX-Stress 

at 

center 

face 

6 

A-3 


9.  PIPE 


A.  TANGENT  (12) 


(  1) 

(PX(I) 

) 

X-Force 

at 

end 

I 

(  2) 

(VY(I) 

) 

Y-Shear 

at 

end 

I 

(  3) 

(VZ(I) 

) 

Z-Shear 

at 

end 

I 

(  4) 

(TX(I) 

) 

X -Torque 

at 

end 

I 

(  5) 

(MY(I) 

) 

Y -Moment 

at 

end 

I 

(  6) 

(MZ(I) 

) 

Z -Moment 

at 

end 

I 

(  7) 

(PX(J) 

) 

X-Force 

at 

end 

J 

(  8) 

(VY( J) 

) 

Y-Shear 

at 

end 

J 

(  9) 

(VZ(J) 

) 

Z-Shear 

at 

end 

J 

(10) 

(TX( J) 

) 

X-Torque 

at 

end 

J 

(ID 

(MY( J) 

) 

Y -Moment 

at 

end 

J 

(12) 

(MZ( J) 

) 

Z -Moment 

at 

end 

J 

(  1) 

(PX(I) 

) 

X-Force 

at 

end 

I 

(  2) 

(VY(I) 

) 

Y-Shear 

at 

end 

I 

(  3) 

(VZ(I) 

) 

Z-Shear 

at 

end 

I 

(  4) 

(TX(I) 

) 

X-Torque 

at 

end 

I 

(  5) 

(MY(I) 

) 

Y -Moment 

at 

end 

I 

(  6) 

(MZ(I) 

) 

Z-Moment 

at 

end 

I 

(  7) 

(PX(C) 

) 

X-Force 

at 

center 

of 

arc 

(  8) 

(VY(C) 

) 

Y-Shear 

at 

center 

of 

arc 

(  9) 

(VZ(C) 

) 

Z-Shear 

at 

center 

of 

arc 

(10) 

(TX(C) 

) 

X-Torque 

at 

center 

of 

arc 

(11) 

(MY(C) 

) 

Y -Moment 

at 

center 

of 

arc 

(12) 

(MZ(C) 

) 

Z -Moment 

at 

center 

of 

arc 

(13) 

(PX(J) 

) 

X-Force 

at 

end 

J 

(14) 

(VY(J) 

) 

Y-Shear 

at 

end 

J 

(15) 

(VZ(J) 

) 

Z-Shear 

at 

end 

J 

(16) 

(TX( J) 

) 

X-Torque 

at 

end 

J 

(17) 

(MY( J) 

) 

Y -Moment 

at 

end 

J 

(18) 

(MZ( J) 

) 

Z-Moment 

at 

end 

J 

B.  BEND 


(18) 


